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Renaissance Data Systems 
P. O. Box 421- V 
Georgetown, CT 06829 
(212) 864-3078 


Books on APL and J and other curiosities 
of merit! 


Renaissance Data Systems announces a change in its 
mailing address. Please note that the telephone number 
remains the same. 


If you would like a copy of our latest catalog, please send us 
a self-addressed legal sized envelope with one first class 
stamp (if in the U.S.). 


Included are such titles as: APL is EASY, APL - An 
Interactive Approach, APL2 at a Glance, APL - the 
Language and its Actuarial Applications, APL as a Tool 
of Thought proceedings, I-APL publications and 
software, Boolean Functions and Techniques, The 
FinnAPL Idiom List, The Toronto APL Toolkit, 
Mathematical Experiments on the Computer, Probability 
in APL, APL - Stat: Do it yourself guide to computational 
statistics, A Source Book in APL - Approximately 80 titles 
in all! 


We also carry J publications, including Programming in J, 
An Implementation in J (structure and source code), 
Arithmetic, and Calculus. 


Shareware and commercial APL interpreters and J2 
interpreters are available as well. 
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Editorial: Is There Anybody Out There? 


by Anthony Camacho 


Does anybody read the notices we send out inserted in Vector? 


APLers who organise meetings, write articles or edit magazines need to know 
what you think and what you want. After a while we will have commissioned all 
the articles we can think of that we really want to read and organised all the 
meetings with the speakers we most want to hear. When that time arrives we 
don't know what to do next unless you tell us. 


When we ask you and not one person suggests a speaker or subject (there was 
one offer to speak!) we get worried, What would you do? 


What Sort of Conference? 


One prominent member of the APL 96 committee has suggested to me that the 
traditional APL conference is dead and that what is needed is something much 
more organised than a tidy arrangement into streams of the papers that fall 
through the letterbox onto the doormat. "What we need is workshops.” 


Is it true that the number of papers submitted for review is steadily declining? Is 
it true that the number of delegates paid for by their employers is steadily 
declining? Is it true that the importance of APL skills is declining relative to those 
skills needed to control the operating system? Should APL conferences contain 
Windows workshops? (Once I thought I was going to spend the rest of my life 
learning a new word processor every six months; now I think I may have to learn 
a new operating system every year.) 


Who can provide the information needed to enable conference planners to do a 
good job? АП the conferences have been run by SigAPL so they ought to be able 
to do this. 


It would also be nice to know, for each conference, how many exhibitors there 
have been (and how much money they contributed), how many nationalities 
attended, how bookings were distributed between the earliest and the latest (this 
is really vital for people doing the financial planning) and a whole lot of other 
things which any intelligent person could suggest. If someone at SigAPL would 
do this we will be glad to print it in Vector. (Or to read it in Quote Quad.) 


If we want to buck the trend, we need to know what the trend has been. 


dyalog КЁ] 


The Definitive APL for Windows" 


LI 


elelsiaialal® 


Dyalog APL/W Version 7.1 includes a sophisticated built-in Grid object, with 
numeric, currency, and date fields as well as combo boxes and button cells. The 
Grid also provides an undo feature, cut and paste facilities (which let you move 
data quickly and easily between APL and your favourite spreadsheet), resizable 
rows, columns and titles, and drag/drop editing. Dyalog APL/W supports Visual 
Basic Custom Controls, ToolBar, StatusBar and TabBar objects, automatic Hints 
and Tips, Metafiles, MDI, 3-D Forms and Controls, a fully customizable Session, 
an ODBC interface, namespaces for encapsulation, and a host of other features; all 
designed.to make it easy to develop fast, responsive and 

attractive Windows applications. 


That is why Dyalog APL/W remains the professional D 
choice. For further details, contact Dyadic or your local : 
istril today. MIROSOF 
distributor today. A 
WINDOWS. 


Dyadic Systems Limited, Riverside View, Basing Road, Old Basing, © COMPATIBLE 
Basingstoke, Hampshire, RG24 7AL, United Kingdom. 
Tel:+44 1256 811125 Fax: +44 1256 811130 Email: sales@dyadic.com. 
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Quick Reference Diary 1995-96 


Date Venue Event 


July 28th — The George Fox Complex, APLS6: Designing the Future 
August 2nd 1996 Lancaster 


British APL Association meetings are normally held in the IEE, Savoy Place. 
Nearest tube outlets: Temple or Embankment. 


BCSNet: Free Registration Offer 
Open to Specialist Group members who join the BCS as new Affiliate Members 


If you join the BCS as an Affiliate Member before the end of October 1995, the 
registration fees for BCSNet Services will be waived. 

BCSNet Gold registration fee is £25* VAT. 

BCSNet Lite registration is £104 VAT. 


For more information contact BCSNet, 1 Sanford Street, Swindon, Wilts, 531 1НЈ 
or ring 01793-417426 or FAX 01793-480270 or email netadminGbcs.org.uk 


Dates for Future Issues of VECTOR 


Vol.12 Vol12 Vol.13 
No.3 No.4 No.1 
Copy date 1st Dec 95 1st March 96 24th May 96 
Ad booking 8th Dec 95 10th March 96 31st May 96 
Ad Copy 15th Dec 95 20th March 96 7th June 96 


Distribution January 96 April 96 July (at APL96) 
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APL96 
LANCASTER 


July 28th — August 2nd 1996 


Key Dates for Contributors 


This Vector should be mailed with a Cal! for Papers enclosed, but just in case we 
fail to make it on time, here are the dates which contributors will need to meet: 


Abstracts as soon as possible, and definitely by the end of November. 
Nothing formal is needed, just a quick Email to me {Compuserve 100331,644) 
or Phil (benkardGaol.com) to let us know that something is on the way. 


31st December. Draft papers received by Programme Committee; we intend 
to have the last draft paper out into the review process by 15th January, with 
pressure on the reviewers to get everything back by 31st January. 


31st December. Workshop and panel proposals (again an informal Email is all 
we need} to one of the Programme chairs. 


9th February. Programme Committee meets; selected abstracts ready for the 
printer by 15th Feb to begin work on the Invitation 


21st March. Final copy of all papers to the Proceedings editor, who will aim to 
have this away to the printers by the end of April. 


If you have potential material which directly supports the major themes of the 
workshops, please let us know; there is no need for this to go through the 
traditional ‘double-blind’ reviewing process (unless you want it to), but we 
would like to have something sufficiently polished to include in the Proceedings. 


For those who have lost the leaflet that was enclosed with Vector 12.1, the major 
themes will be: Effective use of Nested Arrays; Designing for Windows 95; 
Communication among Co-operating Systems; Algorithms and other Tools of 
the Trade. As always, good, relevant papers are always welcome on any topic 
under the sun — we are budgeting on 24 slots, but this is not an upper limit! 
More is better, and more variety is better still. 


Adrian Smith 
APL96 Programme Co-Chair 
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CORRESPONDENCE 


My Computers Are APL Machines 


160 — 255 Снн) 


From: ВШ Chang Recd at APL95 
T x o£ o¥ 4 y Proposed Common APL Coding 
v А 1 1 id 2 M r 1 сап view, enter, edit, print APL from just 
РХ about any application, including the command 
Ц à line; the first version of this note was posted to 
7 * 1 ^ i € 4 v complangaplin March 1994. 
41. OOF To 
e Boe [o $ v w Нее my (revised) common APL coding. Гуе 
> f ¢ 4 © 0 + F implemented it on the Mac and Sun (xterm). 
à á А Z à a x s Туе found that WWW’s HTML file transfer 
г e ё € Î í % Y protocol passes 8-bit characters (as opposed to 
» й ò á 8 8 8 æ &#NNN) just fine, in accordance with 
# ù ú û à z ~ T documentation. You should see a table of ISO 
Latin (8859-1) characters. Since I already have 
û i «< f « «  » АРІ fonts installed on my machines, I see APL 
^&gg «H8 @ characters, both through news (rn) and the web 
doe rc (mosaic and also lynx, a portable character- 
Е з > ә # і & _ based web browser that works great). It works 
“yt for e-mail too, using the MIME protocol. 
128—159 Chad’) 
Rationale 


In the October "94 issue of Vector (11.2 р. 105) Adrian Smith proposed a common 
APL encoding for Microsoft Windows that includes (most) lower-case national 
characters but not upper-case. My major critique of this encoding (compatible 
with APL*PLUS) is that it uses the range 128-159 indispensably. While this is 
fine in Microsoft Windows or the Macintosh (which can even use most of 0-31), 
many network gateways and modem / terminal programs either treat 128-159 as 
0-31 control characters or mis-handle them in other ways. For example, xterm 
(X windows) refuses to display 0-31 and 128-159. The well-established ISO 
Latin1 (8859-1) encoding avoids 128-159, Almost all modern computers already 
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support this 8-bit standard, and mail/news/ web have come very close as well. 
As a result, we can expect 8-bit characters (at least 160-255) to be freely usable 
and transmittable. 1 believe APL can hitch a ride with very little effort. 1 have 
created a new encoding that places "essential" APL glyphs outside the "dead 
characters" 128-159, as much as possible. 


Software Availability 
Ican provide the following on request: 


+ PC/Windows fonts. (DDE keyboard that works in all applications is in 
progress.) 


Macintosh screen font and custom keyboard mapping (drag-and-drop into 
system) that work in all applications. 


The popular shareware terminal program Zterm 0.9, patched to use APL font. 


Instructions for making APL the system font (instead of Monaco), so the Mac 
becomes an “APL machine". 


Sun (X windows) screen font (BDF format) and keyboard mapping (xterm 
resources) that can be used in all xterm-based programs, such as the command 
Shell, unix utilities, and text editors. 


Simple installation instructions. 


A re-encoding of Adrian Smith's APL2741 (Type 3) PostScript font. 


Printing instructions and utilities. 
(Some are not yet finished but are trivial.) 


The custom keyboard is bona fide APL — Caps Lock toggles between "unified" 
and "standard" APL keyboards. The font originated in 1987 as a hand-edited, 
highly optimised version of Courier/ APL*PLUS that has the same bitmap size as 
the Mac's system font (Monaco 9). Over the years it has been revised and re- 
encoded for APL2 (RS/6000), VAXAPL, and APL.68000. I'm also working on 
finishing a larger bitmap font and possibly a TrueType font, as time allows. 


Please comment. Can someone port this to the PC? 


Bill Chang (wchang@acm.org) 
Tel: +1 (616) 367-8866 
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Some Comments 

by Adrian Smith 

I very much want to be able to cut and paste directly between my APL session 
(or editor), my word processor, and my Web browser. I can use Bill's suggested 
layout in Dyalog/W with no problems, but in APL*PLUS III I am stuffed, as 
there is no way of getting at their mapping table. Does this matter? Readers 
please advise! 


A couple of small quibbles — Bill has put « in the ‘bad’ area along with Ө and #. 
Tcan live without I and ¥ (would you ever send someone a locked function over 
the net?!) and possibly even [ (sorry, SigAPL) but I make heavy use of 9 and 
increasing use of є. Some kind of straw poll looks a good idea before we nail this 
one down for good by offering the fonts and .DOTs on an APL Web page. 


With these provisos, I would be delighted to offer the APL2741 TrueType design 
to the agreed encoding, and also to provide a bitmap session font (and 
appropriate web.dot) for Dyalog APL. 


Causeway 


Support and maintenance of Causeway installations 


You already have the software ... 
... but would you bet your business on it? 


We Do! 


If you want security, training, and a shoulder to cry on 


if it all goes wrong ... you should call us ... 
Now! 
Causeway Graphical Systems Led 


S The Maltings Tel: +44 (0)1653 696760 
Castlegare, Fax: +44 (0)1653 697719 
Malton, North Yorks 

YOI7 ODP 100265.1564@compuserve.com 
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British APL Association News 
from Sylvia Camacho 


Members may be interested to know what the BAA is doing about APL 96. As it 
is to be held in Lancaster one would expect the local APL organisation to be 
involved. Since the very first international APL conference these conferences 
have always been controlled by and usually financed by SigAPL, which is a 
special interest group of the Association for Computing Machinery (usually 
known as the ACM), the United States equivalent of the British Computer 
Society. 


Below is the text of the reply from our chairman, sent on 22 September to Bob 
Brown the SigAPL conference co-ordinator. 


Thank you for your email requesting commitment from the British APL Association to 
APL96 at Lancaster. 


Several members of the BAA returned from АРІ95 with the understanding that a 
budget for the proposed APL96 would be received in early July. The BAA committee 
met on 7th July, but no budget had materialised and so no decision could be made, as 
(like ACM) the BCS does not commit funds without a formal budget and proposal. 


We did, however, give our backing to the programme proposed by Adrian Smith and 
resolved to give financial backing to a conference based on that programme as soon as 
SIGAPL had produced the required budget. 


As you may be aware, a year ago the BAA declined an invitation to organise APL96. 
Those of our members who have the necessary skills are heavily committed elsewhere 
at this time and so we cannot take the lead in this project. 


Financially, we need a budget that breaks even with a relatively small number of 
delegates, as we believe that a conference in the UK with such a short lead time is 
unlikely to attract large numbers. Ifa suitable budget is produced, and Adrian is still 
willing to organise the programme, we will make BAA funds available as soon as 
SIGAPL does the same. 


Alan Mayer 
Chairman, British APL Association. 
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News from Sustaining Members 
Compiled by Gill Smith 


Insight Systems 


Insight Systems is proud to present two new Client/Server products, the SQAPL 
ODBC Server and APL Pipes, which make it possible to integrate application 
servers written in APL with other development tools, using industry standard 
interfaces and communications products: 


+ The SQAPL ODBC Server allows your APL application or database system to 
be accessed as an ODBC data source. You can present your APL-based data or 
the results of analyses as a set of relational tables which can be used as input 
into other reporting packages. It is also possible to use application generators 
such as Borland Delphi or Microsoft Access to create front-ends to APL 
applications via ODBC. The APL application must run under Unix or 
Windows NT, but becomes accessible from almost any end-user application or 
development tool on virtually any platform. 


With APL Pipes, you can write distributed APL applications without the 
overhead of a relational interface such as ODBC or SequeLink. APL Pipes is a 
TCP/IP based system which is similar to Shared Variables, except that it 
supports sharing of data and function calls between APL2, APL*PLUS П or Ш, 
Dyalog APL and SHARP APL systems running under Windows, Windows 
NT, 05/2 and Unix. Before the end of the year we plan to add Microsoft 
Visual Basic and Borland Delphi to the list of APL Pipes clients. This will mean 
that developers in these environments will be able to effectively share 
variables and make function calls to almost any APL application, whether this 
application is running on the same machine or accessible over the network. 


Until the end of the year, Insight Systems will be very busy with the completion 
of the Windows version of the KPS System for Adaytum KPS Software Ltd. KPS 
is a multi-dimensional planning system written in APL, and is currently the UK. 
market leader for budgeting software. We are looking forward to developing the 
Client/Server versions of KPS scheduled to follow the Windows product, which 
will allow us to practise our Client/Server preaching and enable APL-based 
solutions to compete head on with mainstream products in a highly competitive 
market. 
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Dyadic Systems Limited 

Announcing Dyalog APL/W Version 8 

What is Version 8? 

Dyalog APL/W Version 8 is a completely new implementation of Dyalog 
APL/W that is designed specifically for Windows 95 and for Windows NT 
Version 3.51. It is fully compatible with the current release of Version 7 (7.1 
Release 2), but contains many enhancements to support the Windows 95 user 
interface and other features. 


Version 8 is based on Microsoft WIN32 in place of the equivalent Watcom 
software used for Version 7. WIN32 provides improved memory management 
under NT and is a key requirement for the new Windows 95 logo for which the 
product is intended to qualify. 


Version 8 includes support for most of the new GUI objects and features 
provided by Windows 95 and Windows NT 3.51. These include ListView, 
ProgressBar, PropertySheet, RichEdit, Spinner, TrackBar and TreeView objects. All 
new and existing GUI controls have also been enhanced to support the drag- 
drop feature of Windows 95. This permits the user to drag-drop file icons into a 
Dyalog APL/W application. It is intended that the new version of Dyalog 
APL/W will also support OLE2 and .OCX custom controls. 


Availability 

Dyadic expects to be in a position to ship Version 8 by the end of 1995. However, 
an early release will be available by the end of September that includes much (if 
not all) of the additional functionality. Customers may obtain this release by 
enrolling in the Version 8 Preview Program (see below) 


What is the policy for 3.1 and 95? 

For the forseeable future, it is likely that many customers will continue to use 
Windows 3.1. Others will migrate immediately to Windows 95. To cater for both 
requirements, Dyadic intends to maintain two separate versions of Dyalog 
APL/W for as long as the requirement to do so exists. 


Version 7 will continue to be maintained for Windows 3.1 and OS/2 users. 
Version 7 is also supported under Windows 95, but it is NOT supported under 
NT. Customers who need to run applications under Windows 3.1 or OS/2, 
should remain with Version 7. Applications based upon Version 7 will also run 
under Windows 95, but will not be able to take advantages of the new Windows 
95 features. 


12 


VECTOR Vol.12 No.2 


Version 8 is Dyadic's product for Windows 95 and NT 3.51 and includes support 
for many of the special facilities provided by these systems. It is not supported 
under Windows 3.1. Version 8 is intended to be the state-of-the-art APL for 
Windows and will be the basis for ongoing developments in the future. 


Prices 


Version 7 and Version 8 will be marketed and supported as separate products. 
Version 8 is not provided as an upgrade from Version 7. However, there will be a 
reduced price for customers purchasing both Versions and for customers who 
already have Version 7 and purchase Version 8 in addition. Please contact 
Dyadic or your local Dyalog APL distributor for details. 


Version 8 Preview Program 


The Version 8 preview program is intended to allow customers to begin to 
develop Windows 95 (and NT 3.51) applications immediately, without having to 
wait for the final release of the Version 8 product. It also provides a basis for 
customers to influence development. This program is open only to existing 
Version 7.1 users and to any customers who purchase Version 7.1 prior to the 
final Version 8 release until further notice. The Version 8 Preview Program is 
chargeable but the price includes a copy of the final release which will be 
supplied as soon as it becomes available. Enrolling in the Preview Program is the 
cheapest way for existing Version 7.1 users to obtain Version 8, Contact Dyadic 
or your local Dyalog APL distributor for details. 


During the Preview Program, Dyadic will supply a stream of test releases as new 
functionality is added and as bugs in the existing code are fixed. In addition to 
bug reports, subscribers may submit suggestions and requests for enhancement 
to the new Windows 95 and NT features as they are developed. Dyadic does not 
guarantee to incorporate all of the requests in the final version, but will 
endeavour to meet popular demands. 


Summary of New GUI Features 


Version 8 provides 10 new objects, over 30 new properties and over 40 new 
events. These new features are summarised below. 


+ The PropertySheet object is a dialog box that allows the user to view and edit 
the properties of an item. A Standard PropertySheet contains one or more 
overlapping child windows represented by PropertyPage objects, each 
containing controls for setting a group of related properties. Each 
PropertyPage has a tab that the user can select to bring the page to the 
foreground of the PropertySheet. The Display Properties dialog (obtained by 
clicking the right mouse button on the desktop and choosing Properties) is an 
example of a standard PropertySheet object. A Wizard PropertySheet consists 
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of a sequence of dialog boxes that guide the user through the steps of an 
operation. in a Wizard PropertySheet, the PropertyPages do not have tabs, 
and only one page is visible at a time. Also, instead of having Ok and Apply 
Now buttons, a Wizard PropertySheet has a Back button, a Next or Finish 
button, and a Cancel button. Windows 95 Setup is an example of an 
application that uses a Wizard PropertySheet. 


The ListView object displays items as a set of icons, in a list format or a report 
format. ‘My Computer’ is an example of an application based upon this object. 


The TreeView control is a window that displays a hierarchical list of items, 
such as the headings in a document, the entries in an index, or the files and 
directories on a disk. Each item consists of a label and an optional bitmapped 
image, and each item can have a list of subitems associated with it. By clicking 
an item, the user can expand and collapse the associated list of subitems. 
“Windows 95 Explorer’ is an example of an application based upon this 
control. 


The RichEdit control is a multi-line text editor that provides extensive word- 
processing facilities including variable text and paragraph alignment, fonts 
and colours. This control allows Dyalog APL/W applications to exchange 
complete documents with Microsoft Word and other word-processors. Data 
may be imported and exported using -RTF files or via the clipboard. The object 
also includes direct printer support. 


The ProgressBar is a standard output control that indicates the progress of a 
time-consuming operation. 


A Trackbar control is a window that contains a slider and optional tick marks. 
When the user moves the slider, using either the mouse or the direction keys, 
the trackbar generates events to indicate the change. 


The UpDown control is a pair of arrow buttons that the user can click to 
increment or decrement a value, such as a scroll position or a number 
displayed in a companion control. 


The Spinner object is a special Dyalog APL object that combines an UpDown 
with an Edit. The Spinner permits the user to select from a range of numerical 
values or text strings. 


The Grid object has been extended to allow Spinner and Trackbar objects to be 
used to edit cell values. 


The List object has been extended. It has a new MultiColumn property that 
allows you to display items in several columns. In addition, the List supports 
the facility for the user to drag-drop items within it. 


Most controls support the new AcceptFiles property. If this is enabled, an 
object will generate a DropFiles event when the user drag-drops one or more 
file icons onto it from another application, such as ‘My Computer’ and 
“Windows Explorer’. 
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Causeway Graphical Systems Ltd 


Our main item of news is a very simple one: around 3 weeks ago, SAP started 
cutting our software onto CD to issue with their R/3 3.0a release. From SAP's 
point of view, 3.0a is a strategic release and is, to the best of our knowledge, the 
first truly distributed business management system to hit the world market. 
Perhaps the most encouraging thing about this is that a major player in the world 
software market has been willing to rely on APL for a critical configuration tool. 
This has to be a good thing for the future of our language. 


This leads naturally to the second item: Causeway — the Next Generation (working 
title: CTNG). For the past two weeks, we have had time to relax and work on the 
future of Causeway rather than hitting the impossible deadlines set by our 
clients. Since we first released the Causeway utilities in November 1993 (at the 
British APL Association's BAA-GUI workshop) we have seen a massive increase 
in the functionality of the Dyalog interpreter. In particular it can now handle 
resize and reposition (via the ATTACH property) in the way that we need it, and 
it has Namespaces. Anyone who has looked through the Gui. xx functions will 
have noticed that a significant slab of code in Gui, make is no longer necessary, 
and may have wondered what all the complex use of ПОК and DSHADOW is 
doinginGui exec. 


Rather than attempting a piecemeal upgrade, we decided to take the bull by the 
horns and do a top-down rewrite. It was important to preserve the dialogue-box 
structure, and it is clearly essential that CTNG behaves in exactly the same way 
ав Causeway, to the calling application. However the internal structure of the 
class table will be very different (and will support managed replication of 
functions between classes), and the implementation of localisation and event- 
handling has been built entirely around the namespace model. 


The main visible difference is a huge reduction in the lines of APL code needed 
to support form-level localisation of variables, and local functions on objects. The 
entire set of Gui. xx functions written so far can be listed on 3 sides of A4, and 
most of these will never appear outside & . Gu i. In particular, all the Causeway 
‘modules’ (such as the Rain PostScript interpreter) simply become namespaces, 
and the code in them can be executed much more simply with statements such as 
ps.Fuli$creenreplacingps Gui do'FullScreen' in application code. 


We are on course to have the first release of CTNG included with Dyalog 8; 
however we also intend to support it under Dyalog 7.1.2 for users who prefer to 
stay with Windows 3.11 for the moment. 
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APL for the Apple Macinto 


$һ 


Since 1985 MicroAPL has pioneered the use of APL in graphical environments. 
Our latest version of APL.68000 Level Il for the Apple Macintosh is now available, 


offering dramatically enhanced GUI programming facilities. 


MOVIES 
PPL 68828 Level 11 Copyright (С) 1990-1994 MicroAPL Ltd. in 


IS Size-3601K, Version 3.04 — met 
Weeds s Н. (SCIES Movies from APL ЫШ 


Justo 
99 POVIES 
Denotiovie 


MOVIES: DemoMovie 


(01  DemoMovie;FM;z;FILE 
UIJ — FMeNextForw 

{21  FILE«'Micro&PL:NOUIES:QUICKTIME LOGO MOVIE’ 
СЭ] zeFM Dei ‘New! 'Form'Ccaplion' ‘Movies from 
a) А 

[S] а Now create a movie. 


16] ZEN," anol Dei "Mew "Movie! shore! 13007 
зол QuickTime 
19] лава a quit button 

191  z«(FM,'.Quit)Dwi ‘Nev! ‘Button’ 'uhere 13 1 — | 


[101 я Force the user to exit the form 


261 OWE “1 


APL.68000 Level Il for the Macintosh includes the folowing features: 


. Runs on all models of Apple Macintosh 

. Native version for the Power Macintosh 

. Conforms ciosely to the APL2 specification 

Е Uses the standard Mac interface 

L Object-based GUI programming via OHI 

LI Full event handling via APL callbacks 

* Free runtime version with application packager 


Accelerated for 


London, SE1 6LN, UK 


Voice: 0171 922 8866 
Fax: 0171 928 1006 
Applelink: microapl 


Internet: microapl @microapi.demon.co.uk 


Power Macintosh 


File Edit Interrupt Terminal Special Window D Ww) 


со 


MieroAPl. 


APL 68000 Level 


MicroAPL Limited South Bank Technopark, 90 London Road, 


Bloomsbury Software 


report some very recent news on the sale of the APL*PLUS 
microproducts by Manugistics and future development plans 


Huge New Investment in APL*PLUS 


Manugistics Inc. have announced the sale of all their rights and interests 
in APL PC software, including APL*PLUS IL, APL*PLUS III for 
Windows, and APL*PLUS H for UNIX, to the Financial Reporting 
Software & Consulting Company (FRS) of Georgia, USA. The sale 
includes the transfer of the existing development and support personnel 
to FRS 


Development team to double in size 


Eric Baelen, Chief Executive of FRS, says that his company has been 
using APL*PLUS programming environments for their financial 
reporting products for many years. He adds that they intend to make a 
significant investment in the continued development of the APL*PLUS 
products, starting with a doubling in size of the current development 
team 


Informed sources at Manugistics, whilst saddened at the departure of the 
APL*PLUS microproduct range, view the move as a very positive boost 
for the product itself 


Windows *95 support soon 
It is expected that FRS will move rapidly to produce a Windows '95 
compatible release for APL*PLUS III for Windows for release first 
quarter 1996, and also a significantly enhanced product release 
supporting OLE2 and OCX controls by the middle of 1996 


/Over 


Price increase soon 


However FRS will be seeking to recoup their expenditure with 
significant product price rises in the near future and we suggest that 
anyone with plans to purchase any APL*PLUS products would be 
advised to place orders with us before the end of October to ensure 
availability of current prices 


Mainframe APL*PLUS development stays at Manugistics 


Manugistics confirm that only their APL microproducts have been sold 
and that their range of mainframe APL interpreters will continue to be 
developed, supported and marketed by Manugistics themselves. 


Bloomsbury Software welcomes the news 


We at Bloomsbury Software welcome this news without reservation. 
The future of the APL*PLUS product had been under a cloud but having 
spoken with Eric Baelen we have renewed confidence that APL*PLUS 
Ш Windows will continue to be developed for the new operating system 
to give developers the advantages of APL*PLUS and Windows ‘95 


We look forward to working with FRS and to continuing to support UK 
users of the APL*PLUS products through supply and support of the 
software, training, application support and consultancy. 


The Bloomsbury Software Co. Ltd. 
- formerly Cocking & Drury (Software) Ltd - 
3-6 Alfred Place, London, WCIE 7EB 


Phone 0171 436 9481 Fax 0171 436 0524 
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Editorial 
Raising Standards — or Flying Kites? 


If you are one of those people who believe that reading newspapers, like tarot 
cards, confers insight into reality, you cannot fail to have noticed that 
something’s the matter with our universities this year. They're being swamped 
by a tide of mediocrity. 


The trouble started with the A-level examination results. Good grades were 
achieved by too high a proportion of school-leavers, which entitled an excessive 
number to apply for university admission. Since the human genome doesn’t 
change that much in the space of a few years, this could only come about, 
according to some, by a general relaxation of standards, 


There is another explanation, which the government has been surprisingly slow 
to latch onto. It is that the painful re-engineering of public education and welfare 
over the last few years has at last borne fruit in a cohort of better educated 
students, additionally motivated by a shot of genuine fear for their future 
prospects, Students, indeed, who are more able to jump the traditional hurdles 
designed to filter the candidates for tertiary education. 


Perhaps the government doesn’t believe this (a cynical standpoint, in view of 
their efforts) but suspects the hurdles have been lowered out of self-preservation 
by unscrupulous academics. But whom exactly do the hurdles filter out, and 
why? 


The traditional concept of matriculation and its successor, the A-level 
examination, is founded on the premise that only a proportion (from a tenth to a 
half) of people who put themselves forward for higher education ought to be 
accepted for it. Otherwise “standards fall”. 


There are two views about why standards should fall. The first is that resources 
ate too limited to allow everybody who wants tertiary education to get it, which 
invites calls for an expansion of places on degree courses. The second is that you 
don’t want to over-educate the duffers. This invites people to believe that there is 
some rational criterion being applied to the acceptance and rejection of 
candidates, a belief which doesn’t bear close examination. 


Of course it is open to the more sought-after universities (whatever they imagine 
they are doing) to select those candidates who can tolerate old-fashioned 
teaching of patchy quality — in effect teach themselves. This frees the tutors to 
concentrate on research, consultancy, authorship, cutting a dash in public and 
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other more rewarding activities than the prime purpose of paying them out of 
the public purse — to train the nation's youth. Other establishments ape their 
betters, hoping to avoid the slur that their only role in life is to scavenge the 
rejects. 


But do rational rejection criteria really exist, let alone are they being applied? A 
longitudinal study of student performance at my old place uncovered no 
significant correlation (in fact zero) between students’ grades on admission and 
their eventual class of degree. Granted there are many caveats surrounding this 
Kind of study, it did raise an alarming question. When we selected candidates on 
the basis of their grades, just what did we imagine we were doing? 


Are there no predictors, then, of a candidate's eventual performance? No, the 
study didn't conclude that. Indeed certain syndromes were associated with what 
might be called (again with a legion of caveats): "losers". These are the hard- 
do'ers and the satisficers, who are capable of putting their tutors, counsellors, 
colleagues — indeed the whole institution — to enormous trouble on their 
behalf, at small cost to themselves. The sort of person you don't get to help you 
carry a ladder, or you end up carrying it yourself, plus him or her dragging 
along on the other end. Purely on cost-saving grounds, these students amply 
repay weeding out at the earliest stage. 


But that accounts for no more than 5*6 of the intake. For the rest, it's not so much 
a creaming-off process, more a thinning out, like seedling in a tray. Gardeners 
among us will know the rational basis for this. It is done like decimating a 
Roman legion, because its purpose is not to select the "best" (most of the 
seedlings look alike) but to relieve overcrowding. 


What is the evidence for overcrowding in today’s educational market? It all 
depends on the course. Romantic ones like Veterinary Science, or lucrative ones 
like Law, are heavily oversubscribed. For Mathematics, Physics, Chemistry and 
Engineering, all the signs point to massive overcapacity. In the Tyne-Tees area 
alone there are now five universities, whereas purely to serve the region's needs 
for people educated to degree level you'd be hard-put to make a case for one. 
Higher education has simply replaced coal and ships as a commodity for export, 
not to mention occupying the unemployed masses. 


Yet we read [Independent, Thursday 31st August, p4] that “7 thousand students 
are admitted to University each year through foundation courses which do not 
always require them to have A-levels. Mrs Gillian Shepherd said applicants must 
be told they had to meet rigorous standards to ... embark on a 3-year degree”. 
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This is revealing of an attitude to University entrance devoid of rational basis. To 
the youth of today, Mathematics, Physics, Chemistry and Engineering sound like 
a lot of hard work for pretty small rewards. Is it any wonder that courses are 
under-subscribed and departments face closure? It’s a matter of public concern 
(or scandal, depending on how you look at it) that a high proportion of 
engineering degree students need remedial Mathematics at the how-to-balance- 
brackets level. But it's only being realistic to admit such ill-trained candidates, in 
circumstances that have arisen over the years through public neglect or even 
deliberate policy. When there's been a frost in the greenhouse, you don't thin out 
the seed-tray. You coddle whatever germinates. 


Some Articles in the Current Vector 


GENERAL ARTICLES 

The Ball Clock Problem Roger Hui 

Is APL a Team Sport? Douglas Bohrer 

Multiprecision Arithmetic (Part 2) John Sullivan 
The Random Vector 

Computing Clopper-Pearson Confidence 

Limits by the Illinois Method Dietrich Trenkler 

The Incomplete Elliptic Integrals and APL Joseph De Kerf 
TECHNICAL SECTION 

At Work and Play with J: Gene McDonnell 

The Bauer-Mengelberg Problem 

Elegant Programming Chris Burke 

A Fractal Verb in J Richard Oates 

J Locales Richard Oates 
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J-ottings 7 


by Norman Thomson 


With the release of the J Release 2 interpreter as Freeware (hereafter referred to 
as JFW which conveniently stands for either J FreeWare or J For Windows), J has 
reached a crossroads in more ways than one. In the first place, the amateur / 
education / non-commercial J user who it is assumed might be reading, this 
column, must decide whether to opt for J7 (that is the last release of DOS J) or for 
JFW. If he/she cannot, or does not, run Microsoft Windows, then the choice is of 
course made. Where the choice is available, my vote is for JFW, primarily 
because it includes control structures. However, DOS J fanciers should not stop 
reading; a template for the equivalent form of a simple loop in a dyadic defined 
verb is the following: 


E 37 
fne.3 :0 fn=.0 :0 

1 

“initialisation “inttialisation 

vhile. condition do. body lab) $.-.»( condition )(end;$. 
end. result body C $.2.1ab 

) end) result 


) 


As an example, here are two equivalent defined verbs which compute continued 
fractions iteratively, e.g. 2 + 1/(2 + 1/(2 +... 100 times. The left argument is the 
number of iterations, the right argument is the start value for the fraction, so that 
for this example x. = 100 andy. = 2. 


El a7 
fn-.3 :0 fn=.0 :0 
"n tr 
rey. C 1=.0 сеу. D feet 
while. i«x. lab) $.-.»(x.»i». 1«1) (end;$. 
do. r=. (+%)r r=. CX). 
1=.1+1 $.-.1ab 
end. г end)r 


) ) 
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The JEW version shows what might be regarded as a “nice” layout for the loop 
— for those who like maximum compression, the loop body could be written 


while.i<x.do.r=,(+%)r(Ci=.i+i епа, 


— shades of APL one-liners, aided by the use of the verb С as statement 
separator. It is stressed that the above verb is written to demonstrate a loop — if 
the result alone were the goal, then the natural J way to write it would be to use 
the power adverb 


(+%)^:100{2) 


To test whether the benefit of control structures is paid for by performance 
overhead, an “empty loop" (monadic) was implemented under both interpreters. 


ч 37 
fnz.3 :0 
i=.1 
while. icy. 
$.=.>{у.>1=.1+1){епб;$, 

do. l=.i+1 $.-.lab 
end, end)i 


) d 


(Note: In J7, it is necessary to define the null dyadic option explicitly, but this is 
not the case in JEW.) 


The JFW version runs almost twice as fast as the J7 version, indicating that the 
control structures per se give a performance improvement. However, with some 
algorithms, J7 is faster than JFW — all that can be said with certainty is that the 
two interpreters are different. 


Apart from control structures, the other major difference between J7 and JFW is 
the definition of amend. This is illustrated using one of the simplest numerical 
algorithms which cannot be expressed without a loop — even in J! Choleski 
computes the “square root” of a square symmetric matrix, that is, for a given 
matrix M, the Choleski matrix C has the property that C'C = M where dash 
denotes transpose. 
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y 
[E] NB. transpose 
+n NB. inner product 
ipte NB. "square" of matrix, i.e. M'M 
є©зд=.[ ( ({з4) NB. inner prod of one column of matrix 
rde. (e: m) }.( NB. drop 1+1 elements from ith rov - 


NB. rd stands 


or "right of diagonal” 


choleskis.0 : 0 
i=. a Û zm, (5y.)50 NB. Initialization 

lab) $.*, >((#y.)=1=.1+1){$.;end NB. Start of loop 

(QC (Co. 1,1) (2-1 esq 2)(1 »: y. ))z NB. Diagonal element 

(CE rd y.)-1 rd sq z)Kt(z)(1 rd 1.$y.))2 NB. Assign to rt of diag 
slab NB. Back to start of loop 
end)z NB. Publish result 


NB. (i rd sq M) is products of col i of M with all cols to its right 


tmn-.3 381024253439 
Jum-.choleski tm 
3.16228 0.632456 1.26491 
0 2.14476 1.02576 
о 0 2.51949 
tm-:(tr ош) їр um 


JEW 

Instead of rd define 
each=.2.> 
rdiz.(>:@C }.,each)1.ee NB. rdi = "right of diag" scatter inds 
choleskiz.3 : 0 

1e. 1 [ z-.(Sy.)$0 NB. Initialization 

while. (#y.)>i=.i+1 do. NB. Loop through rovs 

z=. (%:(t{y.)-1 esq 2) (t=.<1,1) )2 NB. Amend diagonal element 
if. 1-:_1+#у. do. NB. Except for last one, 


amend 
z=. ({{(p{y.)-pisq z)%t{z) (р=.1 rdi z) }2 NB. elements right of diag 
end. 

end. 

) 


NB. (1.0)(1.0))2 is valid in J7 but NOT in JFW, which provides an 
opportunity for introducing the if. control word in the above. 
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Compare the third and fourth lines in J7 with the third and fifth lines in JFW, all 
of which update z using amend. The overall structure of amend (see J-ottings 6) 
is the same in both JFW and J7, namely 


(data) (indices) ) (data object) 


However in J7 "indices" are interpreted as “linear indices", whereas in JFW 
"indices" means a sequence of boxed scatter indices, thus making indices 
consistent across both { and } . 


The verbs rd and rdi share the same fork structure: 
CoD 2. +) (where v = { in J7, „each in JFW) 


The composition »:€L means "add one to the left argument", which in 
composition with >. (drop) identifies elements to the right of a matrix diagonal. 


Tn the fourth line of the J7 version, rd is used repeatedly with different matrices. 
In particular, within the indices component of amend, it is used with 1. to 
produce linear indices. In JFW гаї produces the corresponding scatter indices. 
Both t and p are assigned within indices components, and then used for selection 
in the corresponding data components, thereby demonstrating index 
consistency. 


The arrival of JFW raises broader issues concerning the educational value of J as 
an executable notation. In the days when APL was trying to make its way in the 
computing world, the smart way to condemn it without going to the trouble of 
finding out anything about it was to repeat the glib assertions: 


"Needs special equipment; 
"Costs a lot"; and 
“Hasn't any control structures". 


JFW has removed all of these objections, and it will be interesting to observe the 
excuses which the computer science fraternity will doubtless bring forward for 
ignoring J in its turn. 


The potential of APL as an executable notation for conveying mathematical and 
more general data-structuring ideas has long been proclaimed by small bands of 
enthusiasts around the world. Now that the above objections have been 
removed, JFW can provide even greater generality and power for achieving the 
above goals. Howard Peelle at the University of Massachusetts, has for the 
purpose of in-service teacher training produced admirable workshop material 
for Teaching Mathematics with J. Briefly the principle is to introduce a minimal 


24 


The Education Vector Vol.12 No.2 


amount of J, and use the combination of experimentation, discussion and 
reflection to focus attention on mathematical concepts. However, efforts of this 
sort in curriculum development need encouragement and funding to sustain 
enthusiasm. 


In many ways the ground is less fertile now than it was in the APL days. First, 
computer algebra systems such as Maple and Derive have achieved considerably 
greater popularity than executable notations. Unquestionably these have brought 
many conceptually simple but computationally heavy problems within reach of 
students, which is fine if results as opposed to methods are the primary goal. 
Secondly, the pendulum of educational fashion seems to be swinging away from 
knowledge towards concepts such as "Vocational Qualifications" and "Generic 
Skills". (Arguably data-structuring and data-manipulation are themselves 
generic skills — given the laxity of spelling these days it might even be possible 
to promote them as J-eneric skills!) 


Thirdly, in the UK, teaching by administrative overload seems to be taking over 
at all levels from teaching by inspiration. Consistent uniformity, and customer 
delight elicited by satisfaction surveys are of course to be applauded when mass- 
producing manufactured items, and such qualities are rightly recognised as 
indicators of excellence in this domain. However, it is disturbing to see the same 
indicators being accepted in an uncritical way in the education process, where 
excellence ought to mean something of a quite different sort. Cambridge 
University, for example, has long enjoyed a reputation for providing an 
environment in which budding excellence in mathematics can flower. However, 
in the words of the Professor of Applied Mathematics and Theoretical Physics 
there, there has been since 1990 a steady but quite rapid deterioration in the 
levels of preparedness of first year students which means that it is no longer 
possible to teach all the material. We must, I suppose, be grateful that the Greeks 
never hit on Total Quality Management, which if in place, would presumably 
have put paid at an early stage to the careers of Archimedes, Pythagoras and 
their like! 


Such matters heighten the significance of the crossroads at which J has now 
arrived. Who and what is J for, and what should its place be in the broader 
spectrum of science? — these are questions which are now ripe for debate. Please 
send your comments to the editor. 
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Two Numerical Algorithms in J 


by Antje Muller, Tineke van Woudenberg and Alistair Young 


Polynomial Interpolation 


Interpolation means the process of using a table of values of an independent 
variable x with corresponding function values f(x), to determine a value of f(1) at 
a non-tabular value of x. For many functions, a polynomial which passes through 
the tabulated points gives an adequate representation near the non-tabular value. 
Such a polynomial is then used to interpolate values of f(x) at other values of x 
without explicit calculation of its coefficient values. 


There is a variety of techniques for doing this, among them Lagrangian 
polynomials and Neville's algorithm. The technique used here is known as 
Newton’s Divided Differences, and consists of taking successive differences of 
f-values which are divided at each stage by the matching differences in the 
column of x-values. The values of the nth column of f-differences are divided by 
the differences of the corresponding x-values which are n apart, for example: 


x t ist 2nd эга 
diff diff diff 
o 4 
2 
1 6 0.2 
3 -0.15 
s 18 -1 
-4 
a 6 


In the 1st diff column 18-6 is divided by 5-1, and 6-18 by 8-5. In the 2nd diff 
column, 3-2 is divided by 5-0 and -4-3 by 8-1, and so on. 


The result of the verb divdifs is x laminated with the values at the head of the 
columns, viz. 

015 

4 20.2 0.15 


from which the original values of f could, if required, be reconstructed by 
reversing the arithmetic. 
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The J realisation which follows is in the style of a Pascal program, which is a 
natural way of expressing an iterative algorithm of this sort. Indeed this 
algorithm illustrates the ease with which transition is possible between J and 
more familiar programming languages. 


NB. Obtain leading values in divided difference columns 
divdifs =. 0 : 0 


f 2.0 С f=. as. y. [ o=.#x, NB. initialization 


loop) $.-.»(n > is.1e1) ( end;$. NB. for i = 1 to n-i do 
Je. 1-1 NB. start of inner loop 
loopi) $. =. »(n > ј=.]+1) ( endi;$. NB. for j = i to n-i do 
toc. GU Ја (RG -= (]-104х.) NB. alj1:= 

a tja NB. ((fCj1-f0j-21)/(xC]]-xC]-11) 
$.-.100pi NB. end of inner loop 
endi) f=. а NB. update f 

$.».1o0p NB. end of outer loop 
end) x. NB. result is x,:leading diff'nces. 


) 


Interpolation is performed by an algorithm which uses the leading divided 
differences. A merit of this technique is that as the algorithm progresses through 
its iterations the result vector grows by catenating interpolated values of 
successively higher polynomial degree as more points are brought into 
consideration from the right of the data. 


NB. Program to perform polynomial interpolations of increasing order 
Interpol =. 0 : 0 
" 


foc. iix. [Ox =. Ox. NB. f - function values, x - x-values 


iit =. x divdifs f NB. а = leading divided differences 
пх, #а Гі =, Of p=. O(a КВ. p ca(0), first estimate 
sol =. a NB. sol is result vec (see last NB.) 
loop) $. =. >(п>1 =.i+1} ( end;s. NB. for i = 1 to n-i do 
ESER з NB. start of Inner loop 
loopi) $. =. >(_1<j=.j-1) ( endi;$. NB. for j = 1-1 downto 0 do 
toc. t*(y. -](х) NB. t z= te(y. -х[]]) 
$. =. loopt NB. end of inner loop 
endi) p =. pet«(CD (а) NB. p:sp + twati] 
sol =. p 1 } sol NB. sol(i]:-p 
$.=. loop NB. end of outer loop 
end) sol NB. sol is successive polyn'l interpolns. 


) 
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X -.8 9 9.5 11 
f =, 2.079482 2.197225 2.251292 2.397895 
xf =. ху: 
xf interpol 9.2 

2.079442 2.22078 2.21924 2.21921 


The last line shows that the cubic interpolation is 2.21921, although this is little 
different from the value obtained by quadratic interpolation. 


It is not necessary that the x values be in ascending order. For example, using the 
first example given above the successive polynomial interpolations are 


4 10 11.2 13 


If the points are introduced in x-value order 0 8 1 5, the polynomial 
interpolations are 


4 8,75 8.5 13 


which reflects the fact that the straight line approximation will vary greatly 
according to when the point (5,18) is introduced. 


Choleski Decomposition ("Square Root") of a Matrix 
This is a method of factorisation of a symmetric matrix 
A=LL’ 


where L is a lower triangular matrix. If A is positive definite the elements of L 
will be real, otherwise some will be complex. The mathematical equations which 
govern the method are: 


ia 
ma s non] k+1,....,nik 
Er 
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These lead to the definition of some subsidiary verbs as follows: 


951) и «81 ) 


ssqbutis. +/8( (~ NB. 2 ssqbuti & 2 7 3 5 ylelde the 
mum of squares of à 2 3 and 5 
vb 2 з 4 yields 
sqrt(2), S/sqrt(2), w/sqrt(2) 
performs vb on first column 
of a matrix 
column yields 0 n 2n зп ...(n-1)n 
NB. vhere nesy. (y. is a matrix) 
NB. x. fillmatrix y. ylelds matrix of 
NB. same size as y. filled with x. 
' (нїарї y.) (column y.) ) (0 flitmatrix у.)' 

NB. 1 fillMi m yields zero 
matrix vith first column 
МВ. containing the result of stapi 


узе. осе) 
atapis. (vb. )R.|: 


columns. (+{.) 


fillmatrix ». $-$ 


"зм. tt 


As with interpol a Pascal program is a natural way of expressing the 
essentially iterative algorithm, whose development, given the above verbs is 
now straightforward: 


NB. PROGRAM cholemki 
ch =, 0:0 
ke. 0 C n=.#Me. 1 fIDDMi y. NB. initialisation 


100р1) $. =. »(n»ke.ke1)( endii$. 
te. X: KIRBY. (к ssqbuti (KOD) 
Me. t((ne2)eR) 2 M 

j= k 

loop2) $. =. »(m»je.je1) {end2;$. 
te. GURY. (HM) +/ e jM DX t 
Me. t Uj n) + кури 

в. =. loop? 

епаг) 

$. *. loopi 

endi)k 


) 


NB. if k »- n go to endi 
NB. t = sqrt(AKk - sum(Mksa2) 
кк: t 

. start of inner loop 

NB. if j>= n go to end2 


WB. ts=(AjK - sum MjRMks)/Mkk 
NB. Mjk:=t 

NB. end of inner loop 

NB. jen 


NB. end of outer loop 
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Two examples are: 


А 1х=.сһ A xe/ om [eX 
4 2 18 2 00 4 2 16 
217.5 1 40 217.5 
14 _5 83 7.55 14 _5 83 
а cha 
146 1 0 0 
403 4 ojs o 
632 6 0j5.25 0]2.53722 


The authors of this article ате post-graduate mathematics students at the 
Universities of Kaiserslautern, Eindhoven and Strathclyde respectively, and 
were introduced to ] version 7, and to each other, in the course of ECMI 
(European Consortium for Mathematics in Industry) exchanges at the 
University of Strathclyde. 
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Technical Note on Matrix Decomposition 
by Norman Thomson 


J contains, as a primitive, 1281:0, the so-called QR decomposition of a square 
matrix. This note describes how to obtain this in APL, and also how to obtain the 
LU decomposition where L and U stand for Lower and Upper triangular 
matrices, in the latter case with 1s on the diagonal. L and U can if required be 
merged compactly into a matrix of the same shape as the orginal matrix by 
overlaying all those cells which, by definition, must be either 0 or 1. 


QR decomposition is also known as Gram-Schmidt orthogonalisation. The 
comments on the functions should make the procedures clear, particularly if 
used alongside the sort of algebraic descriptions of the algorithms which can be 
found in many texts on Numerical Analysis such as "Matrix Computations” by 
Golub and Van Loan, or the “Numerical Recipes" series. Perhaps surprisingly, it is 
in the nature of the algorithms that they are not particularly helped by the 
availability of the array operations of APL. 


¥ OLeQRORITiUiIiN 

[1] a R is square. Re>+.x/Z, cols of +2 are orthonormal, 221 is upper ar 

E2)  NetoR о Ze20c(H.0)pI-0 è initialize, Z is Built col by col 

Сај LirUFI-161)/0 ¢ T-RLI ^ begin loop. WORM is sqrt «/x«2 

[41 ZeiE 6 CLEAN 2,” (cU+HORM 0), НЕТ, NORM Ues/ (1, -TeTe .Ix+2)x[2]7,%2 ¢ +L1 
E 


Y Z-BORM R 
[1] — Ze(+/Re2) 00.5 » Euclidean norm 
9 


¥ ZL CLEAN R 
[1) — ZeRxLs|R a sets small numbers to zero 
' 


9 Т+Ш/ОВ:Т;НТ;Е:0 
[1J а #2 15 lover Аг, 252 15 upper ûr with diag is. Re>+.x/LU 1 
[22 28 о UsiWetpR о I«K«0 initialize for 1oopL[3-6) 
Сај Lire UPsIMI01)/L3 o Tell o +(I=1)/L2 first row 15 special case 
Гај KeZL,Iii1-1]8.»2(i1-1;7] K is for use in [51 
[51 L2:2UE; T] (RE 15 f -X) 20151] set Ith row to rt of diag 
Сеј {(1;1+11+#{71:1+1]-17;11)+.%2[11;1+1] о Li a diagebelov of (I+1)th col 
C7] £8: Ze (Ge, 20) &Z) (009.0) 4 (Un eU) x2) omit [7] for compact form 
+ 
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и 
8218 
1518 
103248 
9727 
QR M 
0.5101 70.2921 70.5667 0.5774 15.68 7.268 2.997 8.507 
0.06376 0.7761 0.2556 0.5774 о 5.846 0.8929 7.602 
0.6376 70.2796 0.7179 0 [] ° 0.%723 0.8879 
0.5738 0,485 70.3212 70.5774 0 o 0 2.309 
LUM 
во 0 0 1 0.25 0.125 0.375 
1 4.75 0 0 01 0.1842 1.605 
10 0.5 0.6579 0 00 1 70.84 
9 4.75 0 тъ 00 ° 1 


Vector Back Numbers 


Back numbers of Vector are available from: 


British APL Association, 
c/o Gill Smith, 

Brook House, Gilling East, 
YORK YO6 4JJ 


Price in UK: £10 per complete volume (4 issues); 
£12 (overseas); £16 (airmail) including postage. 
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APL Product Guide — Full 


compiled by Gill Smith 


VECTOR's exclusive APL Product Guide aims to provide readers with useful 
information about sources of APL hardware, software and services. We welcome 
any comments readers may have on its usefulness and any suggestions for 
improvements. 


Pressure on space sometimes prevents us from printing the complete guide, 
however updates will always be listed. We do depend on the alacrity of vendors 
to keep us informed about their products. Anyone who is not included in the 
Guide should contact me to get their free entry — see address below. If we have 
your Email address listed, we will mail you a reminder in good time for the next 
issue. 


We reserve the right to edit material supplied for reasons of space or to ensure a 
fair market coverage. The listings are not restricted to UK companies and 
international suppliers are welcome to take advantage of these pages. 


For convenience to readers, the product list has been divided into the folowing 
groups ('poa' indicates ‘price on application’): 


+ Complete APL Systems (Hardware & Software) 
- APL Interpreters 

+ APL-based Packages 

+ APL Consultancy 

+ Other Products 

+ Overseas Associations 

+ Vendor Addresses 


Every effort has been made to avoid errors in these listings but no responsibility 
can be taken by the working group for mistakes or omissions. 


We also welcome information on APL clubs and groups throughout the world. 


All contributions and updates to the APL Product Guide should be sent to Gill 
Smith, at Brook House, Gilling East, York, YO6 4JJ. Tel: 01439-788385, Email: 
100331.644@Compuserve.com 
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COMPLETE APL SYSTEMS 


COMPANY PRODUCT PRICES(£) 
Dyadic 18м RSj6000 MDG20 1179 
IBM RS/6000 мозго 13817 
IBM RS/6000 моз2о 22658 
IBM 5/6000 MD520 37114 
IBM RS/6000 мо5з0 72054 
IBM нзубооо MD540 122842 
Interprocess 
Systems, APL2 Devt Workstation pos 
Optima IBM Compatible poa 
APL INTERPRETERS 
COMPANY PRODUCT PRICES(£) 
APL Sottwere APL"Plus/PC Release 10 450 
Runtime poa 
APLPlus i 1295 
Runtime. 
Dyalog APL 1000-10 000 
АРШРС pos 
‘Atiantis Software Analyte Platform (К) poa 
Beautiful Systems  Dyalog APLW for Windows poa. 
Dyalog API for Unix pos 
The Bloomsbury Software Company 
APL'PLUS PC Rei 11 250 
APLPLUS Ill Windows эш 
APL'PLUS Ii for DOS 750 


DETAILS 


APL POWERstation (Greyscale) 27.5 MIPS, 7.4 Mflops RISC 
Processor ВМЬ RAM, 120M Disk 

197 1280x1024 Greyscale Graph Display AIX, OSF Mott, Dyalog 
APL (Huser) 

APL POWERstation (Colour) 27.5 MIPS. 7.4 Mfops RISC 
Processor Bhe AAM, 72040 Disk 

16° 1280x1024 Colour Graphics Display AIX, OSF МОМ, Dyalog 
APL (taser) 

‘Advanced APL POWERstation 27 5 MIPS, 7.4 Морз RISC 
Processor 16МЬ RAM, 320M Disk, 150Me Tape 

187 1280x1024 Colour Graphics Display AIX. CSF Mott, Dyalog 
APL (user) 

APL POWERsystem (G-users) 27 5 MIPS, 7.4 Мора RISC 
Processor 16М RAM, З20М0 Disk, 150Mb Tape CD-ROM 
Drive, 16 Ports 

AIX, Dyaiog APL (2-8 user licence) 

APL POWERsyatem (16-usors) 34.5 MIPS, 10.9 Mflops RISC 
Processor 32МЫ RAM, 1.3465 Disk, 2 300 Tape CD-ROM 
Drive. 18 Ports 

AIX, Dyalog APL (B+ user licence) 

APL POWERsystam (22-users) 41 MIPS, 13 Mops RISC 
Processor ВАМЬ RAM, 1.78 Disk, 2.3Gb Tape CD-ROM Drive, 
32 Ports 

AIX, Dyalog APL (8+ user loonce) 


Mainframe APL2 supported on a PS/2 via а co-processor card 
with 16Mb of memory running VM/ESA (370 mode). A complete. 
‘system Inciudes а Р5/2, a P/370 coprocessor card, and 
‘oftware licenses lor VWESA, APL2, GDDM and tho full ine of 
Inarprocess APL2 enhancements. 


‘Complete PC-based station. APL interpreters & sil support eq 


DETAILS 


БТС APL for iBM PCs & compatibies. Upgrades trom earlier 
releases also avajlabio. 


Closed version of APL'PIUNPC which prevents user exposure to 
APL 


AJ! the features of maintrame APL*Plus for your 385PC! 


2nd generation APL for Unix systems 
IBM's APL 2 for the PC. 

Kis an APL-iko language 

US Distributor of Dyalog APL products trom Dyadic. 

‘See Dyadic listing for product detaiis, 

'STSC’s tuii featured APL for IMs and compatibles - Version 11 
gives free runtime. 


The new 32-bit native Windows APLPLUS. Develop in 
Windows. and distribute APL applications with no runtime. 
charges. Reasonable migration charges from APL'PLUSIPC and 
APL"PLUS It. 


Now that APL*PLUS Ill for Windows is available, the facil lor 


creating Windows applications In PLUS І has been removed, 
and the prios reduced. 
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APL'PLUS Il for UNIX poe 
APL'PLUS VMS poe 
APL'PLUS Mairtrama. poa 
Dyadic Dyalog APL for 005/388 9% 


Dyalog APLW for Windows 995 


Dyalog APL for Unix 996-12,000 


JAC/Hurman Inerta 
LAPUMac 13 
FAPUPC or clones 8 
I-APLBBC Master 8 

в 


LAPYArchimedes 
Strand Software Inc 


APL Lid 


IBM APL Products TryAPL2 


АР\2 PC (US Version) 


Pg 


АРІ2 PC (European Version) 


АРІ2 for 05/2 Envy Editon — $185 
 APL2 for OSI2 Advanced Edition $650 


 APL2 for Sun Solaris ЕТ 
АРІ2 tor AIX 8000 
APL2 Version 2 
 APL2 Application Envt Vn2 
APL*PLUSIPC 
Dyelog APL 
АРС 
Iverson Software Inc. J Professional (nc runtime) 
J Persona! Editon 
J Personal (diska only) 


8 


Insight Systema 


APUWIN 
APL Reference Manuel 
System Kt 
J Source Собе 

MasterWork Sofware Manugistics Products end ISI 


3:008 SZETERERE 


'STSC's 2nd generation APL for ай major Spare end Risc Unix. 
workstations, 


2nd generation APL for DEG VAX computers under VMS. 


Enhances VS APL with many high performance, high 
productivity features. For VM/CMS and MVS/TSO offers simple. 
upgrade from VS APL. 


Second generation APL for DOS Runs in 32-bit mode, supports 
very large workspaces, Unique "window-based" APL 
Development Environment and Screen Manager. Requires. 
386/486 based PC or PS/2, at least 2МЫ РАМ, EGA or VGA, 
DOS 3.3 or later. 


As above, plus object-based GU! development tools. Requires 
Windows 3.0 ог Inter. 


‘Second generation APL for Unix systems. Available for Altos 
Apolo, Bull, Dec, HP, IBM 6150, IBM RS/6000, Masscomp, 
Pyramid, NCA, Sun and Unisys machines, and for PC» end. 
РС/23 running Xenix or AIX. Oracle interface availabla for IBM, 
‘Sun and Xenix versions. 


Macintosh version of LAPL 
150 conforming interpreter. Supplied only with manual (see 
"Other Products’ for accompanying books). 

As above 

Аз above 

‘Strand Software Inc has the sole selling rights to Iverson. 
‘Software Inc products. APL stocks a few of these (mainly 
АРИМ and the personal J products and books), bulis no 
longer an agent. 


 APL2 for educational or demonstration use. Write, fax or Emak ko 
APL Products, specify disk size desired, 


Product No. 5799-PGG. РРО Number R0411. 
Order trom t-800-IBM-CALL 


Product No. 5804-260. Part number 38F 1759, 
From ail IBM dealers, including MicroAPL. 


Part No B9G1556. 


Pan No 6961867. Contains all cities of the Entry Editon plus: 
DR? Interface: co-operative processing TCP/IP Interface; tools 
tor writing APs; TIME ону 


Product No, 5848-065. 
Product No. 5765012. 
Product No. 5888-228. Full APLZ system for 5/370 and S/390 
Product No. 5888.229. Runtime environment tor APL2 packages 
APL systems markoted and supported 

from: Dyadic. Manugistics, IBM 

ungar: Windows, OS2 end Unix 


The text of the manuals is available аз Windows Help, so the 
paper copies ara no longer а necessity. 


For S80/PC under Windows 3.1 
Documentation for ай the above. 

362 diskette with manual "Introduction and Dictionary" 
Ful С souros code pius 100-page book 

New Zealand distributor 
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MieroAPL. ‘APL 88000 Level i 2000 
АРІ 50000 Level 11 2500 
Арі 68000 1500-8000 
‘APL 8800 Lovet i 
Mac, ST. Amiga. ar 
Mac, Amiga 200 
‘APL 88000 Level (i 
sr 170 
Amiga 200 
Mec 520 
APL'PLUS Rel 10 450 
APL'PLUS II V 40 1395 
Optima APLYPLUSIPC 389 
APUPLUS I 950 
APL'PLUS II PC Developers К! poa 
Dyslog APL 990 
RE Time Tracker Оу APL'PLUSIPC pou 
APL'PLUS IVDOS 
APLIPLUS IWIN 
APL'PLUSIUNIX 
Solon Associates SHARP APL tor MVS Poa 
SHARP APL for Unix pon 
Strand Software Canada 
AX APL'PLUS Products poa 
Dyadic and ISI products pon 
USA 
Dyadic and IS! products poa 
Uniware APL'PLUSIPC 495 
Run-Time et 
APL'PLUSIUNX cat 
APL'PLUS It can 
APL PACKAGES 
COMPANY PRODUCT PRICES(£) 
AdeptabeSysems FLAIR poa 
APL-385  APL-385 for APL"PLUSIPC 50 
FSM-385 
DRAW-385 
08-385 
GENES, 


Firet generation APL with numerous enhancements. Multi-user 
version (Unix, Miraga, MCS}. 


‘Second generation APL. Nested arrays, user defined operators, 
Selective specication etc. Muhi-user version (Unix, Mirage. 
NCS) 


Secona-generation APL. Nested arrays. user defined operators, 
‘Selective specification, atc. Ми. изат AIX version with lui 
OSF/Mot support 


First generation APL. Single user. Iul windowing Interface, 
software floating point support 


First generation APL. Single user. ful windowing Interface. 
hardware floating point 


Second generation APL. Ful windowing interface, software 
floating point support 


Second generation APL. Full windowing interface. Hardware and 
‘software floating point support 


Second generation APL. Full windowing interlace Hardware and 
software Поабе point suppor. 


Complete APL'PLUS and Statgraphics product range and user 
support for Finland 


Tor IBM MVS maintrames. 
or IBM 95/6000 and Sun SPARC 


AN APLPLUS products Including upgrades and educational. 


STSC's full faature APL for IBM PCIXTIAT. Compaq, Otvetti 


Closed varsion of APL*PLUSIPC which prevents user exposure 
10 APL. 

STSC's full fonture APL for UNIX based computers 

STSC's ful feature APL tor 386 machines. 


DETAILS 


Finito loader and interactive rescheduler. Customisable fui- 
function scheduling system. (Avallable outside Australis by 
special arrangement only.) 


including 
Screen development 
‘Screen design. 
Feintional W.S. 
Miscellangous Utiles 
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The APL Group 


APL Software Lis 
(maintrame) 


{microcomputer} 


Вели Systems 


Quaiedi 


ROS 
Lt 
REGGPAK 


POWERTOOLS 


REGGPAK 
ROS 
ASF FILE 


NAT, FILE 


Ове РШЕ 


The Bloomsbury Software Company 


(for VSAPL) 


(or APL2) 


Causeway 


Cinerea AB 


CODEWORK 


CYBEX AB 


HMW. 


HRH Systems 


Enhancements & Sharefle 
Compiler 

SherefiijAP 

‘Causeway lor Dyalog 
Causeway for APLPLUS Ш 


Rain Graphics Workspace. 


ORCHART 


HELM 


APL GraljPC 

APL Grat РС 

Liity Functions APL2 
ишу Functions РС 
AXTRA 


 APL'PLUS Utilities 


= 


ТТУ 


H 


з sc ig 


H 


E 


323228888 


Electronic Data interchange (EDI) translation software for the 
PC, with strict compliance checking. 


Raiation Dala Base System 
Project Management System. 
Regression Analysis Package. 


 Assombler written replacement lunction for commonly used 
GPU-consuming APL functions, Indudes a Forms Processor. 
Regression Analysis Package 

Relational Database System 


Dyeiog APLIW auxiliary processor for access to APL*PLUS/PC 
APL component fles (ASF), 


Dyalog APL/W auxilary processor which emulates the 
APL'PLUSJPC quac-N native fle subsystem for access to the 
DOS the system. 

Оувю APLIW auxiliary processor fo efficient block mode 
coss to dBASE format fles. Designed to get large amounts of 
dete in and ovt of dBASE. Not suited for random access ю smal 
‘emounts cf deta (t does not handia keys). 


Component fles, quad-tunctions & nested arrays for VSAPL 
under VM/CMS & MVS/TSO 


The Fist APL complet 


STSC's shared access component fle system tor APL2. 
‘Comparable to all APL*PLUS fie systems. multi-user storage ot 
APLZ arrays wit efficient disk usage. 


Manuals and Ciass-managemeni utes for the Causeway 
platform suppled free with Dyalog АРУУ 


Software, ciass-management utiles and printed documentation 
Tor Causeway under APL'PLUS Il! 


Full on-tne documentation (Windows Help) end handy-reterence: 
card for the Fain business and statistical graphics workspace 
Supplied with Dysiog/W. 


‘Organization chart package for IBM APL2/PC. Ful & heavy 
commented source code cluded - ree iniegraion nto other 
"applications. NB: ASCII output wit tne-crawing (sec graphic) 
charactars tor boxes 


Decision Support system for top management. Developed in Italy 
over 7 years. Requires APL mainframe or APL"PLUSIII. Optional 
modules: EIS, Excel interface, ОТР output via LATEX, output on 

map background, 


Presentation graphics for APL*PLUSYPC (CGI) 
Prosentaton graphics for APL*PLLUS IPC (CGI) 

For APL mainframe; Inci. a very fast search. 

‘Same package for APL'PLUS IVPC. 

Front-end Foreign Exchange dealing / pos keeping 

Arbitrage modeling 

Baskoi currency modeling 

pull-down menu tor APL"PLUSPC 

‘Software to vanster workspacas between APL'PLUS and Sharp, 
‘and between APL‘PLUS and I-APL. Software lo impor IBM ATF 
Мез to APLPLUS. 


Public domain software, unlock locked tns, a user-riendly 
aliematvo to locking, fne of mathematical physics, menus, and 
others 
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\ACiHuman Interteceel АСга 5 

Смех 15 

LAPL Lid Educational workspaces 5 

IBM APL Products A Graphical Statistical System $250 

(sss) E 

$2500 

Impetus Lid imo ке 
INFOSTROY  APL'PLUS/Xbase interface 

01306 Version 2) $188 

(PC Version 2) эв 

(DLL Version 1) $198 

Insight Systems  IUTILSXP 2095 

лз LI 

Wincom os 

20220.X2X pos 

САРІ Client soa 

SOAPL Server pos 


Interprocess Systems APLZ Development Workstation ров 


от 53000-5000. 
(valntrame) AFM $15300 
Enhanced Format $2575 
PowerCode ‘$2000 
мова $1500 
JAD Software JADSMS 150-500 
Linge Allegro. FRESCO Business Graphies ров 
API28/PC. ров 
APIPTIPC poe 
APIIgPC pos 
FACS poa 
TOPR pos 
Mercia Lost s2 poa 


E 
8 
Н 


Graph plottng for LAPL Mac 
Spoken APL characters for APUMac: 

PC format disks with the examples trom: Thomson. Espinasse 
(ts 1-4). Kromberg, Jizba & FinnAPL All the exemples to save 
your fingers 

for DOS, Product Number 5784-009 

tor Workstations (05/2, Aix, Solaris), Product Number 8784-082 
for CMS, Product Number 5764011 


Corporate Modeling and Reporting System. 


‘Complete package written in C. Comparable with the data, Index 
& memo fles of FoxPro, ABASE, & Clipper. Multiuser support. 
No DBMS license required. 


As above for APLPLUSIPC. 


The same in a DLL form! Gives your Windows applications. 
advantages of DLLs. 

Cross-platform tity Nbrery including simple OS calis (DIR, 
COPY, DEL, RENAME) end DATE functions. For APL'PLUS Il, 
 APL2 and Dyalog APL under Windows, OS/2 and Unix. 


APL Spreadsheet rtarfaco. "Device-independeni" spreadsheet 
driver supporting Excel, 123 and Quatro-Pro for Oyalog APLW 
‘Asynchronous comms package for Dyalog APLW 

Advanced APL syntax analysis and conversion packages trom 
‘Sharp and APL2 to Dyalog, and between any two APL 


Intertace from APLPLUS І, APL2 and Oyalog (Windows, О5/2 
or Unix) to most SQL databases over most networks, 


Makes APL'PLUS Il, APLZ or Oyalog APL (Unix) avaliable as. 
‘SequeLink servers. Can be called from SQAPL cilents or other. 
‘applications such as Excel, C++, Smaltak, Visual Basic. 


Full screen APL2 editor with Immediate APL execution, nnd tull- 
screen debugger 


High performance component and keyed Mie system (VS APL 
and APL2) 


А QuacFMT data formatter for VS APL and APL2 
External functions for APL2 
Full-screen Workspace Organizer tor APL2. 


‘Software management system for APL"PLUS Il based on 
hierarchical databases: includes full-screen interface and stand- 
‘alone functions. Price depends on number of users. 


Fast and Easy Business Graphics DLL 
ФООМ iniertace tor Dyalog APLIW 

ODBC intertace for Dyalog APLIW 

TCP/IP imterfoco (or Dyalog APL/W 

ЕММА-4кэ зносе to DR2 or ODBC databases 

APL Code end Appilcation Managemen for Dyalog АРУУ 


Logistics management system for 386/486 & RISC computers. 
Sales Forecasting, Inventory Management, Master Scheduling, 
Distribution Requirements Planning, Sales & Operations 
Planning. 

А modular library of tools to teach and explore state-of-the-art 
materiais managament concepts. Developed by R.G. Brown. 
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RE Time Tracker Oy ИТУ poa 
ов poe 
Solon Associates 10005 poa 
MAILBOX pos 
VIEWPOINT: poa 
UNIWARE 

бог mainirame) STSC's ENHANCEMENTS poa 
STSC's SHAREFILE pou 
TOOLS & UTIUTIES poa 
EXECUCALC poa 
(forAPL'PLUS/PC) APL Debugger 2.1 FF1950 
FF9750 
Menus 30 FF2450 
FF1250 
ETATGEN 20 FF1950 
FF9750 
UNITAB 20 L2] 
кешю 
UNIASM 3.0(sitelicense) FF4650 
UNISTAT 51 FF2900 
MorAPL'PLUS II)  UNIWARE Toaikt 41 FF29000 
APL Debugger I1 2-1 кгз 
ШЗ 
Menus i 4.0 FF3950 
FF1970 
ETATGEN ll 20 FF2050 
FF14750 
UNITABII 20 FF6950 
FF34750 
UNISTAT Pus 5.2 FFa300 
Warwick University BATS 250 
Faa tree 
Zark APL Tutor (PC) оз 
APL Tutor (MF) $5000. 
Zark ACE Ld 
APL Advanced Techniques... $58.95 
Communications $200 pe, $500 mt 

APL CONSULTANCY 
COMPANY PRODUCT PRICES(£) 
Active Workspace APL Programming poe 


TMT-Team Оу User Intertace Toolkh for APLYPLUS Il and 
PLUS Ill under Windows. Comprehensive spreadsheets, 
repiiceted leis, special field typos, etc, 


TMT-Team Oy's database interface for APL*PLUS N & PLUS IN 
under Windows. Interfaces Yo almost twenty diferent databases. 


Application Development Environment 
Electronic Malt 
Report generator with interfaces io DB? and MVS data 


‘Quad-tunetions & nested arrays for IBM VSAPL 
‘component fles for IBM VSAPL and for IAM APL2 


FILEPRINT, FILESORT, FILECONVERT 
FILEMANAGER(EMMA) STSC's database package 


Maintrame spreadsheet compatible with VISICALC and part of 
LOTUS 1-2-3 under VSAPL{VM or TSO) 


‘A visual APL debugger to help develop applications 
(sit license) 


Complete sel of hierarchical menu utiles 
(she license) 


Ал APL*PLUS sproadsheet-ike data entry and validation system 
(she license] 


Assembler uites 10 speed up APL'PLUSIPC applications 
Date analysis add-on module for Statgraphics 


(she license only). Relational database system end complete set 
of lies for APLPLUS II development 

A visual APL debugger to help develop applications 

(she icense) 


Complete set of hierarchical mouse-criven menu utilisa 
(site license) 


Page layout report generator 
(aho license) 


An APL'PLUS spreadsheet.ke data entry and velldation system 
{айе canse) 


(Data analysis add-on module for Strigraphics. 


Menu driven system for time series analysis und forecasting 


using Bayesian Dynamic modelling. Prios is reduced to £35 for 
academic institutions. 


Training program for Ihe above. 


APL computer-based raining. Avaliable for APL'PLUS PC & 
 APL"PLUS II. Demo disk $10. 


 Meintrare version. 


APL continuing education. APL tutor news and hotine phone. 
‘support. 


&88рр. book, (ISBN 0-9619087-07) including 2-disk set of ину 
tunctons (APL'PLUS PC format) 


Move workspaces or es between APL environments 


DETAILS 
‘Short or long-term consultant available 
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жне Consultancy E 
Andrews Consuiiancy poe 
APL People Consutancy poa 
Bloomsbury Soltwere Consuttancy 300-750« VAT 
Camacho Consutaney poa 
Ray Cannon Consultancy poe 
Causeway Consultancy and Training poe 
Paul Chapman  Consutancy pos 
Daud Crossley Consultancy юе 
Poter Cyriax Consultancy 100-150 
120-200 
160-300 
Dogon Research Consultancy poe 
Dyadic Consultancy pon 
eas Consultancy poa 
Evestic AB. Consultancy pos 
General Software Consultancy from 120 
Greymande Assoc — Consulting poa 
нмм ‘Consuttancy poa 
Micheal Hughes Consuancy poa 
1AC/Human Intertaces Consultancy 350 
Documentation 100-200 
Training poa 
INFOSTROY Consultancy poa 
meghtSyseme  Conauhancy poa 
intaigent Programs Consultancy poa 
Documentation. E 
Training poa 


‘Development, maintenance, conversion, migration, 
‘documentation, of APL products in ай APL environments 


APL programming and analysis, specialises in tree-procassing 
gore. 

Consultants avaiable at al ievols Expertise in APL system 
design. project management, prototyping, Hnanclal applications, 
‘decision support systems, MIS, Inks to non-APL systems, 
documentstion, өк. 


Manuais: fessibiity reports and estimates; analysis and 
programming, APL and MS Windows applicetions; Sharp, 151 
APL, APL*PLUS, APL2/PC and other АР\ з spoken. Fixed price 
‘systems a specialty 

APL, C, Assembler, Windows, Graphics: PC and mainframe, 
Management and upkeep of the Causeway development 
environment for Individuals and large APL sites. On-site traning 
tor Causeway/Dyaiog and Causeway/Plus Il 

241v programmer: APL. C, Assembler, Graphics; PC, minl, 
тамгат and network 


Broad experience in many APL environments 


Junior Consuttant 
Consultant 

‘Senior Consultant 

APL Systems consultancy, design, implementation, suppor, 
documentation end maintenance, Al dlalects with special 
‘emphasis on APL2 and Dyalog АРАЛАУ. 

APL and Unix system design, consultancy, programming and 
vanning. 

Systm prototyping: ak ypes of information system, engineering. 
softwere, graphics and decision support systems APL'PLUSIPC, 
АРІ2, Dyalog APL 

Excellent treck record from 10+ years of APL applications in 
‘banking, insurance, and education services. Al dialects, 
Platforms and project phases. SQL expertise 


Company reporting, business graphics, Windows applications 
мап Dyalog APLIW. 


Synem design consultancy, programming. HMW specialize kn 
белеп and prototyping work. 


Consultant with 10+ years experience with various APL 
Intarprutars and C. 


APL on Macintosh & PC. НСІ design. VDU ergonomics: 
ECiHeahth & Sefety compliance. 


On-Ine assistance, product demos & mock-ups, manual wring: 
orsign languege software locakzaton, 


Using -APL tor courseware & distance leaming materials; Mac 
programming In С, APL & HyperCard, 


APL"PLUS & Windows consutancy. Porting of software written 
in C into APL'PLUS. 


Experts in APL conversions between any combination of: 
 APL*PLUS. APLZ. Dyalog APL and Sharp APL We ara also 
experienced right-sizers, comfortable with networks and 
relational databases (that also means when NOT to use SQL) 
and cient/server development in APL, C nd Visual Basic. 


‘Systems development. enhancements, support 
Preparation of new manuais, rewriting of existing materials, 
Training for APL experts through to non-technical system users. 
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JAD Sotware Consuttancy poa 
кази ‘Consultancy pos 
Ungo Allegro USA Consultancy poa 
MasterWork Software Consunancy юа 
MieroAPL Consultancy poa 
Enis Morgan Consultancy 250-500 
Optima Constancy poa 
QB On-Line Consutancy эю 
RE Time Trucker Oy Consultancy pon 
Rex Swain Consultancy pon 
Rochester Group  Consutancy 
Snake Isand Research Inc 

Consultancy. pon 
Stand Sofware ‘Consultancy poa 
Sykes Systems Inc Consultancy pos 
Шотага Consuttancy (Senior) FFjdeySO00 

Consultancy (Senior) — FF/day 7500 

Training FF10000 
Wickie Computer Consultancy poa 
OTHER PRODUCTS 
COMPANY PRODUCT PRICES(£) 
Ados Employmont poa 
APL Prope Employment agency pon 
Bloomsbury Software Training pon 
сотор Comte Logger s25 95«p&p 
HMW Employment. pos 
HRH Syetems APL lessons 

The BESAPL: ape 
LAPL Las An APL Tutorial 


‘An Encyclopaedia of APL (2d Ed) 
APL in Social Studies. 
LAPL instruction Manual (2d Ed) 


Vol.12 No.2 


Systems design and development, project management, 
technical manuais, financial and actuarial expertise in APL. 


AN APLs, all environments. Design, anaiysis, coding. 
maintenance, documentation, training, interfacing. 


General APL consulting: 
Migration and Downal2ing; Performance Tuning. 


Consuhancy 

Technical & applications consultancy. 

‘Business Forecasting & APL Systems. 

А range of consultants with 3-15 yrs APL PC and mt experience, 
Specialising in Banking, Financial & Planning Systems. 
‘Specialised in compranensive APL Windows user interfaces, 
APL Multimedia, APL 1 AP! level interfacing for Windows, 
Windows applications, DLLs & databases. 


Independent consultant, 20 years experience. Custom software. 
development & training, PC and/or mainframe. 


Specitilse in MIS using Sharp APL 


APL interpreter and compiler enhancements, Intinsic functions, 
performance consulting. APL compiler for sec and parallel 
systems now under est 

Advice on migrating to end trom all favours of APL ала 
hardware piatiocme. Fullscreen interface Implementation, APL 
ties, benchmarking. efficiency analysis, actuarial software, 
system development bois, valuation, pricing end modeling 
systems. 

Completa APL services speciakting in sudit optimisation and 
Conversion of APL systems. Excellent design skits. All dialect 
and platforms. 17-23 years experience. 

Consultancy from people with at east 8 years APL experience, 
Advion and brining In Windows programming with APL*PLUS II 
5-day class on Windows programming with PLUS И version 4 0 


System design, consukancy, programming and documentation, 
Especially project management and decision support systems 


DETAILS 

Contractors and permanent employees 
‘Employees placed at all levels. 

Contact the company for detalis. 


APL*PLUS Il comic-book inventory system. Sharewara version 
available on America OnLine. 


Contractors and permanent employees placed. 


On-screen interactive APL lessons for APLYPLUS, ТАРС, 
‘Sharp and APL — in English or French. 


703-528-7617, 1200-14400b, N-8-1, 24 hours. APL educational 
material is downloadable frae. An additonal 30 megs of APL 
software for APL'PLUS, PLUS Il, IBM, Sharp & APL is 
‘avaliable to subscribers (cost is $24/yr). Selection avallable on 
disk for $15 post-paid. Free on-disk catalogue. 

45pp by Alvord & Thomson. 

228pp by Helzer 

38pp by Traberman 

55рр by Camacho & Ziemann 
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APL Programs for he Mathematics 
Classroom (Springer-Verlag) 10 185ррһу Thomson 

Programming in J 10 75pp by Ken Iverson 

Апттезс 12 Ttüpp by Kan Iverson 

Tangible math $ 36pp by Ken Iverson 

Shap APL ایا وی‎ 

APL Press Books А comprehensive selection o! early APL tterature 


Piasso noie areas packing charge of £9 per ower 


Көзгө Employment poa Permanent and contract, home and abroad. From individual 
Placarnent ko supply of completa project wams. 
Software Library poa Low-cost software distribution service. сай ог detalis. 
Renaissance 
Data Systems Booksellers The widest range о! APL books avaliable anywhere. Seo Vector 
advertisements, 
Solon Associates — MVSUNK poa Interface trom Sharp APL (Unix & MVS) to non-APL data and 
software In the MVS environment 
эзш. poa High-performance D82 Interface tor Sharp APL (Unix and МУЗ). 
OVERSEAS ASSOCIATIONS. 
GROUP LOCATION JOURNAL OTHER SERVICES Ann.Sub. 
ACM/SIGAPL Intemational Quote Quad. 
APL Вау Area USAN. Calfomia APLAUG Monthly Moetings (2nd Monday) 915 
APL Club Austria Austria З Quarterly Meetings. Z200AS(indiv), 1000AS corp) 
APL Club Germany Germany APL Journal ‘Semi-annual meetings. омео 
APL interest Group South Africa. : š 


Ass. Francophone pour 


la promoton АРЫ. France Las Nouvelles ГАРЬ. 
BACUS Beigum APLCAM Conferences & Seminars £18 ($30) 
Chicago SIG Chicago 

Hartord Group Harford. Connecteut, USA 

Capta РСС  Washingion, DC. Monhor Monthly meetings, occasional classes hoe 
Danish SIG Denmark 

июл APL Assoc. Holland Min congress, APL ShareWare Initiative 

Finnart. Helsinki. Finland — FInnAPL Newsietter Seminars on APL OOF IM (private). 3O(student). 1000 (Co) 
Japan APL Assoc. Tokyo. Japan 

Melbourne APLUG Melbourne. Australia Quarterly meetings free 
New Yok SIG New York USA 

Potomac SiG Washington DC, USA Free monthly meetings 

Rochester APL Rochester New York 

RomeltaySIG Вота, my 

SE APL Users Grp Айал, Georgia SEAPL Newsletr 

SOCAL Soutem Cantoria Seminars $15 ($5 students) 
SovAPL Obninsk. Russia 

SwoaAPL Sweden SwedAPLNyt — Semkennual meetings, seminars ЗЕК 75 
SWAP. Texas. USA SWAPL se 
Swiss APL (SAUG) Bom Part of Qty 5100 SFBO (SI) + SF20 (SAUG) 
Sydney APLUG Sydney. Australa Epsilon Monthly Meetings 

Toronto SIG Toronto, Canada Gimme Arrays! Monthly Meetings, APL skills database, J SIG, Toronto Toolkit — $25 
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VENDOR ADDRESSES 
COMPANY CONTACT 
ACMISIGAPL Donne Baglio 
‘Active Workspace Lid Ross D Ranson 
Adaptable Systems. Lols & Richard Hit 
Adaytum Systems 

лове Bernard Smoor 
Andrews Dr Anne D Wilson 
APL-385 Adrian Smith 

APL Bay Area Users Group 

APLBUG Lewis H. Robinson (Sec) 
APL Club Austria Erich Galt 

APL Club Germany Dieter Lattermann 


The APL Group Inc ‘Stuart Sawabini 


APL Interest Group, South Africa 


Mike Montgomery 
APL People / Softwere ЈИ Moss 
‘Association Francophone pour 
la promotion d APL Dr. Gérard Langlet 
Atlantis Software Arthur Whitney 
BACUS Joseph de Kerf 
Beautful Systems, Inc. Jim Gott 
The Bloomsbury Software Co Lid 

Peter Day 
Camacho Anthony Camacho 
Ray Cannon 
Paul Chapman 
Causeway Graphical Systems Ltd 

Adrian Smith, 

Duncan Pearson 
Chicago SIG Lary Mysz 
Cinerea AB. Во! Komemark. 
CODEWORK Mauro Guazzo 
ComLog Software Jeff Podnesu 
cPcuo Lynne Зал 
David Crossley 
CYBEX AB Lars Wentzel 
Peter Cyriax Systems Petar Cyriax 
Danish User Group Per Gjeriot 


ADDRESS & TELEPHONE No. 

АСМ, 1515 Brosway, New York, NY 10038 USA Tel +1 (212) 626-0606 
Emali: baglio @ acm org 

Moulsham МИ Centre. Parkwsy, Chelmsford, Essex, CM2 7PX, UK. 

Tei: 01245-49647; Fax: 01245-490046, 


49 First Street, Black Rock 3193, Australis. 
Tei: «613 589 5578 Fax: +61 3 580 3220 


13 Great George Street BRISTOL BS1 SRA UK Tol: 0117-921 5555 


50, 4128 BZ Lexmond, Netherlands. 
Тө +31-3474-2337, Fax. *31-3474-2942 


23, The Green, Acomb, YORK YO2 SLL, UK Tet: 01904-792670 


Brook House, Giling East, York YOO 4 UK. Tel 01439-788385 
Fax: 01439.789194 Emai: 100331.544@ compuserve. com 


1100 Gough St, Ap! 14А, San Francisco, CA $4109, USA 
Tet +1 (415) 928-2058 Email: trgp?1a6 prodigy.com 


IBM Osterreich, Obere Donaustrasse 95, A-1020 Wien, Austria 


Fheinstrasse 23, 0-69190 Walidod, Germany. 
Tet +49 6227-53469 Compuserve: 100332, 1451 


644 Danbury Road, WILTON, CT 06897 USA. 
Tet +1 (203) 782-3033 Fax: +1 (203) 762-2108 


Price Bog X11. Rivonia 2128, South Africa: 
Tet +27 (11) 803.7200 Fax: «27 (11) 603-9194 
Emai: тетот! spl oo za 


The Old Malthouse, Clarence SI, BATH. BA1 SNS UK. Tel: 01225-462602 


SOM, C.E Saclay, F-S119-GIf sur Yvette, France. Fax:+33 1 69-08-79-63 
1105 Harker Avenue, Palo Alto, CA 94301 USA 
Racinberg 72, В-2570 Duffel, Beiglum. Tel: «32 15 31 47 24 


308 Оа York Road, Sune 5, Jenkintown. PA 12046. USA 
Tol +1 (218) 886-2636; Fax: +1 (215) 855-4868 


3-6 Aired Place. Bloomsbury, London WCTE 7EB UK. 
Tet 0171-438 9481; Fax: 0171-436 0524; CompuServe: 10010,1487. 


11 Auburn Road, Redland, Bristol BSS SLS UK. Tel: 0117-9730038. 
emai: ecernecho @ cix compulink.couk. Reuternet (Sharp): ACAM 


21 Woodbridge Rd, Blackwater, Camberiey, Surrey GU17 OBS UK 
Tet 01252-874897 Emal: 100430.740 8 compuserve.com 


518 Lambs Conduit Street, London WC1N 3NB UK. 
Tet 0171-404 5401. Compuserve: 100843,3210 


5 The Maltings, Castlegate, MALTON, North Yorks YO17 ООР UK 
Tol: 01853.896760 Fax: 01653.697719 Compuserve: 100285,1564 


83 Highland Drive, Chicago Heights, IL 60411 USA C'servo:73040,3032 
‘Skyttegaten 25, S-193 00 Sigwna, Sweden. 


Corso Cairo 32, 10123 Torino, aly. 
Tol: #39 11 885168 Fax. «39 11 812 2052 


РО Box 5570, Derwood, MD 20855 USA 
Tel: +1 (301) 990-7083 Emall: ot @sottmed com 


Capital PC User Group, 51 Monroe Street, Suite PE-2, Rockville, 
Maryland 20850-2421, USA. Toi: +1 (301) 782-9372 Fax: (801) 762-9975, 


187 Le Tour du Pont, Quarter Le Mourra, 84210 ST DIDIER, France 
Tet +33 90-68-08-87 


Gruvgetan 958, S-421 30 V. Frölunda, Sweden. 
Tel +48 31-46 37 40, Fax: «46 31-45 24 23 


22 Horeford Road, London W2 4AA UK Tet: 0171.229 5344 
Ета! gierper@ inet ипк 
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Datatrade Lis тал Tomiin 
Dogon Research Dick Bowman 


Dutch APL Association Bernard Smoor (Soc) 


Dyadic Systems Lis Poter Donnelly 
Е A S Associates Frank Evans 
Evestic AB Ole Evoro 
FinnAPL 

бета Sofware Lid ME. Martin 


Greymante Associates ^ George MacLeod 
Hartford CT Group Bob Pomeroy 


H.M W Trading Systems Lid Stan Wilkinson 


HRH Systems Dick Holt 
Michael Hughes 
lAC/Humanintefaces lan A Clark 
LAPL Lad Anthony Camacho 
(lor queries, order forms) 
JC Business Services 
(for pre-paid orders only) 
IBM APL Products Nancy Wheeler 
impetus Lid Себе Heddle 
INFOSTROY Alexei Miroshnikov 
Insight Systems ApS Morten Kromberg 


Inveligent Programs Lig Mike Bucknall 
Interprocess Systems inc. Stella Chembertein 


Iverson Software inc. Eric Iverson 
JAD Software Оема Crossiey 
Japan APL Association 

Kestrel Consulting Mark Haris 


Lingo Allegro USA Ine. Vere H. FF 
MasterWork Software Lid Fraser Jackson 
Melboume APL Group Harvey Davies 

Mercia Sowere Lid. Gareth Brentnall 


MieroAPL Lid. David Eastwood 


1 & 2 Storing Business Park, Salthousa Road, Brackrnils, Northampton. 
NNA OEX UK. Tel: 01804-76041 


2 Dean Gardens, London £17 JQP UK Tet 0181-520 6234 
Email: bowman api demon co.uk 


Postbus 1941, 3430BH Nituwegein, Netherlands. Tel: +31 3474-2337 


Riverside View, Basing Road, Oki Basing, Basingstoke, 
Hants ОДА OAL UK Tel 01256811125 Fax 01256 211120 


19 Homesdale Road, Orpington, Kent BRS 1S UK Tet 01680824741 
Bertellusvagen 12А, S-146 38 Tullnge, Sweden зн»: «48 778 4410 


‘Suomen APL-Yhdistys RY FinnAPL RF. PL 1005, 00101 Meisinia 10, 
Finland 


22 Russell Rond. Nochnhoh, Middx, UBS 405 UK. Tet 0181-884 9537 
Bartum House, Ravens Lane, Berkhamsted, Hens, HP24 2DY UK 

Tel 01442.873065 Emañ: 100412,1305 @compuserws com. 

Mass Mutual Life, 1295 State St, Mallorop F455, Springliela. 

MA 01111 USA Tol: «t (412) 759.841 12838. 

Haminon House, 1 Temple Avenue, Victoria Erbankment. 

London ЕСАҮ ОНА UK. Tei: 0171-353 8900: Fax 0171-353 3325: 

Emal: 100020 26826 compuservs com 

3802 N Richmond St, Suite 271, Aringion, VA 22207 USA 

Tel +1 (709) 528-7824; Emaii: dick Not @ aem org 

28 Rushton Road, Wiibarston, Market Harborough, Laics. 1Е16 8OL UK 
Ter 01536-770098 

9 Н End, Frosierey, Bishop Auckland, Co. Durham DL1325X UK 

Tel 01388-527190. Compuserve: 100021 3073 

11 Aubum Road, Rediand, Bristol BSB BLS UK Tel 0117-9700038 
email acamecho@ сіх compullnk.co ик Heutemet (Sharp): АСАМ 

56 The Crescent. Milton, Weston-euper-Mare, Avon, 8522 BDU UK 


APL Producta, IBM Santa Teresa. эрі M4012, 555 Bally Avenue. 
бал Jose CA 95141, USA. Tel: +1 (408) A&S-APL (#2752) 

Fax +1 (406) 489-4488 Ema APL2@ vnetibm com Cserve: GO IBMAPL2 
Rusper. Sandy Lano, Ivy Hatch, SEVENOAKS, Kent TN15 OPO UK 

Tel: 01732805428 

ЭЗ Tulenin Lane, St. Petersburg 191196 Russia. 

Tol-+7 812-312-2673 Fax:+7 812-311-2184 Emal: aim@intostroy.spb. eu 
Nordre Strandvej 119A, DK-3150 Hellebæk. Denmark 
Tol-+45 42 10 70 22 Fax: «à5 42 10 75 74 Ema: nsightO net uni. ж 

9 Gun What, 130 Wapping High St. London E1 SNH Tet 0171-285 1420 
11860 Alpharetta Highway. Suite 455, Roswell, Georgia 30078, USA 

Tel: +1 (404) 410-1700, Fax: +1 (404) 410-1773 Cserve: 70373.2576 

33 Major Street, Toronto, Ontario, Canada MSS 2К9 

Tol: +1 (416) 825-0096. Fax: +1 (416) 408-7559 

580 Eyer Orive, #31 Pickering, Ontario, Canada L1W 387 

Tel: «1 (905) 897-1895 Fax:+1 (005) 831-5172 

29-2902 Hiromichi, Adachi, Tokyo 120, Japan 

Business & Technology Cente, Bessemer Drive, Stevenage, 

Herta 501 2X UK Tol: 01438-310155 Fax 01438-310131 

113 McHenry Road, Suite 181, Buffalo Grove, IL 60089 USA 

тее! (708) 450-7829 Рах +7 (708) 459-8501 Emai: nfo @ilngo com 
PO Box 56-096, Tawa, Wellington. New Zealand 

Tel пла Fax: +64 (4) 232-4440 Emak: 100242.2535 @ compuserve.com 
CSIRO Div Aim Ras, Private Bag No.1, Mordialloc, Vicioria 3195, Auszala 
Tel +613 566 7574 Fax «613 586 7800 Emak NGO dar сзко au 

Holt Court North, Heneage Steet West, Aston Science Park. Birmingham 
B7 АХ UK, Tel: 0121-950 5006, Fax: 0121-259 0375 

‘South Bank Technopark, 90 London Road. LONDON SE1 SLN UK 

Tel 0171-922 8808 Fax: 0171-028 1008 

Email: MercAPL O microapl demon. co uk 
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Elis Morgan 


New York SiG APL 
Optima Systems Lis 


Potomac APL SIG 


Ов On-Line Systems. 


Ells Morgan 


Nestor Nelson 
Paul Grosvenor 


John Martin 


Philip Bulmer 


Renaissance Dela Systems Ed Shaw 


RE Time Tracker Oy 
Rochester APL 


Tho Rochester Group Inc. 
Romesitely SIG 


ЗЕ APL Users Group. 
Shendel Systems Lis. 
‘Snake Isiand Research Inc 
SOCAL (South Caltoria) 


Soliton Associates 


SovAPL 
‘Strand Software Inc 


Flex Swain 


МАРЫ 
SwedAPL 
‘Swiss APL User Group 


Sydney APLUG 
Sykes Systems inc 


Toronto SIG 
Uniware 


Wickiffe Computer Lid 
Warwick University 


Zark Incorporated 


Richard Eller 
Gary Dennis 


Robert Pullman 
Marko Sacco 


John Manges 
Maurice Shanahan 
Bob Bernecky 
Bay Sykes Jr 


Laurie Howard 


Alexander Skomorokhov. 
Anne Faust 


Rex Swain 


Эше Увца 


(Christer инет 


Rob Hodgkinson 
Roy Sykes Jr 


Bon Best 
Eric Lescnsso 


Nick Tetor 
Prof. Jett Harrison 


Gary А. Bergquist 


Myrte Farm, Winchester Road. Stroud, Petersfieid, Hants UK. 
тө 01730-263843 


PO Box 188, NY 10185-0002, USA. 


Alport House, Purley Way, Croydon, Surrey CRO OXY UK. 
Tol: 0181-781 1812 Fax: 0181-781 1999 


61 Monroe Street, Plaza East Two, Rockville MA 20850-2421 USA 
Tel: +1 (901) 762-9372 Fax: +1 (301) 762-9375 Emailjam@acm org 


5 Surrey House, Portsmouth Rd.. Camberley, Surrey, GUIS 108 UK. 
Tel: 01276-955880 Fax: 01278-855301 


PO Вох 421, Georgetown, CT 08982, USA Tel: +1 (212) 864-3078 
РО Box 363, FIN-00101 Helsinió, Finland, Tel «358-0-400 2777 


оноп Associntas, 1100 University Avenue, Rochester, NY 14607 USA 
Emal: дві ipsalab tor soliton com 


50 $,Union St, Hochester NY 14807-1828, USA. 
Tok +1 (718) 454-4360. Fax: +1 (716) 454-5430 


Casella Postale 14343, 00100-Roma Trelo, пау 
Еле! marsac @ упо: бт сот 


991 Creokdalo Drive, Clarkston. GA30021 USA 
Chiltom House, High Street, Chafont St Gies, Bucks HPB 40H UK. 

18 РАА Straat, Ward's Island, Toronto, Ontario М5) 289 Canada 

Tel: +1 (618) 203.0054 Fax: +1 (415) 203-8099 

Emall: bernecky @ eecg toronto.edu 

Sykes Systems Inc, 4649 Wilens Ave, Woodland Hils. 

СА 91384-3812 USA Tol: +1 (818) 222.2759 Fax: «1 (818) 222.9250 
Soiton Associates Lid, Groot Blankenberg 53, 

1082 AC Amsterdam, Netherlands 

Tol: +31 20 648 4475 Fax. +31 20 544 1206 Ema ip @ soliton com 

РО Box 5061. Obrinsk-5, Kaluga Region, Russia 

Emalieskom api obrinsk su 

19235 Covington Cour, Shorewood MN 55331 USA 

Tel: +1 612] 470-7345 Emai: amlaust@ aol com 

8 South Street, Washington. CT 05798 USA 

Tel +1 (880) 888.0131 Fax +1 (860) 858-9970. 

Ела rswain @ ix netcom com CompuServe: 70062,2303 

PO Box 210967, Bedford, Texas 76095. USA Tel: +1 (817) 577-0165 
Gompusarve: 73700,2545 

Novator Consulting Group AB, Sv&rdvágen 11C, S-182 33 Danderyd 
Sweden Tel: +48 8 622 83 50 Fax: «46 8 622 63 51 CServe: 100341 404 
Suas APL User Group, CH-3001, Bem 1, Switzerland 

Emak: іб тато ch 

PO Box 1511, Macquarie Centre, NSW 2112, Ausbralla Tel:+01 2 257 6913 
4649 Wilens Avo., Woodland Hils, CA 91964, USA 

тө. +1 (818) 222-2759 Fax: +1 (818) 222.9250 

РО Box 85, Adelaide St Post Office, Toronto Ontario MSC 2/6, Canada. 
Email benbestQlo org 

Tour Neptune, Cedex 20, 92088 Paris la Datense 1, France. 

Tol: +33 (1) 47-78-78-00. Fax: +33 (1) 40-90-04-11 

78 Victoria Ва, Whltohavon, Cumbria, СА28 BJD UK. Tot: 01946-692588 
Dept ot Statstica, University of Warwick, Coventry, CVA TAL. UK 

Tel 01203-523360 


29 Ketchixook Lane, Ellington CT 08029, USA. Tel: +1 (203) 872-7006 
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Dyalog APL for Windows 95 


first impressions, by Adrian Smith 


Disclaimer 
Dyalog for Windows 95 arrived at the weekend, and is still not fully documented, 


so I cannot claim that this is in any way a full and fair evaluation. I have spent 
most of an afternoon playing with it - that is all. 


Background 


This is a full 32-bit native Windows program, and uses the Win32 dynamic link 
libraries, rather than the old 16-bit code. On the surface there is little visible 
change, but it loads quickly and offers easy access to all the new controls that 
form part of the ‘95’ look and feel. The only conversion required is to ONA calls 
which must be changed to use the new 32-bit libraries. Here you may need to do 
some research, for example I have yet to find the equivalent call to SndPlaySnd 
(in mnsystem.dll) which I have been using to annoy Duncan by making silly 
noises from within APL. 


A major step towards Windows 95 integration is the addition of a DropFiles 
event to Forms and many of the ‘List box’ objects. The basic interface of Windows 
95 is very much more Mac-ish, and users will expect to scatter directories (sorry, 
folders) all over the screen (sorry, Desktop) and be able to drop them on to 
application short-cuts or running applications. If you set the AcceptFiles 
property of your form to 1, it reports the name of the dropped file for you, which 
is just what you need to go and open it (having checked the extension to ensure it 
is one of yours!) without the user ever touching a menu or toolbar button. 


The rest of this short review concentrates on the new Windows controls, hopefully 
giving you some clue as to where you might find them useful. 


Some of the New Controls 


Here are a few of the things I was able to turn over this afternoon. All of them 
worked first time, and I only blew the interpreter out of the water once. Even 
here, I can report a major improvement in behaviour - it shut down without 
causing any damage to Windows, and I was able to restart Dyalog immediately. 
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The Tree Viewer 


This supersedes the msoutlin 
control from Visual Basic, but is 
used in a very similar way to 
build indented lists which the 
user can roll up and uncurl as 
required. As you can see, it is 
very easy indeed to put pictures 
into it, and you can also have 
more than one item at the top 
level (an annoying limitation of 
the old tree). I tried it with some very big lists (easily enough to blast the old .VBX 
to kingdom come) and some ludicrous indents (up to 1000, after which the control 
has had enough and keels over, but fair enough I think). 


In summary, it looks good, is dead easy to use and I can think of lots of places I 
could put it to work. 1/10 to Bill Gates for this one. 


The List Viewer 


E] Here is what it took to make the 
pictures (also used in the tree 
shown above) and then create a 
window to manage them as a 
list: 


'il' Dwe 'ImageList! (32 32) 

'11.* Owe 'Icon' 'picsVou.so* 

'11.* Омс 'Icon' 'pics\ou-whse' 

'11.* Owe 'Icon' !picsVou-plant' 

'11.* Dwe 'Icon' 'picsVou-co! 

'11.* Dwe 'Icon' 'picsVou-dv' 

"Ef. 1v"OWC'ListView'('Posn' 16 16)(120 400) 
'£f.lv'ÜWS'Items' ('Sales' ‘Warehouse’ ...) 
£f.lv'ÜWS C Imagelist! '11')('ImageIndex' 0 1 2 3 ч) 
"££. 1v'OWS' View! ‘List! 


This strings a bunch of icons into the new Image ist object and then sets this 
list as the property of the ListView. You can re-use the icons as many times as 
you like, as you specify which picture to associate with each item by a separate 
Image Index property. 1 was in 0=0 here, by the way. 
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Once you have your list, you can get all the different views that Windows itself 
provides from Explorer. If you choose to tabulate the items, the object will let you 
set the headings, column widths and so on. The user can also drag the items 
around, and (if you set AutoArrange to 1) the object will line them all up 
nicely in the new sequence. 


In summary, some imagination will be needed to make good use of this one, not 
to mention an artistic streak to make decent icons (both 32 by 32 and 16 by 16 
sizes required). Let's give ВШ 5/1 for effort. 


The TrackBar 


Two basic styles, either of which 
can be used horizontally or 
vertically (by setting е 
Hscroll and Yscroll 
properties - yuk) to set the 
volume on your sound card, 
record the results of a 
questionnaire (the wishy-washy qualitative kind), play with the fitness factors for 
your Genetic Algorithms, and so on. The range and tick values can all be set, and 
the range is no longer limited to values from 1 upward, mirabile dictu. 


The major improvement here (apart from the visuals) is that it gives you a stream 
of events as it tracks, so you can steer the aeroplane, zoom the view, or rotate the 
3D barchart as the user pulls the slider, rather than waiting for him to drop it. 
PLUS III already had this on scroll bars - but this is better for the job it does. In 
summary, another winner. 10/19 again. 


The Rich Edit Field 


This is going to take some 
The cat sat on the mat (again) Hare is some realy getting used tof It is more or 
Interesting rich text in a nch edit held less a complete word pro- 
cessor in a box, with font 
and formatting properties 
which apply to the currently 
selected text, rather than the 
whole content of the object. It will read and display existing .RTF files (such as 
those saved by Word) or you can save text (such as the above) to file, and read it 
in Word with no problems. More interesting from a programming angle - you can 
get at the Rt f Text property directly. This is always a simple vector and has all 


We get things in alics too, I seems 


Adnan 


48 


VECTOR Vol.12 No.2 


the formatting codes embedded. I can see that for applications which must build 
and print form letters (with the punter’s name in italics and the special offer in a 
fancy bold font), this could be a god-send. I think I had better reserve judgement 
on this one, but 5/10 or better is likely. Approach it cautiously, with manual in 
hand. 


Property Page and Property Sheet 


This one is lovely! Having tangled with the tabbed subforms in the current 
Dyalog release, it is wonderful to throw out all that horrible stuff with TabBtns 
(which never quite fitted right as soon as you changed the defaults) and just let 
Windows do it for you. 


The Property Sheet is a top-level object, which comes ready supplied with buttons 
for OK, Cancel, Apply and (optionally) Help. As you add Property Pages to it, 
they appear with their titles in the tabs, correctly aligned and ready to use. You 
then simply add objects to each page, and everything else is done for you. 


'ps' Owe 'PropertySheet' ‘Properties! (0 0)(120 400) 
(coord! 'pixel') 

'ps.ppi' Owe 'PropertyPage' 'Group One! 

'ps.pp2! Dwc 'PropertyPage' ‘Group Two! 

'ps.pp3! Dwc 'PropertyPage' ‘Group Three' 

'ps.pps! Dwc 'PropertyPage' 'Miscellaneous' 

'ps.pp3.ck'ÜWC Button! 'Oheck Me! (32 32) (2% 120)'Check' 


mm NENNEN] if 1 could give ВШ 1/10, 1 would. A lot of 
‘Sra Ona Í Беър Two Group Time |Мьсебгеон | messy Causeway code goes out of the 
window, replaced by a ‘Properties’ object 
реа ial which does the same job faster and better. 


Conclusions 


I have not had time for a decent go at the 
spinner or progress bar - again both replace 
existing Causeway objects but will be 
quicker to build and more responsive. If 
you are starting from raw Dyalog, the Spinner will save you a lot of time and 
some quite subtle Gui programming. 


Basically, this all came straight off CompuServe, we unzipped it and it worked. 
The new objects are excellent, but don’t expect any startling change in the rest of 
the product until later in the preview cycle. An essential upgrade for the serious 
Windows developer. 
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RECENT MEETINGS: DDE Workshop 


Presenter's Notes by Duncan Pearson 


In this article 1 will not cover all of the material from the July talk but will 
concentrate on those areas that have not already been covered in Vector articles. I 
will explain some of the terms used in talking about DDE and how they relate to 
the implementation of DDE in Dyalog APL/W and APL*Plus III. All of the 
previous Vector articles on this subject have covered DDE from and to 
Dyalog APL/W and so I will concentrate on DDE from the point of view of 
APL*Plus Ш. The examples that I use will be from the July talk, but updated 
where appropriate to use Excel 5 or above. I will show two APL workspaces 
sharing a “common” edit field, an APL workspace sharing data with an Excel 
spreadsheet and finally a small application in which a Word document can 
contain fields linked to live data in APL via an intermediate Excel spreadsheet. 


Some Terms Explained 


DDE was designed to allow applications to share their data and so the terms used 
in DDE are really just a common way of specifying what data one application 
wants to share with another. The three main elements of this specification are 
" Application", “Topic” and "Item". 
* Application. This is the name of a Windows application. In general it is the 
name of the .EXE file used to start the application, although there are 
exceptions. 


Topic. In document-centred applications this is generally the name of the 
document from which the data is to be taken. In the case of Dyalog APL/W it 
is the name of the workspace, in the case of APL*Plus it is the name of а form. 


Item. This specifies precisely what data is required. In Dyalog APL/W it is the 
name of a variable and in APL*Plus III it is the name of a label or edit field on 
a form. In Excel itis a cell range specification or a defined range name and in 
Word it is a bookmark defining a vector of text. 


Differences in Approach 


In the meanings attached by the different interpreters to the DDE terms we see 
the main difference of approach. In Dyalog APL/W it is the variables of the 
workspace that are linked with other applications. In APL*Plus III it is the text 
elements of the GUI that the APL application has built that are linked. The means 
by which the interpreter notifies the application of changes to the watched data is 
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also different. In Dyalog APL/W when a linked variable changes as a result of 
activity in another application then a DDE event occurs on the root object. This 
message does not specify which variable has changed whereas in APL'Plus Ш 
each individual GUI element that changes will receive a "Changed" event which 
can result in the execution of code in the normal manner. This is more specific and 
with more than one element linked this distinction is important. 


In the shared variable paradigm with which many APLers will be familiar, the 
processes on each side of the share are considered equal. This is not the case with 
DDE. The process requesting the share is considered the “Client” and the process 
that accepts the share and thus provides the data is considered the “Server”, 


Setting Up a DDE Client 


In APL*Plus III you link either a label or an edit field as the client of another 
application. Both the edit field and the label have specific properties where you 
can specify the application, topic and item of the data you wish to link. 


If the label is called “form1.label1” then: 


'formi.labeli'Dwi'ddeApplication' ‘Excel’ 
‘form1,label1'Qwi'ddeTopic’ '[bookilsheeti' 
'formi.labeli'Dwi'ddeItem' 'R1iCu:R3Cu! 


will define it as linked to a three-cell selection D1:D3 in sheet] of book! in Excel. 


There are two ways in which the link can be made. If the link is “cold” then the 
client application is not told of changes to the data with which it is linked. It must 
request the data. If the link is “hot” then the server application should notify the 
client application whenever the data changes. Some applications, notably Excel, 
are a little over-zealous in the notification; if any cell in a sheet changes then 
applications linked to any of the cells in that sheet are notified, even if the linked 
cells have not changed at all. This problem will crop up later. 


The following line defines the link as cold. 
'formi.labeli'Qwi'ddeMode' ‘cold! 


It is at this point that APL*Phus III actually tries to initiate the link. Only now will 
you find out whether the names specified for the application, topic and item are 
correct. If not APL will generate an error: 


OWI DDE ERROR: Unable to connect with ... 
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Finally we must request from Excel the data for the label. The DdeReguest 
method on the label does this and retums its success or failure: 


'formi.labeii'(wi'DdeReguest' 


Setting Up a DDE Server 


If APL*Plus is the server, then the topic and item requested of it must specify a 
label or edit field. Thus the topic identifies the form and the item identifies the 
label or edit field. However they are not the names of the form and label. The 
names by which they are identified are given in the "ddeTopic" property of the 
form and the "ddeltem" property of the label or edit. 


'form2'Üwi'ddeTopic' 'ServerForm' 


will set “form2” up as a DDE server. АРІ. will at this point register itself with 
Windows as a DDE server. 


'form2.editi'Üwi'ddeltem' 'ItemString' 


will give "editl" the DDE item name "ItemString". Now in a cell in Excel if you 
were to type 


=APLW|ServerForm! [temString 


the contents of the edit field would be copied into the cell, and any further 
changes to the edit field would be reflected in the cell. In Word the corresponding 
mechanism for inserting a DDE link is the DDEAUTO field. It takes as parameters 
the application, topic and server separated by spaces. It is again a hot link and 
will update automatically. Unfortunately Microsoft have removed DDEAUTO 
from the “Insert Field” dialogue box and so you have to insert the field manually. 


An Example Application 


At the BAA talk a small application was shown which demonstrates some of the 
possibilities for systems interacting via DDE. It was a combination of a Word 
template, an Excel spreadsheet and an APL workspace. The user would insert fields 
into a Word document that would be linked to data in the APL workspace. 
Associated with each expression would be a “natural language” APL expression 
that would evaluate to a statistic. The field could be moved around the document 
and formatted in any way desired. Furthermore if the data underlying the statistic 
changed then the field would be updated. 
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The Excel Spreadsheet 
This is at the heart of the system and consists of three columns and one extra cell. 


In the first column are the expressions to execute, in the second are the results and 
the third column displays the values in either the first or the second column 
depending on the value of the single cell which is Boolean (TRUE/FALSE). It is 
this third column to which the fields in the Word document are linked and this 
allows a facility in the document to display either the expressions or the results of 
them. A macro within the Word template sets the value of the single cell to TRUE 
in order to see the results and to FALSE to see the expressions. A name is defined 
for the range of the first column and one also for the second. This allows the APL 
system which is linked to the columns to manage without knowing how many 
expressions have been entered. 


The APL System 


The APL system consists of a form on which are a label and an edit field. The 
label is “hot-linked” as client to the expression column of the spreadsheet and the 
edit field is offered as a server item to which the second column of the 
spreadsheet is linked as described above. The following function sets up the form. 


Setup 
a Set up the server form to calculate the results of the expressions in an 
a Excel spreadsheet and to place them in a field accessible by Excel 


^ Initialise the record of the previous contents of the expression cells 
olddata" 


a Reset the GUI 
"si Dei'Reset" 


a Create the form 
Deselfe'form'Üvi'New' 'Form' ('vhere! 1 1 15 30) 


я Set its DDE topic. At this point APLW declares to Windows that it is а DDE 
a server. 
Dei 'ddeTopic' ‘Calculate! 


a Create a multi-line edit field and give it the DDE item name “results” 
Deself«Qui ':res.New' 'Edit' ('style! w)('where' 1 13 11 10) 
Dwi 'ddeItem' 'results' 


a Create a label and link it to the "expressions" range in the spreadsheet 
a "[booki]sheeti" 


OwselfeQwi ':expr.New' Label! ('vhere' 1 i 11 10) 
Ovi 'dde4pplication! ‘excel’ 
Owi 'ddeTopic' '[bookilsheeti' 


Dwi 'ddeltem' 'expressions! 
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a Make the link "hot" so that 4 change in excel wiil automatically reflect 
^ in the label 
Owi 'ddeMode' 'hot' 


a Set it to run the function "UpdateResults" whenever it changes 
Qui ‘onchange’ 'UpdateResults' 


a Request Excel that it provide the data to populate the iabei 
Qui 'DdeRequest' 


When the label changes (as a result of a new expression being entered or an 
existing expression amended) then the UpdateResults is run. This calculates 
the results of the expressions, formats them and places them in the edit field. 


UpdateResults:newdata 
a Calculate the results of the expressions in the caption of [wself 
^ and place them їп the field "res" of the seme form. 

a Called on a change to the label expr" 


a Get the expressions 
newdata-Duitcaption! 


a Excel sometimes notifies a change when none has occured so we must check 
tif ~newdatawolddata 
olddatarnewdata 


e Convert the newline delimited text vector to nested form 
newdatart+”(1,Otcnienewdata)Dpenclose ' ',newda 


^ Execute each subvector and place the concatenation of the results 
a in the result field "res 

Це: 'ires,text' (iseÜtcnl, “ve newdata} 
endif 


The only complication is the check to see that the data has in fact changed. Excel, 
as described above, notifies a change to all cells when any cell is changed; when 
we change the edit field and thereby the results column of the spreadsheet we are 
notified of a change in the expressions, causing us to update the edit field again ... 


The functionality of the system is provided by functions and variables in the 
workspace so that a flexible query language is implemented. 


The Word Template 
This contains the macros that manage the contents of the expressions column in 
the spreadsheet and make sure that the names defined there cover the full range 
of expressions. There is also a macro to change between "result" and “expression” 
view by (again using DDE) poking a Boolean value into the cell D1. 
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The system is shown in action below. The APL form is shown in the top left 
comer but obviously it would normally be hidden, preventing the user from 
knowing that he was using APL at all. 


Microsoft Excel - Booki vit 
=| form IL =| File Edit View Insert Format Tools 
Date Window Help ls. 
eben Е К КЛИШЕ! RET 
vesday А B c Im 
usns |43 3 Max sun 3 3 
temp 2 {day With Max sun Tuesday Tuesday 
3 |Median temp 143 143 
4 
=| Ele Edit View Insert Format Jools Table Window Help $ 


[fer ee aea EE 
The week has been cloudy with a maximum recorded sunshine of 3.4 hours on 
Tuesday and a median temperature of 14.3 degrees} 


The Future 


As you can see from Microsoft's reluctance to document the DDE hooks in their 
applications, DDE is on the way out, It has been replaced by OLE, Object Linking 
and Embedding, which provides much greater capabilities, but which has some 
strange ideas underlying it. The embedding aspect of OLE is about a document 
from one application being embedded within a document from another. It is not 
clear to me what relevance this has to the majority of APL systems. 


Central to Microsoft's thinking seems to be that users work primarily with 
documents, and this has influenced the interface to a great extent. It ignores the 
task-orientated nature of a great deal of client/server computing. How the 
Windows APL interpreters implement OLE will determine whether we are 
funnelled down this route or can use OLE to link applications in innovative ways 
undreamed of by Microsoft. 
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The Ball Clock Problem 


by Roger K.W. Hui 


This article discusses Problem #1 in the Finals of the 1995 ACM International 
Collegiate Programming Contest sponsored by Microsoft. The problem statement 
is verbatim from the WWW page http://www.acm.org/~contest/clock.html; 
earlier versions of the text appeared in the Internet news group comp.lang.apl 
(АСМ [1995], Hui [1995], and Weigang [1995]). 


The Problem 
Tempus et mobilius Tempus est mensura motus rerum mobilium. 
Time and motion Time is the measure of movement. 


— Auctoritates Aristotelis 


.. and movement has long been used to measure time. For example, the ball 
clock is a simple device which keeps track of the passing minutes by moving ball- 
bearings. Each minute, a rotating arm removes a ball bearing from the queue at 
the bottom, raises it to the top of the clock and deposits it on a track leading to 
indicators displaying minutes, five-minutes and hours. These indicators display 
the time between 1:00 and 12:59, but without "a.m." or "p.m." indicators. Thus 2 
balls in the minute indicator, 6 balls in the five-minute indicator and 5 balls in the 
hour indicator displays the time 5:32. 


Unfortunately, most commerciaily available bail clocks do not incorporate a date 
indication, although this would be simple to do with the addition of further carry 
and indicator tracks. However, all is not lost! As the balis migrate through the 
mechanism of the clock, they change their relative ordering in a predictable way. 
Careful study of these orderings will therefore yield the time elapsed since the 
clock had some specific ordering. The length of time which can be measured is 
limited because the orderings of the balls eventually begin to repeat. Your 
program must compute the time before repetition, which varies according to the 
total number of balls present. 


Operation of the Ball Clock 


Every minute, the least recently used ball is removed from the queue of balls at 
the bottom of the clock, elevated, then deposited on the minute indicator track, 
which is able to hold four balls. When a fifth ball rolls on to the minute indicator 
track, its weight causes the track to tilt. The four balls already on the track run 
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back down to join the queue of balls waiting at the bottom in reverse order of 
their original addition to the minutes track. The fifth ball, which caused the tilt, 
rolls on down to the five-minute indicator track. This track holds eleven balls. The 
twelfth ball carried over from the minutes causes the five-minute track to tilt, 
returning the eleven balls to the queue, again in reverse order of their addition. 
The twelfth ball rolls down to the hour indicator. The hour indicator also holds 
eleven balls, but has one extra fixed ball which is always present so that counting 
the balls in the hour indicator will yield an hour in the range one to twelve. The 
twelfth ball carried over from the five-minute indicator causes the hour indicator 
to tilt, returning the eleven free balls to the queue, in reverse order, before the 
twelfth ball itself also returns to the queue. 


Input 

The input defines a succession of ball clocks. Each clock operates as described 
above. The clocks differ only in the number of balls present in the queue at one 
o'clock when all the clocks start. This number is given for each clock, one per line 
and does not include the fixed ball on the hours indicator. Valid numbers are in 
the range 27 to 127. A zero signifies the end of input. 


Output 


For each clock described in the input, your program should report the number of 
balls given in the input and the number of days (24-hour periods) which elapse 
before the clock returns to its initial ordering. 


Sample Input 


30 
45 
0 


Output for the Sample Input 


30 balls cycle after 15 days. 
45 balls cycle after 378 days. 
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A Solution In J 


The balls are assumed to be labelled with the integers і. n. The clock period, the 
number of elapsed days before the clock repeats, can be computed as follows: 


m =. »0(08() 
v =. >@(18{) 
h >0(28{) 
q ›@(3&{) 


(m, {.@qu); v; 
(z: (v, {:0m); © (fim^:5) 
(z:zi(h,(:0v); (qu, retGv)) @ (f5m^:12) 
(z:z:2:Q0,retGh, {:@ћ) 9 (f1h*:12) 


qu © f12h 6 init 
*./ @ (#8>"_) OC. 
. -: € ord @ perm 


The basic data structure is a 4-element list of boxes (m; v ; h ; qu) of the balls in the 
minute, 5-minute, and hour tracks and in the queue. (The fixed ball in the hour 
track is ignored.) Verb init initializes the clock: all tracks are empty and all balls 
are in the queue. Verb f 1m models the clock action every minute; f 5m models the 
clock action every 5 minutes (including the action every minute); f 1h models the 
clock action every hour; and f 1 2h models the clock action every 12 hours. 


At the end of 12 hours, all tracks are empty and all balls are in the queue; 
therefore, the action of the clock is a permutation of the balls, computed by the 
verb perm. The order of this permutation is the LCM of the cycle lengths of the 
permutation, representing the number of 12-hour periods to return to the identity, 
and the clock period is half this number. 


The following examples illustrate the internal workings of the algorithm: 


days 45 
378 
Г ss.init 45 Initial state for 45 balls (mn: v:h:qu) 
0123985678 36 37 38 39 чо 41 42 аз 44 
fim s After 1 minute 


129656789... зе 37 зв 39 ко ы ча из se] 
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fom s After 5 minutes 
|S 6 78 9 10 1132 ... 25 vous uz ua чазаа о 
fih s After 1 hour 


Е 15 23 22 21 20 28 27 26 25 33 32 31 30 38 37 ... 


f5m^:6 s After 30 minutes 


4 9 1з 19 24 29| 30 31 32 33 34 35 36 37 38 39 ... 


fim^:2 fSm^:6 {1һ^:5 s After 5 hours and 32 minutes. 


21 2|8 26 4 5 43 Tm 36 22 1 23132 41 33 17 26 ... 


[ p =. perm 45 The queue after 12 hours for 45 balls 
25 30 0 24 35 2 44 33 15 19 13 6 29 43 B 26 40 31 


Ue Display in 5 columns 
25 30 0 2% 35 
2 44 33 15 19 
13 6 29 u3 8 
26 40 31 12 7 
38 28 3 42 32 
«1 їз 5 3? 39 
4 21 20 11 17 
34 18 27 10 23 
122 36 16 9 


C. p Cycles of this permutation 
26 1% 8 15/42 36 18 ...|u3 16 40 ...|uu 9 19 7 33 11 6 
S0. р Column display of cycles 


26 1% 8 15 


%2 36 18 12 29 39 23 


43 16 40 t 30 4 35 34 17 31 21 28 37 27 5 2 0 25 41 


uu 9 18 7 33 11 6 


#%›С. р Cycle lengths 
9727 

к/а 7277 LCM of cycle lengths 
756 
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days 45 # days is half the # of 12-hour periods 
378 
d-.days"0[27*i.101 Clock periods for 27 to 127 balls 
S2 [25 EN (2 d$! р ЖЕ ы ЕЕ; 
23 76 102 
15 85 65 138 91 
12 117 120 345 35 
37 217 EU 110 105 
378 126 6270 12 513 
1116 770 86 51 
693 180 330 559 858 
495 11067 718 455 378 
au 105 12 236 7095 
255 60 4524 3888 1530 
1955 268 714 25389 1155 
9360 2006 805 4446 1122 
540 36 570 1026 11033 
1218 6580 312 301 3367 
42780 2550 550 1365 6630 
105 861 451% 291 3960 
1464 — 1577 4284 5187 126 
17094 15330 1785 53890 25872 
5762 3325 204 3420 78120 
1228 776 273 108855 17% 


14030 80080 2415 


Permutation Power and Log 


Given the current reading (т; v ; h ; qu) of the clock, can one determine the elapsed 
time since the initial operation of the clock? (Assuming that the clock has not yet 
repeated.) 


If p is a permutation and k is an integer, the phrase q- .p&(^:k 1.#p applies p 
to the identity permutation К times, obtaining ч. By analogy with ordinary 
multiplication, q is the k-th power of p and (ord. p) |X is the log of q relative to 
р. Determining the elapsed time reduces to finding k-.p log q, where p is the 
generator permutation of the clock (the state of the queue after 12 hours) and q is 
the current permutation (the state of the queue at the next even 12-hour). We 
proceed as follows: 


First, compute the residue of q in each of the cycles of p. For example, 1 2 3 4 
5 0 consists of a single cycle and 2 3 4 5 0 1 is at 2 relative to that cycle. In 
general, the residue of q relative to a cycle c іѕ _1+ш-( (>с) {9)і.{:>с Г 
m=.#>¢; moreover, if q is a power of p, the result of q indexed by the cycle, 
(>c) (q, must be a rotation of the cycle. 
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For example: 


2345678190 
6385970124 
P 
71359028568 
[ co-.0(C. p l c1=.1{C. р 
(>с0) (9 Lerig 
1357 "683202 
{:>c0 {:>с1 
135871. 5 4689021. в 
2 2 
[ mo-.s»co [ mi-.s»ci 
4 6 
_1+mo-((>c0){q)i.{:>co _1+ш1-((2с1){9)1.{:>с1 
1 з 


The preceding computations can be encapsulated as follows: 


asserts, 


resi 
chke [: assert { -: resi |. [ 
res =, resi { chkr 
mr =. #&›@[ ‚. (reso <) 

(»c0) ces q Residue in cycle 0 

(»c1) res q Residue in cycle 1 
E 

(C. p) mr q Moduli and residues 
ud 
53 

assert 1 Return 1 if the argument is 1 
1 

assert 0 Signal error if the argument is 0 
assertion failure 

assert 0 

(C. p) mr i.-10 Not a power of p 
assertion failure 

(С.р) тг і.-10 
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The verb mr produces a 2-column table of moduli (cycle lengths) and residues. p 
log q obtains from this table by application of the GCD algorithm discovered by 
Fuclid in 300 B.C. and the Chinese Remainder Theorem by Sun Tsu in 350 A.D. 
(Graham, Knuth, and Patashnik [1988]; Iverson [1995]). The Euclidean algorithm 
produces a and b such that (m+.n)=(a+m)+b*n, and the Chinese Remainder 
Theorem specifies that d=. (n« .n) ! (ms. n)X- (r«bsn)«s«avm satisfies r-mid 
and в =п |d, for moduli m and n and residues г and s obtained from a power of p. 
1 cr is a verb such that (m,r)cr(n,s) produces (m«.n),d, then the answer 
that we seek is k- .(:cr/t. 


As indicated, the GCD is computed as a linear combination of the arguments; the 
algorithm operates by repeated remaindering. Thus: 


=@1.@2: 


{2,3 E {: = <.6%8{. 
Q.9C) @ (Gt^:( :)^:.) 9 go 
32 gcd ча GCD as coefficients. 
43 
+/ ا‎ 3 * 32 чь The actual GCD 
4 
{ a-.32 go + Initial argument for GCD 
3210 
44 0 1 
ite One iteration 
uu 0 1 
32 1 0 
it ita Two iterations. 
32 10 
iT 3f 4 
<"2 it*:(i.6) a All iterations; stop when 0=lower 
left corner 


The verb it applies to a 2-by-3 matrix: column 0 are the two current remainders; 
columns 1 and 2 are the corresponding coefficients in terms of the original 
arguments. At each step, a new remainder is computed by using row 1 as the 
divisor, and the iteration stops when the divisor is zero. 
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The version of Chinese Remainder used here rejects residues not obtainable from 


a power of p. Thus: 


| Gecri 


ab .@(gcd/ + [ X *./19C,8(.) 
cri Ce PINS SES , .8{: +/ . ab 
chkc (: assert .8(: -: ,‚&{. 

ег cri [ chkc 


a 1 cr 63 


12 9 
ШЕ 
1 
619 
3 
«1.4 
‹+1.6 


[ р*.С. со, сі 
1230567894 


[ 9=.с0 C. 1.10 
12308567288 


(C. p) mr q 
1 
60 


«1 cci 60 
12 3 


V 1 cr 6 0 
assertion failure 
1 cr 60 


Applies to (modulus, residue) pairs 
Produces (LCM, remainder) 


Verify residue 0 
Verify residue 1 

A 4-cycle 

A disjoint 6-cycle 

The perm. whose cycles are со and c1 


One application of cycle co 


Moduli and residues. 


Chinese remainder says answer is 3, 
but 1-:413 and 0-:6|3 


Chinese remainder with built-in check 


The power and log of a permutation are now defined, together with examples 


which illustrate their workings: 
pow =. 1.@#0( C.- (43»9C.6[1]) # C.Ol 
leg =. (: @ (cr/) @ (С.@[ mr )) 
pe.4 17 7 18 10 11 0 13 8 16 9 6 15 19 12 16 3 1 2 5 
р log 1.20 
P log p 
P log р{р 
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p log p{p(p 
p log p pow 3 


[ с».с. р The cycles of p 


8|15 1% 12|17 1|19 5 11 6 0 4 10 9 16 3 18 2 7 13 


ord p The order of p; LCM of the cycle lengths 
42 

p log /:p The log of the inverse is 1 less than the order 
Pr 

[ q=. p pow 38 A power of p 
19194 5 243 16 8 6 11 7 10 3 15 12 0 17 10 18 

р log 9 
зв 

c nr q Moduli and residues of q in each cycle 
1 0 
з 2 
2 0 
14 10 

cr/ e me q Order & log by repeated Chinese Remainder 
42 зв 

40 сг e mr q Select last item 
38 

к -: p log^i p pow"1 0 ks.i.ord p Verify all powers. 
1 

p log ?-#р A random perm. їз unlikely to be a power 
assertion failure (probability (ord p)%!#p) 

P log?-#p 


The verb pow exploits the dyad x C. y, which permutes y by the cycles x. 
Although the definition is less direct than p&(^:k i.#p, the time required is 
exponentially less. Thus: 


91:1 [7^5 Set random seed to 16807 
A random permutation of length 300 
The cycles of p 

#&› с Cycle lengths 


17 8176 89 121 6 
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ord p Order (LCM of cycle lengths) 
3862320 
E ke.?ord p A random power 
2908096 
(r&»c) |k The residues of k modulo the cycles 
0 2 0 0 32 16 103 4 
(p ром k) -: (i.#p) C.~ 0 2 0 0 32 16 103 usc 
1 Apply each cycle the residue # of times 
(ip) -: р pow i Works on all integer powers 


(i.#p) -: p pow 0 


(p pow j) -: p pow nsjs.?nes.ord p 


powi-.(^:(] (i. 080[)) Alternative definition pê (^:k i.#p 
timers. 6!:2 J} 2.06 under Windows on 80486/50 
timer 'q0-.p pow k' 

0.05 
timer 'qis.p powi k' 

938.35 
qo -: qt 


That is, 0.05 seconds versus approximately 15.5 minutes. Finally, to give a sense 
of the relative times required for pow and log: 


timer 'j-.p log qo’ 
0.55 


= 
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Appendix: Collected Definitions 


Clock Period 
m »9080 
Y E TETTE) 
h >@(28{) 
qu >0(з3&{) 
z 
ret 
init 
fin vihi).Gqu 
fsm )ibiqu.retGm) 9 (fim^:5) 
fih :9v);(qu,retOv)) @ (f5m^:12) 
fi2h (2:2:ziqu,ret85, {:@Һ) © (f1h^:12) 
perm qu Ә f12h Ф init 
ord *./ € (#&>"_) Ө C. 
days =. -: @ ord © perm 


Moduli (Cycle Lengths) and Residues 
assert=, 13!:80(12" )^:-. 

«:940[ - { i. (100 

f: assert ( -: resi |. [ 

resi [ chkr 

+ FEXO[ .. (res&» <) 


GCD as a Linear Combination 


go i.92: 
1с {. = {: + <.0n8{. 
ged .Je(it^:(0(.0(: M _)@go 


Chinese Remainder 
!.@(gcd/ « E% +./)®(,8{.) 
ON eti &{: * ab 
[: assert ,&{: 1 {:@cra 
cr =. ert [ chkc 


Permutation Power and Log 
pow =. 1.@8@[ C.- (w8»9C.G6[1]) # C.Ol 
log =. {1 9 (cr/) © (C.O[ mr ]) 
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Is APL a Team Sport? 


by Douglas R. Bohrer, Kemper Securities Inc. 


Abstract 


The article discusses ways in which APL programming teams can be managed to 
complete large, time-sensitive projects using parallel development. All aspects of 
the team effort are considered including design considerations, tool-kit use and 
staffing and organisation. The concepts described were used with success during 
the author's five-year consulting engagement with a large retailer. 


Do We Really Need a Team? 


Many APLers view APL as so powerful that only one or two of us are needed to 
do any project, no matter how large or complex. However, 1 think that a 
combination of complexity, lack of project definition and extreme time pressure 
can make a project impossible without a team. Let me give you some 
background on the business situation so you can understand my experience with 
APL teams. 


The projects I had to do were all merchandise planning systems with a four to six 
month development window. They all started life with a one or two page 
management directive requiring a new system or major enhancement in very 
general terms. The only clear requirement was the deadline. The main objective 
was to make it as painless as possible for company management in hundreds of 
stores, with over forty categories of merchandise, to plan sales, inventory, profits 
or losses. 


When I started, I was hired to get a creaky ADRS (A Departmental Reporting 
System) based annual planning system through one last planning cycle while all 
of the previous programmers were being transferred to other areas. The PL/1 
development of a replacement system was running late. I got through it with a 
lot of help from the folks waiting to be transferred and a lot of overtime. The job 
was just too big for me alone, or even for two people. Possibly Superman could 
have done it without help, but he wasn’t available. We had to have a team. 


The second effort was a major redevelopment with a six-person team. The PL/1 
development was going to be even later, so the current system had to be used yet 
again one more time. Our success in this effort led management to cancel the 
PL/1 project. Instead, we would develop the added functionality in the old 
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system. Once that was completed, we were given a lot more assignments. The 
APL team got as big as 21 people doing multiple projects. Personally, 1 ended up 
with a five-year extension of a two-month contract. As we did more and more, 1 
saw a pattern to our successes, which I describe below. 


Design Considerations 


The first step to each of our projects was to get the project defined. “What's it 
gotta do?” is the first question to be answered. Once the desired behaviour was 
described, the APL team had to answer "How ya gonna make it work?" 


1 found that the best way to get the users to describe what they wanted was to 
work first on what they wanted to see. We could mock up menus and menu 
paths first. Then we would design the interactive screens and the calculations 
that supported them. The next step was the reports and the calculations needed 
for them. Once we had what the user wanted to see, we could then go through 
each screen and report to ask where the data would come from. The data sources 
usually would lead us to reading files from other systems. Finally, we would ask 
what data had to be sent to other systems. 


Once there was a fairly complete definition of what the system should do, we 
had to talk about deadlines. Using the menu paths, we would define which 
options needed to be available immediately and which could wait. For example, 
reading last year's actual sales had to be ready the day the system was installed, 
but exporting the plan data to accounting systems could wait. Because of the 
time pressure, development effort had to focus on the immediate needs, with 
development of delayed features while the users were using the initial product. 


Specification Changes 

Planning systems seem to suffer from a lack of planning, possibly for the same 
reason that a shoemaker's children lack shoes. Specification changes during 
development were frequent. I found that the best way to handle them was first to 
ask often about future directions, and second to discourage changes enough to 
make sure proposed changes were worth fighting for. 


The users were reluctant to reveal what new features might be required until the 
planning methodology was approved at the highest levels. However, over time I 
was able to get them to reveal ‘hypothetical’ directions so that I could allow for 
them in the design. In exchange, I would provide ‘ball park’ estimates on how 
hard it would be to make the ‘hypothetical’ changes. The team also kept their 
mouths shut about the information. 
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Discouraging last minute changes requires a light touch. 1 always mentioned our 
motto, “We're in a hurry!". I did this so often that the users would quote it to me 
before they asked for a change. I would then tell them how much the change 
would cost, particularly in terms of time. The easy stuff I would do quickly. The 
hard stuff I would try to negotiate into something functionally adequate that was 
easier to program. In particular, 1 would try to keep changes within the 
capabilities of the existing tool-kit. 


Technical Design 


The primary design consideration is to maximise parallel development. This 
means that dependencies between functions have to be minimised by design. 
The key design concept to minimise dependencies is the hub for data flows. Just 
as airlines use a hub airport to move their passengers, the APL team should use 
pre-defined global workspace variables as a hub for all data flow. Restrict all 
data movement to movement into or out of the hub variables. In summary, to go 
to hell you have to change planes in Atlanta. 


I mean all data movement goes through the hub variables: file I/O, screen entry 
and display, reports and even data imports and exports. Make as few exceptions 
as possible. In effect, the hub variables and the utility functions that get and set 
the hub allow object-oriented programming. 


Let me illustrate by describing the data flow in a typical interactive screen. The 
screen driver program uses file I/O functions to fill the hub variables. It calls 
common calculation functions which use and set hub variables. Toolkit functions 
get data from the hub variables for special calculations using local variables. The 
screen driver then puts the results back into the hub using tool-kit functions. A 
screen-handling utility displays hub contents, and updates the hub from screen 
entries. 


Benefits of Hub Concept 


Using a hub maximises parallel development by making each piece of a project 
independent of every other piece. Each feature depends only on the hub, and is 
specified in terms of its interaction with the hub. Screen development, for 
example, can begin before there are any files because the screen driver program 
interacts with the hub. In the actual projects, screen development often started 
before there were any test files. 


Using a hub also improves reliability. Testing can be done on each piece of code 
because each is independent. Results are verified in observed effects on the hub 


69 


VECTOR Vol.12 No.2 


variables. Once verified, more code is re-usable because it works on the hub in 
the same way whenever it's used. For example, a gross profit calculation using 
and setting the hub can be called from both screens and reports because they all 
use the hub. 


Tool Kit Construction 


My basic philosophy on tool-kit functions is that they have to be easy-to-use 
black boxes. It has to be easy to learn how to use them or they won't get used. 
Over half of the people 1 worked with were short-term consultants hired 
specifically for the four to six month project window. They needed to become 
productive as fast as possible, so the tools had to be easy to learn. (Remember 
our motto.) 


I started with the functions provided by ADRS. The hub started as the global 
variable structure of ADRS and the functions ADRS used for the structure. I 
gradually replaced the file system, the screen handlers and the reporting system 
by liberally borrowing from whatever was lying around. Improvements were 
based on customer requests for new features or persistent performance 
bottlenecks. While the specifics are probably obsolete, an example is worth 
discussing to show how the evolution was eased by having the hub to work 
with. 


File 1/0 in ADRS, with VSAPL under MVS/XA, was slow because there was far 
more data than ADRS was ever designed for. In addition, allowing for multiple 
users had been accomplished by splitting the data into multiple small files which 
users could lock exclusively for update. АЙ the files were read in a loop for 
combined reports. Starting with the functions in library 2 VAPLFILE using 
AP121, I replaced the multiple ADRS files with component pseudo-files in 
VSAM clusters. File locking was based on exclusive control of QSAM "shadow" 
files, one "shadow" for each component pseudo-file. 


At the highest level, the new file functions set and used the hub variables, just 
like the ones they replaced. The only differences were in behaviour outside the 
hub. In the old file structure, multiple small files were segregated using hub 
values. In the new file structure, pieces of the component files were locked by 
using the same hub values stored on the file. None of the application code other 
than the file function calls changed at all. 


The new arrangement changed clock time for reports from 40 minutes to three 
minutes. The file functions were designed to make a switch to real component 
files relatively painless. When we finally did switch to a commercial component 
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file server, no code above the file functions changed at all. The hub variables 
were set just as before. True component files made the file-read time for the big 
reports about 15 to 30 seconds. 


The advantages of using а tool-kit are many. Reliability is enhanced because the 
code is debugged thoroughly then re-used a lot. Flexibility in addressing 
performance problems comes from modular replacement, like the file I/O 
evolution described above. Application maintenance is also easier because 
everybody knows what the tools do and can follow the code by following the 
tool use. Coding speed is increased because you only invent the wheel once. 
Training is simplified because initially people only have to learn what the tools 
do and can leave how they do it for later. 


Staffing and Organisation 

Once you have more than a few people, staffing and organisation become almost 
as important as purely technical considerations. You have to interview, train and 
structure a team of notoriously quirky APLers. You have to maintain discipline 
and morale. You also have to handle a lot more external politics because a team 
is much more visible than the usual 2-3 APL people in the back closet. 


The first step is interviewing new prospects. The most important thing here is to 
be yourself. You want the prospect to get a good idea of what working with you 
will be like so he can decide whether to take the job or not. Start by explaining 
the job. It’s a good idea to have several people gang up on the applicant so that 
the applicant gets a more complete picture of who he'll work with and what he'll 
be doing. It also gives several people in your shop an idea of whether or not they 
want to work with the applicant. I always gave a little test, asking for an 
explanation of a few lines of APL. This gave me a surprisingly accurate view of 
how well the prospect knew APL. However, I found that a real burning desire to 
work in an APL shop was probably more important than specific knowledge. 1 
found that good training expanded our choices. 


Т did most of the training personally. I did this to show the newcomer that 
training is important. My involvement allowed me to evaluate potential 
assignments for the trainees, but more importantly, allowed the trainee to get 
used to working with me. Remember that training is both social and technical. 
On their first day, everybody is a trainee even if they have lots of technical 
experience. They don't know your application and they don't know how to work 
with you. 
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At the start of training, 1 had everybody work through APL: An Introduction by 
Howard Peelle. This served as a review, since many of them had not been 
working in APL recently. While the book may be out of print, any similar 
introduction would do. The idea is to get everybody up to speed in APL without 
embarrassing anybody. Working at their own pace, trainees took from 2 days to 
2 weeks to get through the book and some selected problems. I never made a big 
deal about how long it took. I did mention often that once through the book, 
everyone was starting more or less even in terms of his exposure. 


It is important to give realistic training problems so that the trainee knows you 
are not wasting his time. The realism should include incomplete specification of 
the problems. You can do this in a fairly humorous fashion. I always used to tell 
trainees that I was going to simulate our typical users in specifying problems, so 
they had to ask me for clarification if they needed it. Then the first time I caught 
them assuming something I hadn't told them, 1 would get very wide-eyed and 
ask, “Whatever gave you that idea?”. I did this even if the assumption was OK, 
just so they would leam to ask. 


One good training exercise is to have trainees write new tool-kit functions you 
always wanted or re-write functions that are too slow. If the results are good 
enough, and they usually are, you can install the new code. For example, 1 
needed a menu generator which would take a character matrix and put it on the 
screen, making every underscore character an entry field. My boss thought this 
would be too difficult to be practical. (We were using GDDM.) I had two trainees 
do it to my specification, then installed it. It was a big boost for their morale to 
see their "graduation exercise" used by everybody on the team. It also 
demonstrated feasibility without taking any experienced people off time-critical 
projects. 


Team Structure 


When you are structuring an APL team, you have to remember that 1 man- 
month of APL code is a LOT of functionality. APL teams tend to have more 
design and management work per coder than programming teams using other, 
less productive, languages. I found that I needed a foreman for every 3-5 people. 
Foremen came in two varieties, an architect / designer who didn’t code much and 
a general contractor who coded about half-time. The architect handled most of 
the user-contact during development and designed like mad trying to keep 
ahead of the coding. The general contractor answered most programmers’ 
detailed questions during coding, contacting the architect or the users on points 
that had yet to be decided. 
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I found that people tended to develop functional, rather than project, 
specialisation. Some were really good at reports, some were super screen 
specialists and some were great at getting data from other systems. People 
would switch from project to project often, but would usually do the same area 
on each project. 


Functional specialisation seemed to be popular with the team. It also allowed 
rapid redeployment from one project to the next, which was usually a scheduling 
requirement. Team members understood their functional area completely and 
found it easy to apply their knowledge to new projects. The methods in each 
project tended to converge toward a de facto standard. As the same people 
worked on similar problem areas across multiple projects, they found it easier to 
use familiar tools in new situations. 


Discipline and Morale 


APLers are special characters... a sense of humour is required. For example, if 
everybody is writing his own tools instead of looking what's available first, you 
can’t just order everyone to look before writing. You instead suggest that the 
Committee of Public Safety and Programming Practice, fresh from its efforts in 
the recent (1789) French Revolution, is now in control of approving new tool-kit 
functions. The penalty for submitting a duplicate function could be a pain, given 
their history. Well, at least it worked for me. 


It's best to explain why something has to be done, and not just what has to be 
done. Then you can ask for advice, which you should do whenever you have 
time for it. You have to be approachable, particularly if you're the architect/ 
designer. The specification is not going to be complete unless questions are 
asked. 


Everybody has to be approachable on past mistakes. If you have a problem with 
others finding your mistakes, remember our motto. We were in a hurry when 
you messed up. 


Some things will go right, some wrong. When they go right, broadcast 
compliments loudly. When they go wrong, let the punishment fit the crime. 
Remember, you want solved problems not fear and trembling. For example, 
when I got beeped at 3 a.m. about a bug in the overnight job, I didn’t yell at the 
programmer. I just walked into his office the next morning and put the pager on 
his desk. I quietly told him the beeper was his to carry until he got a good night's 
sleep with no interruptions, Two nights later the code worked perfectly, and he 
gave me the beeper back. 
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Politics 

Most programmers I've met fear and loathe politics, probably because they don't 
really understand how politics works. For some reason, Гуе always enjoyed 
politics and have managed to do fairly well at playing them. The key thing to 
remember in the politics of programming is that programming exists to provide 
valuable information with the least possible problems. To find out what's 
valuable, you have to talk to management and your users. You also have to 
minimise their problems in dealing with you and your system. 


APL team political problems come up dealing with three types of people outside 
the team. The most immediate, if not most important, is the team's non-APL 
boss, who can barely spell APL before he gets stuck with a very strange 
collection of subordinates. The second group includes all the poor desperate 
users of your systems. They are stuck with a new or vastly altered system with 
advance billing that reminds them of the last three unnatural systems disasters 
they lived through. The last type, the most important type of all, is the client / 
systems administrator who was selected by senior management to explain their 
grand idea to the ignorant slobs in the systems department. Each type has special 
needs, which I will describe based on my results of trying a lot of things that 
didn't work before stumbling on something that did. 


Your Non-APL Boss 


This individual probably has a background in data processing management and 
is used to a normal pace of development. The speed with ‘which you do things 
will pleasantly surprise him. How you do it will make him uncomfortable. 
Always tell him what you are doing, describing the features as the users will 
view them. Avoid telling him how these things are coded in APL. Use as little 
jargon as possible so you can be sure he understands what you tell him and 
never gets the impression you are trying to baffle him with fancy-sounding 
baloney. Especially avoid APL jargon, since he is unlikely ever to have heard it 
and will begin to wonder if you speak English (or the local equivalent). 


While there are many things about APL that might require special treatment 
from the systems bureaucracy, never ask to be a special case unless there is NO 
WAY you can avoid it. Remember that special cases are special trouble for your 
manager. Save your special requests for things you will die without. 


Your Poor Desperate Users 


These people may never have seen a system that they liked using. They are used 
to rude responses from data processing every time they find a bug in the code 
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and may have given up on reporting them. They will love you if you admit your 
mistakes immediately. You should encourage error reporting. Offer to fix the 
first guy's data for free. You have to test anyway, so you might as well use his 
data. Follow up after you put out the fire and ask for suggestions. Your quick 
response will contrast starkly with what users have seen before and you will see 
а big change in their attitude. 


Using this approach | quickly became very popular. Users became almost 
deferential, saying things like, ^I know you're really busy, but I think I may have 
found something that doesn't work right." Secretaries would say things like, "T 
KNOW he'll talk to you, please wait while 1 find him." I found out that many of 
the managers I talked to routinely were notorious for not taking phone calls from 
mere peons, like my boss's boss. 


Your Client/ Systems Administrator 


[Author's note to the Pronoun Gender Police: female pronouns are used in this section 
because the majority of the people I dealt with in this position were women. If you are 
sensitive about pronouns which reflect reality, 1 apologise.] 


You are wondering how she got this job? She was drafted! She probably had a 
succession of increasingly important jobs with managerial responsibility in the 
user area. Now her success depends on you, and you look weird to her. She 
doesn't know what to expect and the lack of control scares her. The trick here is 
to be reassuring without over-promising. Never make anything look or sound 
easy. Tell her doing the job on time will be tricky, a LOT of hard work, but you 
have been successful before and expect to be this time as well. 


It is better to be early than late, so be careful when making time estimates that 
there is some room for unexpected problems. In scheduling, remember the 
specifications are going to change and make allowances. She has to publish your 
completion schedule and explain to her management when the schedule slips, so 
itis very important to make the schedule realistic. 


If you solicit rumours of changes in requirements, you can provide her with 
rough guesses about how they might change the schedule. Then she can use the 
estimates to force her management to decide if the changes are worth the delay. 
No matter how they decide, everyone involved will fee! more comfortable 
having exchanged the information. 


Encourage her to be as loud about your successes as she is about your problems. 
Constantly remind her, your motto is "We're іп a hurry!". 
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Multiprecision Arithmetic: Part II 
by John Sullivan 


The Story so Far 

In Part I, I explained my motivation for writing a multiprecision arithmetic suite 
in APL, and I gave the code for the addition, subtraction, multiplication and 
division functions, together with some other goodies, such as changing the radix 
of a number and raising a multiprecision number to a small power. I also showed 
the matrix multiplication way of generating Fibonacci numbers using the inner 
product operator with two user-defined functions. 


Now read оп... 


Input and Output 


If we set our radix to 10° and we multiply 0 11 by 0 9090910 we get 0 1 10, which 
is pretty meaningless. We need a function to convert our output to a suitable 
format so that we can read it. Similarly, we need an input function, because it is a 
lot easier to input large numbers in character format, with blanks or commas to 
group the digits in suitable blocks, than it is to input them directly as integers. 


The formatting function is called Ffmt (Float format). It takes a scalar left 
argument, which indicates the size of the blocks into which the number is 
separated. If this is positive the separator is a blank, if negative it is a comma. For 
example: 


73 Ffmt 1 12345 67890000 
12,345.678,9 

3 Ffmt ^3 12348 67890000 
12 355.678 9 

0 Ffmt ^1 12345 67890000 
12345.6789 


Since we are preparing human-readable output, the radix of our operations 
should be a power of 10. For the purposes of this function I am going to change 
the radix to 10" if the radix is not already a power of 10. 
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v ze(h)Ffnt aibixikinipis 
[2] +(0=11р+10еЉаѕе) /л1 
(21 a+(100000000,base)chbase a о pr8 


(31 ai:pelp 


We now strip off the number of radix places, and drop trailing zeros. 


гє) п«а[о] о arita 
[53 аз(-5=(0#фа)11)+а o ntes 


But, if we have an integer with a positive exponent, we put the trailing zeros 
back, or add new zeros if there weren't any. 


Ce] »(n«0)/a2 
[7] а,+по0 9 neo 


Get the sign of the number. If it is zero, then make a quick dash for the exit. 


(81 582:—(0»5ex1*ác((à40) 11) va)/A3 
te1 2+'0' o +0 


Having got the sign we now work with the absolute value of the number. If the 
number is fractional and we need leading zeros then supply them. 


[10] à3:a-la 
[11] — *(20)/a& 
[12] — ae(ni-padta 


Format each part of the number into blank-separated blocks of numeric characters 
with leading zeros. Each block is p characters long, where p is logy. of the radix. If 
we have a fraction then put the decimal point in the appropriate place (this is 
made easy by the blanks between the blocks of numbers. 


[13] а%:т+,* ',('21',®%р)ЦЕМТ a 
[18] — *(n20)/a5 
[15]  z[(ez)enepeile'.! 


Tidy up. If our number is not completely fractional drop the leading blank (line 
16). If the number starts with one or more zeros, drop them (line 17), and replace 
one zero if this leaves the decimal point in the first character (line 18). If the 
number is negative put a minus sign at the start (line 19). If we have a decimal 
point then drop any trailing zeros (lines 20 and 21). 
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{16} as:z-(C! '=2[0]),2 
£17] ze((2z'0')i1)42 

[18] ze(Ct.'ez[01)7'0!),z 
t19] ze((s«1)/!-"),.z 

[20]  -(-'.'ez)/a6 


£21] 2+(-('0'ф2}11)+7 


If no left argument (4) was supplied, or if its value would result in no change to 
the appearance of the number as we have it now, then we have finished. If A is 
minus the size of the radix then we can bypass the code that puts the blanks into 
the right places. 


[22] a6:+(o=QNC'A')/o 
[23] +(h=px1 71)/0,a8 


Get rid of all blanks. If ^ is 0 then we have finished. The base from where the 
separators (blanks or commas) are measured is the decimal point, so find that. 
Calculate a boolean vector that can be used to put blanks in the number: line 27 
does it before the decimal point and line 29 after. Line 28 handles the case where 
we have an integer with no decimal point. 


[25] шке + 
[25] +(O=keth)/0 
[26] пе21'.' 


[27]  s+ġ(lnxbe(k+1)+k}pce(kp:), 0 
[28] +(n=pz)/a7 

[29] ses,1,Cl(toz)-ne1)»b)pc 
[30] a7:2*s\z 


Drop the resulting leading and trailing blanks. If there is a blank between the 
minus sign and the number then remove that also. 


[31] ze’ 
[32]  ze(- 
1331 "Vv. 22tZ)/a8 


[34]  ze'-',24z 


If was negative then we want comma separators, and that's it. 


[35] aB:+(hz0)/0 
[36] 2((2=' ')/1р2]+',' 
v 


The input function is called Fexec (Float execute). If the number is already 
numeric we get out quickly. The same comments about the radix for Ffmt apply 
to this function also. Be careful here with the terminology: we have two meanings 
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of the word exponent: as the number that comes after the E in a number like 
1.234Е5, and as the number defining the position of the radix point in the 
multiprecision array. I have tried to avoid confusion, but I have probably failed! 


v xeFexec х;е;Ё;1;п;р;р1:р10:9;5 
a1 oC tx) e04Y)/0 
(2) -(pi0c0ciipei0ebase)/Ai 
[3] р+8 
Сај ai:pie(pti)ipelp 


It is sometimes convenient to enter numbers in scientific notation, i.e. with an E to 
indicate multiplication by some power of 10. We handle that here. It is an error if 
there is more than one E (lower-case is also acceptable) in the number. Split the 
exponent off, drop any leading plus-sign, check that the exponent is all numbers 
with an optional leading minus, make it numeric (this function will crash here if 
the exponent is just a minus sign: this is deliberate, to allow the calling function to 
handle the error). The value of the exponent is stored in e, which will be 0 if there 
was no exponent. (Note: if your APL does not have ПР you can replace it by 
10123456789"). 


[5] +(0=e++/n+xe'Ee?’)/A2 

[6] DSIGNAL(14e)/11 

[7} qeni 

[81 е(9+1)+х o xegtx 

(9) er('+'=14+e)se 

[10]  DSICNAL(-^/(ecDD) У (ре) +(1+е)є' 7 -')/11 


[11] erse 


Now we basically repeat on the main part of the number what we did to the 
exponent, but we also have to check on the decimal point if it is present (there 
must not be more than one of these), and we split the number into two at the 
decimal point. In this part of the function s is the sign of the number, i is the 
integer part and f is the fractional part. 


[12] a2:xe(xe[D, '.-7')/х 
[13] xe(sexlOle'-"')4x 
[14].  —(0:0х)/0 

[15] ie(qexi'.')tx 

[16]. — *(0spfe(qt1)9x)/A8 
[17]  QüSIGNAL(fv.-'.)/11 


Set up the result (null vector). If there is an integer part, convert it to numerics by 
spacing it out according to the size of the radix, catenate this to the result. Do the 
same with the fractional part, which may have to have added trailing zeros in 
order to ensure that the last block is the correct length. 
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[18] a3:x«e8 

[191 c(0-pi)/au 

[20] Xxex,s(6Clp1xpi)p(p*i)*po1)Ni 

[21] 56 :*(020f)/a5 

[221 [+@ o +46 

[23] a5:fe,s((pixp£)o(pei)*poi)Vf£ef,(pl-of)p'0" 


Put the number all together, multiply all but the first element (which is the 
location of the radix point) by the sign of the number. Drop leading and trailing 
zeros. 


[20] A6:-(0^.7xe(-p£),1 "1(s])xx,£)/0 
[25]  xe(-fe(0s$x)i1)sx 

[26]  fex(0)«f 

(271 Xxef,COxs0) 11 axel ex 


If the radix we used in this function is not the same as the global radix then 
recalculate the number in the global radix. 


(28] +pt0/a7 
[29] x*(base,100000000)chbase x 


Process the exponent if it is non-zero. This is done the lazy way, by generating a 
Character number of the appropriate size and recursively calling this function. 
(This is the reason why this function is in part 2: we needed to define Fmu1 first.) 


[30] 47:*(1*xe)288,0,a9 
[31] A8:xex Fmul Fexec'C.',(( 1-e)9'0'),'1! o +49 
[32] А9:х+х Fmul Fexec'i',ep'0' 

Y 


Square Roots 


There are several algorithms for extracting square roots, but they all seem to have 
disadvantages somewhere along the line. I had to make a choice of algorithm 
here, and the one 1 have chosen seems to be the fastest function around, but it 
uses multiprecision floating-point numbers and we have to be careful to prevent 
the intermediate results blowing up and giving WS FULL. IBM[3] gives an 
algorithm using long division that is only slightly slower than this version, but it 
only works with even radices (because at one point the radix is divided by 2), 
whereas this version works with all radices as defined in part 1. For yet another 
algorithm, which is slower than these two, see [2]. 


As with division there are two functions for square roots, one for integer square 
roots with remainder, and the other for floating square roots, rounded to the 
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nearest value. Again, as with division, the floating function calls the integer 
function and massages the results. 


Isqrt (Integer square root) 
First we make sure our input is multiprecision, then we check that it is not 
negative. It’s also an error in this function if the input is not an integer. 


v zeIsqrt a;bifinigir 
[1] atscalar а 
(2) ÜSIGNAL(Ov.»a)/11 


If our number is 0 then we can exit at once. 


Сз] (07. я1+а)/в1 
[51 z-(0 0)(0 0) о +0 


Use the full precision of the number. 


[5] ài:acfullint a 


If the square is small we can use APU's built-in exponentiation primitive, 
rounding the result, then go to the end for processing of the remainder, 


(6) *(3«pa)/A2 
(73 2+0,1(Ьаѕеіа) «0.5 o +48 


First, we get an accurate starting value. Rather than just looking at the first pair of 
“digits”, we try to take as much as we can of the number for the starting value. 
For example, in base 100, if we try to extract the square root of (10 73 74 18 24). 
our starting value is (3 27 68),о, rather than the Зуу we would get from just 
looking at 10,0. If we decide to work in base 10 (heaven forbid!) the starting value 
is the same, i.e. (3 27 6 8), which is still a lot better than the 3,9 from toho. And 
in both these cases the starting value is the square root, so we can exit straight 
away having saved ourselves a lot of unnecessary effort. The 0 Fadd on line 9 is 
the quickest way to convert the scalar starting value to a multiprecision number. 


[8]  a2:n«(na)l(210a) *2«1[[8z10ebase 
tei 2+(п+1+10.5хра)+0 Faddl(baseinta)e0.5 


Have we hit upon the square root straight away? 


[10] (0 Osr-a Fsub z Fmul z)/a9 
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Now we have the usual Newton-Raphson iteration: 2,,; = 2, + @-f2,))/f'@,) 


[11] a3:bez Fadd(a Fsub z Fmul z)Fdiv т Fmul 2 


If we have more than 2 radix places we should ignore the extra (because they are 
not really needed, they slow down the algorithm, and they can give rise to WS 
FULL as we iterate). 


[12] *(5CO]1« 2)/'b- 2,14 (2*b[0]1) 5" 


If the integer part of this iteration equals the integer part of the previous iteration 
we have finished, otherwise iterate again. 


[13] «(floor b)Fequal floor z)/ay 
[14] Zeb о +43 


Drop the fractional part of the result, calculate the square-root remainder, 
catenate it to the result, and exit. 


[15] au:s(O»2[0])/'zefloor z' 
[16] a8:rea Fsub z Fmul z 
(17] a9:zez г 

y 


The called function #2 оог removes the fractional part of the number (just like |): 


v 2+floor z 
Uu] *(z[0]120)/0 
2) 2+0,1+7[0]+2 
s 


Feqrt (Floating-point square root) 
The floating square root function starts off the same way as the integer version. 


¥ z-Fsqrt a;b;d;r 


til aescalar а 

сә] — QÜSICNAL(OA.»a)/11 
[3] (0v. 1+а)/а1 
[m 2+0 0 o +0 


We now pretend that our numbers are integers, and perform the integer square 
root. The number of radix places must be even, so add an extra trailing zero if 
they are not. 
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[5] &1:-(-2|d-a[0])/a2 
t61 ded-1 o a«a,0 
t7J 582:z r«Isgrt 0,1+а 


1f the remainder is greater than the square root add 1 to the root. Put the radix 
point in the correct place. 


[ai -(-r Fgt z)/as 

гә] тет Fadd 0 1 

[10] au:z(0]ez(0]«10.5xd 
v 


Raising to a Power 


In part 1 I showed how to raise a multiprecision number to a small power. In 
many algorithms we need to raise a small number to a multiprecision power, 
taking the result modulo some prime p (if we do not take the result modulo p we 
could soon get WS FULL). 


The following function takes a scalar or multiprecision number m, and raises it to 
the power of x[0] modulo x[1}. The method is a variant of the small-power 
function shown in part 1, and is described in Ribenboim[1], p.38. 


First we separate the right argument into power and modulus. Convert the power 


into a bit string. If this is null then our power was and the result is 1, so we can 
exit. 


v zem Impow xii;mod;nigir 


[1] — x modex 
(2] 2+1 
(3) — +(O=pnebinary x)/0 


If the power is 1 then the result is the number we started with. We do not reduce 
by the modulus since we assume that the programmer has ensured that the 
number starts off smaller than the modulus. 


f+] zem 
[51 +(1=pn)/0 


Square z. Starting at the left, look at the bits of the power. If the current bit is 1 we 
also multiply z by the number m. Then take the residue modulo mod. Proceed this 
way until the right of the power is reached. 
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Of course, the first bit in the power is always 1, so we can cut some processing, 
and we get the following: 


[6] 1+1 

[7] А1:2+2 Fmul z 

[8]  -(-n(i1J)/a2 

[ei z«z Fmul m 

[10] 42:2«z Imod mod 

[11] -((рл)>»1+1+1)/41 
v 


The two functions called by Impow are as follows: 


binary converts a multiprecision integer to bits by changing its base to a power 
of 2, then converting each element of the resulting number to bits. Note that we 
do not check that the number is a multiprecision integer: the function will fail 
with a DOMAIN ERROR if it is not. 


Y zebinary x 

[1] Z+(1073741824,base)chbase x 

[21 2+(211)+2+,®(30р2)т2+1+42,2[0]00 
v 


Imod retums the residue of one number modulo another. It is, of course, possible 
to extract this result by using the remainder from Idiv, but there is a lot of 
processing in that function that is wasted when you only want the remainder, so 
Imod has been coded to run much faster. 


v atx Imod bibiijigigfis 


[11 xescalar х ¢ bescalar b 
[21 ÜSIGNAL(Ov.»b(0),x(03)/11 
[3] 5+0 9 arx 

C4) +(2<p,befullint b)/ai 


(51 *b11sbsqr)/47 

[6] А1: bi+bl1]+(3+b)(2]+base 
(7) А2:+(-0 0 Fgt a)/a3 

[8] arla о se-s 


If line 9 signals an error then there is a programming error in the function. This 
hasn't happened for me yet! 


[9]  as:DSIGNAL(a Fgt|x)/11 
[10] *(-(a Fgt 0 "1)^b Fgt a)/a5 
[11] +(-sa-a Fequal 0 0)/au 

[12] aeb Fsub a 

[13] a4:+0 

[15] a5:qelg£e(al1]«(35a) [2]15base) tbi 
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[15] jeo 
[16] -(0sg)/56 
[17]  gelgfegfxbase o yet 
[181 a6:aea Fsub д Fmul(a[01«(pa)-J)*b 
[19] +42 
[20] A7:beb[1] o xefullint x 
[21] — aex(11 о xe24x 
[22] a8:aebla 
[23] +(0=px)/ag 
[29] — acbaseia,1*x 
[25] xe1+x 
[26] +as 
[27] а49:+(—-влб=а}/а10 
[28] а+Ь-а 
[29] а10:а+0,а 
Y 


A User-defined Operator 


I suggested last time that the Fibonacci calculation could go much quicker if we 
wrote our own scan operator. This is because the primitive scan operator does far 
too much work for our purposes. As an example, consider the following function: 


v zea plus b 
Uil ieiti 
[21 z-a*b 

v 


If we perform pius reduction on the first 10 numbers, we get the following: 


1+0 
а+р1ив/110 
1 

s 


But scan in the same position is horrifying: 


1+0 
а+ріиѕ\110 
1 

4s 


And this is what is happening in our Fibonacci example. Normally we do not 
mind the extra processing all that much, because it does not take too much time, 
but with the multiprecision suite it takes up far too much. 
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For these circumstances we need a new operator, like this one: 


v 2e(F scan)xiy 


1 Zecyeox 

[2] А1:з(0=рх+1+Х)/0 
гз) z,+cy+y Fox 

C4) +41 


9 


This operator only works оп vectors, and F scan 1 2 3 4 is equivalent to 


1 (1 F 2) ((1 F 2)F 3) (((1 F 2)F 3)F u) 


In particular it does not replace \ (try -\1 2 3 4 and -scan 1 2 3 а). 
However, it will speed up the operation we are doing here considerably. 


Fadd.Fmul scan 10pc2 201 1 1 2 


And there's more ... 


You now have all of the "basic" functions: the "building blocks" of the multi- 
precision suite. If you have Dyalog APL you could have keyed these into a 
namespace called mp in order to keep them out of the way; if your APL does not 
support namespaces the functions will just be available in the workspace. 


In Part III I shall introduce some applications of these functions. 
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THE RANDOM VECTOR 


Computing Clopper-Pearson Confidence 
Limits by the Hlinois Method 


by Dietrich Trenkler (University of Osnabriick) 
Email: trenkler@rols1.oec.uni-osnabrueck.de 


Introduction. 


The method of Clopper/ Pearson [2] to construct confidence bounds for a probability p is 
covered in several text books, see e.g. Bickel/Doksum [1] or Mood/Graybill/ Boes [4]. Ex- 
plicit formulas exist only for marginal cases and their graphical determination from charts 
as the ones published in Pearson/ Hartley [3] is very cumbersome and imprecise. Further- 
more, well-known approximations based on the normal distribution may lead to unreli- 
able results especially when the sample size is small. 

In this note we present programs written in APL to compute these intervals, The first one, 
called ILLINOIS, has proven to be efficient in those cases where the determination of 
a solution of an equation f(x) = 0 is desired without the use of its derivative f. As the 
construction of the Clopper-Pearson intervals amounts to solving such an equation and 
the derivative of the function involved is complicated we resort to this method. 


The Method and the Program 


Let us start with a description of our goals. A Bernoulli trial is an experiment which can 
end in exactly two outcomes F (failure) or 5 (success). It is assumed that these Bernoulli 
trials are independently performed n times (the sample size) such that the number X of 
successes follows a Binomial distribution with probability function 


P(X =x) = (era - py", x 20,42, 


where p = P(S) is the probabilty of S. It is intended to use the observation (X = x) to con- 
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struct a confidence interval jîr, Pu] = [é:(X), bu(X)] for p such that 


PSPS ё) 21-а m 


where 1 — a is a confidence level chosen in advance, usually 0.95 or 0.99. Note that БХ) 
and p, (X) have a discrete distribution as they depend on X. This means that опе сап com- 
pute n+ 1 pairs of confidence limits for a given sample size л. 

If n is large and p not too small we can use a normal approximation and compute well- 
known confidence intervals from 


with û = Х/п and zis the a/2 quantile of the standard normal distribution. This interval 
is not exact as it can happen that P(p € 7) < 1 – a for certain values of p which is especially 
true when nis small or p 2+ 0 or pas 1. On the other hand, it can be shown that inequality 
(1) holds whenever й is chosen as the solution of 


б(р|л,х о/з,о<х<п о 
where 
п) P(x <®\= (о-в © 
Analogously, p, solves 
бр) = 0/2, 0 x «n, o 


cf. Bickel/Doksum [3], pp. 180-2. An explicit solution exists in exactly two cases: If x 
wehave G(p|n,a- 1) = 1- p = 1—a/2 and consequently the confidence interval [ĝ;, 1; 
[(о/2)\/”,1]. The case x = 0 yields the confidence interval (9. A4] = 10,1 — (a/2)/^]. An 
other cases call for application of numerical methods. 

To solve the equations (2) and (4) we could resort to Newtor-Raphson iteration which 
however makes it necessary to program the derivative of the function under considera- 
tion. Instead, we use the Illinois method which is a modification of regula falsi, see Ral- 
ston/ Rabinowitz [6], pp. 338-41. Assume that / is a continuous function defined on an 
interval [xy x] such that f(x) /(xz) < 0. By the intermediate value theorem of calculus 
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there exists an xo € (14,3) such that / (ха) = 0. This value is unique if f is strictly increas- 
ing or decreasing. Regula falsi is a sequence (z;) of values defined by 


#=3,4,5,... 6» 


where y, = f(x;). Notice that x;, , iscomputed from (x, x;..1) if ууу, < 0 otherwise it is com- 
puted from (Xi, x2). This latter case can cause poor convergence and the Illinois method 
is a modification of this procedure by replacing (x;-2,y;-2) for (x;-3,Y;-3/2) in (5) when- 
ever y;y;-, > 0. This algorithm is implemented in the function Axà*4a& ILLINOIS 
AfA where Aad ++ х.х and ûf û is string describing f(x) in terms of X 4 x The 
local û eps û controls the accuracy of & ХА and can be modified as desired. 


Vàxà*àaà ILLINOIS AfA:X:byb:bepsh 
[1] gepsa+0.00001 
[2] X+Aaûl1] 9 AxA*sûfû 
[3] X*Aaû[2] ¢ AXA*-AXA,sAfA 
[4] 
C5] REGFA: 
[61 *(Aepsà» |AxA[2])/STOP 
[7] X«àa&L11-AxAC11x(-/4a4)*-/AxA 
L8] *(0«Ax&AL21xAyaà*e Af A)/ILLI 
[9] Aaà*4a4[2],X > AxA*AxAL2], aya 
[101 +REGFA 
(111 
[121 ILLI: 
[13] &aà*áa4[11,X 
[15] gxd+(0.5x4x601]), дуд 
[15] +REGFA 
(16) 
(17) 5ТОР:4хћ+-1+Х 
v 


Examples 


Example 1. (Ralston/ Rabinowitz [6], Ex. 8.3.7(a), p. 401) Compute the root of f(x) = 
созх — xe*. Using &epSa 17-10 and D+ in line ILLINOISL? 1 we get 
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DPP«10 9 Х0+0 1 ILLINOIS '(20X)-Xx*X* 
0.3156653378 
0.446728081446 
0.5338558719 
0.5167740785 
0.5177442058 
0.514177701487 
0.5177573635 
0.5177573637 

(20X0)-X0x*X0O 
4 .4450892099£715 


Example 2. Arguing for р, it is not hard to see that С(р|л,х) is strictly monotonically 
decreasing in p for 0 < x < n and that G(0|n,x) = 1 and G{1|n,x) = 0. Hence there exists 
a unique solution of (4), that is G(p|n,x) = a/2. The following expressions illustrate the 
combination (n, x) = (10,6) 


Р+0.01х0,1100 0 COEPF+(I+0,16)110 
PLOT P,C1.1) ((Po ,*I)x(1-P)9.*10-I)*.xCOEFF 


G(p|19,6) 


Equipped with the function ILLINOIS we have tool to compute Clopper-Pearson in- 
tervals. All it takes is to compute G(p| 2,2) from equation (3). This is done in the function 
© BINOM which makes use of the recursion Р(Х = 0) = (1— Р)" and 
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P(X = 


(л-х+1)р 
(1 - p) 


xP(X -x-1), 


VR+NX G BINOM P;N;X:Q 
C1) Ne-NXC1] 0 X«iWXt2] 0 Q«i-P 
C2] 2(P=1)/'R+0 9 +0" 
[3] Ret/xN(Q* S) ( Px¥+1-X)4QxX 


* 


х= 1,2,... 


E 
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Example 3. The data matrix WEIGHTS contains weights before - after (in kg) of young 
girls receiving a treatment (family therapy) for anorexia. These are taken from Hand et 
al. [3], p. 229. Notice that the unit kilogram is presumably incorrect because otherwise 
these data are more typical for girls suffering from obesity. Adrian Smith commented that 
the unit ibs seems more appropriate. Nevertheless, we use them to illustrate the applica- 
tion of the Clopper-Pearson method. Let p be the probability that the therapy is a success 
which means that it will lead to a weight increase. If we assume independence then the 
number X of successes follows a Binomial distribution with n=17 +> 1*pWEIGHTS. 
From the following expression we see that there have been 13 weight increases so that 


(X= 13). 


83.8 
83.3 
B6 

82.5 
86.7 
79.6 
76.9 
94.2 
73.4 
80.5 
81.6 


WEIGHTS ,(--/WEIGHTS),(1.1]0>-/WEIGHTS 


95.2 
94,3 
91.5 
91.9 
100.3 
76.7 
76.8 
101.6 
94.9 
75.2 
77.8 
95.5 
90.7 
92.5 
93.8 
91.7 
98 


11.4 
11 
5.5 
9.4 
13.6 
72.9 
70.1 


mMPHPHPHPHOORPHOOPHHHuA 
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To compute a 95% confidence interval for p enter 


0 1 ILLINOIS'0.975-17 12 6G BINON X" 
0.50101 

0 1 ILLINOIS'0.025-17 13 G_BIHOM Х' 
0.93189 


Hence the Clopper-Pearson interval is [0.501,0.932) and we conclude that the therapy is 
effective. 

Example 4. Suppose we want to compute all confidence intervals [û,(x), P, (x)] for x = 
0,1,2,..., for a given n. Then it is only necessary to compute the first m = 1 + [n/2] in- 
tervals (r(x), f, (x)], x = 0,1,2,...,m. Roughly speaking this is because the first intervals 
for p « 0.5 can be regarded as a reflection of those for p > 0.5 To illustrate this point we 
consider the case ¥ ++ n= 10апа 1 — а = 0.95. The variable HALF consists of the first 
m= 1+4 [n/2] = 1 [10/2] = 6 Clopper-Pearson intervals. The rest is computed as follows: 


(0, 1E) ,FULL*HALF,(1]1-e$(([0.5xN) , 2) 4 HALF 
00 0.3085 
1 0.0025286 0.44502 
2 0.02521 0.5561 
3 0.066739 0.65245 
4 0.12155 0.73762 
5 0.18709 0.81291 
6 0.26238 0.87845 
7 0.35755 0.93326 
8 0.4439 0.97479 
9 0.55498 0.99747 
0 


10 0.6915 t 


Example5. To demonstrate formula (1) suppose we take n = 30 and 1~a = 0.95. After 
having computed all 31 intervals [5;, f, ] it is relatively easy to compute the coverage prob- 
ability c(p) := P(f; < p < д.) using a program COVERAGES _BIFON (not displayed 
here). Its output can be used to plot c(p) against p. Notice that c(p) > 0.95 holds through- 
out: 
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Coverage probabilities с(р) for CP-intervals, n = 30 
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Concluding Remarks 


Using the Illinois method it is easy to compute all Clopper-Pearson intervals for the pa- 
rameter p of a Bernoulli distribution as no derivatives have to be computed. In this way 
опе avoids the tedious and inexact task of estimating the bounds from charts. For large 
sample sizes one can resort to an approximation using deMoivre-Laplace's version of the 
Central Limit Theorem. 

Notice that formulas (2) and (4) can be applied to compute lower or upper bounds for p if 
one replaces o/2 by a. For instance, for the data from Example 3 to compute a 95% lower 
bound for p enter 


0 1 ILLINOIS '0.95-17 12 C BINOM X' 
0.53945 


Finally, it should be obvious that the approach can be easily extended to other discrete 
distributions such as the Poisson or negative binomial. 
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The Incomplete Elliptic Integrals and APL 
Joseph De Kerf 
In a previous note [1], it was shown how the common mean concept [2] may be 
used to calculate the complete elliptic integrals of the first and second kind [3, 4], 
illustrated by user-defined functions in APL. In this third note of the series, it is 
shown how the concept may be used to calculate the incomplete elliptic integrals 
of the first and second kind [3, 4]. APL user-defined functions are given. 


The incomplete elliptic integrals of the first kind К(фр) and of the second kind 
Е(ф;р) may be defined as: 


e^ 16 
көп [Fat 
вел) [реште de 


with p! <1 or -1 < p< 1. Both incomplete elliptic integrals are even functions of p 


K($;-p)- К(ф;р) and E(6;-p) = Е(ф;р) 


In addition: 
K(-&p)- -K(é;p) and, as special ages KZP) = Kp) 
Е(-фр) = -Е(өр) "PP Е(П/2;р) = Ely) 


the complete elliptic integrals of the first and second kind respectively. 


The incomplete elliptic integrals may be calculated for instance by numerical 
integration or series expansion (3, 4]. Once again however, much easier to 
program — and this is especially true in APL — are the algorithms based on the 
common mean concept [3, 4]. 


Indeed, the incomplete elliptic integrals of the first kind К(фр) may be calculated 
by application of the relation: 


Ken- d 
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where c(1,4) is the common mean of 1 and 4 = 1= p 


the х; ‘s and y, 's 


and ф,=ф 


Ф) 


(x - y)sin2"e) 


22 


zl Ax, cos (2:4) + y, sin*(2'4,)) 


being the successive arithmetic and geometric means 


respectively, while calculating the common mean (lg). The sequence 6; 
converges toward 0. as rapidly as the x, 's and y, 's toward their common mean 


c(1,9). 


A user-defined function ТЕГІ, for —П/2 x ¢ < П/2, based on this procedure, with 
as left argument ф and as right argument p, may be: 


[1] 
[21 
[31 
Cul 
[51 
[6] 
[7] 
[el 


v. 


vVIEIi[D]v 
ReF IEIi P;I;A;B 


-((12|P)v(ot2)-|F) /LAB1 


+0, R+0300,5xF+00,5 


EAB1:Re(2+Ie71),(1-P*2)*0.5 
EAB2:Ae(-/R)x10Fx2+t+IeI+1 
Be2x+/Rx(2 ioFx2«1)42 


y 


F+F-(~ 304+ 


)+2*1+1 


c(x/Re(0.5x*/R) , (X/R)*0,5) /ГАВ2 


ReFS0.5x+/R 


which for ¢ = 1/6 and = 1/3, p = 0.0(0.1)1.0, and comparison tolerance 1E^ 15 
gives respectively: 


Rie(ot6)IEIi"P«-0,0.1»i10 


R2«(0:3)1EIY"P 


з 1 20 15 20 15»P,R1,[1.5]R2 


0.523598775598299 
0.523825500165390 
0.524508805294440 
0.525658228737263 
0.527290159175087 


1.0%7197551196598 
1.0408738631962163 
1.053430587029900 
1.061%89706726052 
1.073313629047138 
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5 0.529428627051906 1.089550670051885 
6 0.532106525784461 — 1.111233322932338 
7  0.535367802759971 1.140044752769332 
+8  0.539268084090886 — 1.17890229953882u 
9 0,64 3882214161571  1.23358063254523&4u 
0 0.549306144334055 1.316957896924817 


For ф= П/2, p = 0.0(0.1)0.9 and comparison tolerance 1Е-15 the user-defined 
function IEI1 gives: 


Re(o*2)IEI1"P-0,0.1xi19 
3 1 20 159Р,[1.518 
1.570796326794897 
.57%745561517 356 
-586867847454166 
-608048619930513 
-639999865864511 
1.685750354812596 
1.750753802915752 
1.845693998374723 
1.995302777664728 
2.280549138422770 


Been 


while for ф= #П/2 and p = 41, as the integral is infinite (оо), a domain error is 
reported and the user-defined function is suspended. Eventually, the error report 
may be trapped and an appropriate response may be programmed. 


As for ф= П/2 the elliptic integral of the first kind К(фр) is the complete elliptic 
integral of the first kind K(p), the results are the same as those reported for CEIL 
in (1). 


It should be noted that for ф#&П/2 and р= +1, both the limit of ¢, and the 
common mean of c(1,4) are zero, and the procedure fails. However, the integral 
may be evaluated analytically and as such, in lines 1 and 2 of the user-defined 
function ТЕТІ, the result R is explicitly set to its value Intan(@/2 + 11/4). This 
cannot be done in one line because then, for ф= – П/2 a domain error would be 
reported, whether p=+1 or not. During the loop, in lines 4 and 5, the local 
variables A and B are introduced to improve readability, while the overhead is 
rather small. 
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The incomplete elliptic integral of the second kind E(@;p) may be calculated by 
application of the relation: 


нөр) sso - rej -vsin 4.) 


or Е(Ф;р 


and 5= 2 -y) 
T-Y s - y)ein(26,.) 


pri 


where c = (1,4), the 9,'s, and the x, 's and y, 's are defined as for K = (фур). 


A user-defined function IEI2, for -H/2 < ф<П/2, based on this procedure, with 
as left argument ó and as right argument p, may be: 


vIEI?2[D]V 
VReF IEI2 Р;1;5;Т;А;В 
[1) +(4=|P)/0,Re10F 
[21 Re(14S-Te141- 1),(1-P*2)*0.5 
[3] LAB:A-(-/R)xioFx2*1«Ie-Ie1 
[4] — Be2x*/Rx(2 10F»2*I1)*2 
[5] | F«F-(7304:B):2«I*41 
[6]  SeS«(-/R*2)x2«4I 
[7] TeT«(-/R) x10Fx2*I*1 
[8] c(R/R-(Q0.5x*/R) , (x/R)*0.5) /LAB 
[9]  Re(I:2)*(2-S)xFt*/R 
Ж, 


which for ф= П/6, 11/3, 11/2, p =0.0(0.1)1.0, and comparison tolerance 1£~15 
gives respectively: 


Rie )IEI2" P-0,0.1x110 

R2«( )IEI2" P 

R3-(012)IEI2" P 

3 1 19 15 19 15 19 15*P,R1,R2,(1.5]R3 
0 0.523598775598299 1.047197551196598 1.570796326794897 
.1 0,523372224005088 1.045660219705633 1.566861942021668 
2 0.522691528560574 1.041025536968957 1,554966546242529 
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3 0,521553538774118 1,033223451406541 1.534833464923249 
4 9.519952906838045 1.022130132787699 1.5059416123600H0 
5 0.517881934859938 1.007555555144472 1.467462209339427 
6 0.515330355388322 0.989221593549378 1,418083394448724 
.7 0,512284956933147 0,966723133094528 1.355661135571956 
8 0.508729236545024 0.939454803724951 1.276349943169906 
9 0.504642686598563 0,906456986263154 1.171697052781614 
-O 0.500000000000000 0.866025403784439 1.000000000000000 


"oooooooc 


As for ф= П/2 the elliptic integral of the second kind E =(@;p) is the complete 
elliptic integral of the second kind E(p), the results are the same as those reported 
forCEI2 in[l]. 


It should be noticed that for p = +1, both the limit of ф, and the common mean 
с=(1,4) are zero, and the procedure fails. However, the integral may be 
evaluated analytically and as such, in line 1 of the user-defined function IEI2, 
the result R is explicitly set to its value ѕіп(ф). During the loop, in lines 3 and 4, 
the local variables A and B are introduced to improve readability, while the 
overhead is rather small. The local variables S and T however, cannot be omitted, 
as their values are needed in the final evaluation of E(e;p) in line 9. 


For both user-defined functions, the relative accuracy of the common means 
calculated is at least the current value of the comparison tolerance, for the results 
shown 1£~15. In practice, however, this relative accuracy is slightly better, as the 
exact value of this common mean lies between the two iterations calculated and 
the result is set to the arithmetic mean of the final values of those iterations, 


Note: the number of iterations needed and the cpu time depend on the current 
value of the comparison tolerance and оп p, but don't depend on ф. Results of a 
benchmark for 1£~15 and p = 0.1(0.2)0.9 are shown below and compared with 
those for СЕТ1 and CEI2. The cpu times are in milliseconds. 


P 01 03 05 07 09 
number of iterations | 3 3 4 ^4 5 
СЕП P 112 112 138 138 164 
(о+2)1ЕП P 347 м7 44 41 5% 
CED P 209 209 256 256 303 
(0+2) IEI2 Р 478 48 62 62 746 
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Of course, the user-defined functions IEI1 and IEI2 are considerably slower 
than СЁТ1 and CEI2 respectively. This means that, when the user in an 
application only needs the complete elliptic integrals, it is much more efficient to 
refer to the user-defined functions CEI 1 and CEI2. 


The calculations and the benchmark were done on a MicroVAX 2000, with VAX 
APL Version 3.1, under VMS Version 4.7 [5]. The default value for the comparison 
tolerance in VAX APL is 1E7 15. The number of iterations has been measured by 
globalizing the counter 1. The cpu times have been measured with the system 
function DMONI TOR. 


References 
[1] J. De Kerf, The Complete Elliptic Integrals and APL, Vector 12.1, pp 102-105, 1995 


[2] J. De Kerf, The Common Mean and APL, Vector 11.3, pp 21-23, 1995 


[3] M. Abramowitz and LA. Stegun (Eds), Handbook of Mathematica! Functions, 
Applied Mathematics Series no. 55. National Bureau of Standards, Washington 
D.C., 1964 (and subsequent revisions) 


[4] J. Spanier and К. B. Oldham, An Atlas of Functions, Hemisphere Publishing 
Corporation, New York, 1987 


VAX APL Reference Manual — Vols 1 and 2 (Software Version V.3.0), Publication 
Numbers: Al-P142D-TE and AI-GV09B-TE. Digital Equipment Corporation, 
Maynard, Mass., June 1987 


[5 


Joseph De Kerf 
Rooienberg 72 
B-2570 Duffel 
Belgium 


100 


VECTOR Vol.12 No.2 


TECHNICAL SECTION 


This section of VECTOR is aimed principally at those of our readers who already 
know APL. It will contain items to interest people with differing degrees of 


fluency in APL. 
Contents 
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Technical Correspondence 
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Hackers’ Corner: 
Gremlins, Pixels and Brownie Points 


by Ray Cannon 


I thought I might have gremlins in a Dyalog APL/W system I recently produced. 
A letter printed on A4 paper using a standard True Type (TT) font (Times New 
Roman) would be printed at varying sizes, even though the same font size was 
requested. This was not winning me any Brownie points. 


As a test I would produce a standard printout using Windows WRITE.EXE 
selecting the Times font with a size of 36 points (approximately % inch high 
characters) as a control, and then print the same output from Dyalog and 
compare the results. 


Upon investigation, this behaviour appeared to be related to the printer driver 
being used (although now I have found an example where I can produce this 
behaviour by varying the printer setup for handling TT fonts between 
downloading as bit image and downloading as graphics). 


The problem (and solution} was finally revealed after spending a couple of hours 
with John Daintree of Dyadic Systems Ltd. 


There are two (common) ways of 
specifying the height of a font: you can 
specify the “cell height" or the 
"character height". To explain the 
difference, think back to an old- 
fashioned printing press with lead type. 
To vary the gap between lines of type, a 
lead spacer is placed between them. 
This spacer is called “leading” (as in the 
metal. In the simple case, the cell 
height is equal to the character height 
plus the "leading" height. It is possible 
to produce fonts with "built in" 
leading, useful if you want to have line- 
drawing characters able to produce 
(from “Basic Typography — а Design continuous vertical lines (that is to say 


Manual" by James Craig. Published by N 
Watson-Guptill, New York, 1990) the vertical bar character extends to the 


4 


ERA 
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full extent of the character cell). 


Under Windows, True Type fonts’ point sizes specify the “character height" not 
the "cell height”. The parameter specifying the font size that is passed to the 
utility that creates a (logical) font within Windows can define either the "cell 
height" or the "character height". This is done by using a positive number of 
pixels for the "cell height" or a negative number of pixels for the "character 
height". Under Dyalog, the Size property of the Font object supports this 
protocol. (When you interrogate the Size of a Font, Dyalog always returns the 
positive "cell height") 


Armed with this knowledge, I tested a dozen or so printer drivers (there is no 
need to have the physical printer attached) creating fonts specifying the size in 
both positive and negative pixel sizes, and then reading the size returned. The 
results are shown below (Table 1). From this table I surmise that the way the 
printer driver handles TT fonts determines how it calculates the "leading" 
height. In particular when the TT font is downloaded as a bit image, the 
"leading" height is zero. 


The difference between the cell height and the character height is the internal 
leading, as demonstrated by the following diagram: 


- Character Height 


Y r«device Points2Pixels pointsize;res 
[1] а Returns the font "Size" property value required to generate 
[2) в a True Type font vith the specified point size ол the 
[3) a specified device. 


C5) a Calculate the resolution of "device" 
L6) — reses/2721device QKG'DevCaps' a pixels per nm 


(71 — res+25.4xres а pixels per inch 
8] — reseresi72 » pixels per point 
[9] — repofatsizexres а pixels for point size 


[10] ^ Make negative as TT point size із Character height 
[11] a mot Cell height. and round to integer 
[12] г+-{0.5+г 


С19] a NOTS Some devices such as screens, vork in "logical units” 


[15] a (e.g. logical inch) rather then physical units (inches) under 
(16) а MS Windows to cater for their low resolution. 
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Table1 
Driver Cell Height of font returned on 
requesting creation of font size 
(default settings) 50 pixels -50 pixels. 
HP LJ II 50 57 
HP LJ Il 50 50 
HP LJ Ш Postscript 50 59 
HP LJ ЇЇ Postscript 50 59 
HPL] IV 50 57 
HP DeskJet 120C 50 57 
Generic Postscript 50 59 
Epson FX 80 50 59 
Canon BJ10e 50 57 
Xerox 4045 50 50 
IBM ProtoPrint 50 50 
Brother HL-6V 
Download as bit image 50 50 
Download as graphics 50 57 


Note that ALL the fonts created with a size of -50 look the same size. Those 
created with a size of +50 are of varying size. 


Addendum (from Adrian Smith) 


This solves a considerable puzzle for me — in particular I have always had the 
problem with my RAIN graphics printing that the font size was not predictable 
across printers. The basic message is very simple — use -ve sizes always, and 
what you get will be what you want. Thanks Ray! 
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TECHNICAL 
CORRESPONDENCE 


APL-Fortran Calls using quadNA 
APL2/OS2 + DNA + FORTRAN 


From: Professor R.G. Selfridge September 95 


Many APLs in recent years have provided facilities to connect across to compiled 
software in some other language. I have looked at the process in APL2 on PCs, in 
particular under the latest version, APL2/OS2. 


Mainframe APL provides DNA, or ‘name association’, along with a file that must 
describe the elements in the calling sequence. Naturally there are variations to fit 
the ultimate language that will be ‘associated’. The first version of APL2 for PCs 
was constructed for DOS-based hardware, and while it had a connection across 
to compiled routines and followed many of the constructs of the mainframe 
version any resemblance was superficial. An attached support processor was 
used (with the usual calling connections). The Fortran source (as a subroutine) 
was compiled and linked, with one or more entry points. Some provided APL 
functions then took the appropriate listings and generated source for an 
assembler that was then used by a supplied PFORTPAR.EXE file to complete the 
connections. While it sounds complicated, if the instructions were followed 
closely the result usually worked, The assembler source contained data about 
each variable in the calling sequence, type, size, results returned, as might be 
expected. Usage now required connecting to the attached support processor, 
sharing the needed two variables and then sending data over and bringing back 
results. 


The arrival of APL2/2 opened up this connection, and materially simplified its 
usage. ONA was now the interface, a file contained all the data about variables in 
the calling sequence, and the actual compiled Fortran code was stored as a 
»**.DLL file. The conventions for the identifying file become very similar to that 
of the mainframe, i.e. type, dimensions of each variable in the calling sequence, 
with a flag for ‘returning results’. 
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It is, however, reasonably clear that the principal use of this connection, within 
the development group, was to connect to C or C** programs, and Fortran had 
not been used (this is not meant as any condemnation, just a statement of 
apparent fact). Fortran has some problems which can be easily overcome, once 
you have the ‘scheme of things’. Here they are. 


1. APL2/2 runs under OS2, version 2 or later. This is a 32-bit addressing system. 
As of late 1994 there was only one Fortran that was built around 32-bit 
addressing (as far as I or the APL development group could discover). IBM's 
earlier Fortran (ex Microsoft) would compile and link under OS2, but all 
Real*8 operations appear to give 0 as an answer. There are several Fortrans 
that have extenders to access higher memory, but these also probably have 
trouble in the OS2, v.2 environment (1 have not actually tested them, no 
promises). The Fortran that was used came from Watcom, Fortran7732) (for 
the rest of this article I shall just refer to it as Fortran). This Fortran has the 
ability to compile and link into an **.DLL (Dynamic Link Library) subject to a 
few caveats. 


2. Any floating point number in APL2 is based on a high precision 
representation, whether as REAL or COMPLEX. The ОМА association will 
provide for matching in these conversions. As a result most Fortran source 
that is intended for sharing will have been built around declarations of 
REAL*8 and СОМРІ.ЕХ*16. All Fortrans provide for ‘generic’ functions for 
many functions; Watcom 77 is no exception. There are two exceptions, which 
are documented if the user thinks to look for them (I did not at first). Here 
they are: 


If X is COMPLEX*16 and you want to get the real part you must use 
DREAL(X). The choice of REAL(X) will give only single precision. (Surprise, 
IMAG works as you might expect). 


If X and Y are REAL“, the double precision complex requires DCMPLX(X,Y). 
If you use CMPLX(X, Y) you will get a single precision complex. 


3. In order to use INA the calling sequences must follow OS2 conventions 
exactly. Hence there can be no communications that use registers. When 
staying entirely in a Fortran environment, Watcom allows for register 
communicating: APL2 cannot. Thus you must force the compiler to avoid 
registers in communications, and thus the compiler must be given two 
options, /BD which says a DLL is being built, and /SC which effectively says 
not to use registers in communicating. 
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4. 


There are two files that are needed for ОМА, one of which describes the 
contents of the calling sequence, the other is the **.DLL file just created. These 
must be put in the proper libraries. While many APLers may know where 
these are, or can interpret what ‘default’ really means, the documenting is ‘a 
wee tad short’ in this respect. Add to that that the OS2 ‘search and find’ utility 
first located a back-up directory when I went searching, and it is not 
surprising that nothing really worked. The calling sequence file, which might 
typically be inside an *** МАМ file, should be stored inside the BIN sub- 
directory of the directory APL2OS2. There is a sample МАМ file provided, 
but it is suggested that the user not add to that file, it could be changed in a 
later service release. The translated Fortran should be stored as an ***.DLL file 
inside the DLL sub-directory of the directory APL2OS2. 


Example 


Suppose we have a Fortran routine ROOTS(N,X,Y), which takes the complex 
coefficients X of a polynomial of degree N, and returns the solutions in Y. We 
assume the source has been sufficiently well checked that errors are remote. Here 
are the steps. 


1. 


Compile the source: 
CA wfc386 /bd /sc ROOTS 


/ bd forces a .DLL design, /sc controls registers in the calling sequences. There 
may be several ‘entry points’, but there is only one in this example (more entry 
points won't affect this command). 


Link. While this can be done with a ‘linker directive’ file, I prefer (upper case 
is from system, except for control z). 


CA wlink 

WLINK system os2v2 dll — forOS2, as a DLL file 

WLINK file roots assumes there is a file ROOTS.OBJ 

WLINK export roots entry point is roots, as many entries as desired, 
each entry should be ‘exported 

WLINK CILZ ends linker entries 


Your directory will now contain ROOTS.OBJ from the compilation (no need to 
keep), and a file ROOTS.DLL which needs to be moved to the sub-directory 
DLL of the directory APL2OS2. 
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3. Create the communications data lists. This is well described in available 
documention. For this example X is complex, of size determined on use, as is 
Y. Since Y is used to return results it must be created properly in the APL 
using routine (see later). The file, called TEST.NAM, is stored in the BIN sub- 
directory of the directory APL2OS2, and is 


ick ROOTS 
dink SYSTEM 

:lib. ROOTS 

:proc.ROOTS 

:valence 0 1 0 

:rarg(G4 1 3)04 0)016 1 50716 1 5 


The usage of > in the last line is to indicate that this argument is used for 
returning results. There is some case sensitivity, so create the Fortran routine 
with upper case for the routine name, and keep all names upper case in 
TEST.NAM. The APL program that calls this routine is now (shortened to show 
only the critical parts) 


Y +0075 X;N;ROOTS 

N e'«TEST.NAM»! 11 QNA'ROOTS' 

aComment:TEST.NAM is name of characteristic file 
Y +(Ne(-1)4tpX)p 071 

aComment: create Y of correct length and type 
ROOTS NW 'Y' 


Comment: The argument for ROOTS must be a nested list, so по“, ". 


Since results аге in Y, it is passed as a name. N is generated as a scalar with + 
since it is defined as a scalar in TEST. NAM. The results will be in the variable Y. 
I note that since the DLL file has the entry name ROOTS, it is important to make 
ROOTS local to the APL function. If enough different names are used this 
localizing can be avoided. 


The primary use of ONA that I have made is this one application, but I believe it 
should cover nearly all uses of Fortran. 1 have not attempted using a Fortran 
function routine (i.e. one that transmits a result back in-line) with DNA. 


I also offer many thanks to the APL2 development group, who offered major 
support in tracking all this down. 


Ralph Selfridge 

300 CSE 

University of Florida 

Gainesville, FL 32611, USA 

Email: selfrid@nervm.nerdc.ufledu 


108 


VECTOR Vol.12 No.2 


Calculating nCr 
From: Norman Thomson 28th August 1995 


Alan Sykes rightly makes the point in Calculating Probabilities for Elementary 
Distributions (Vector Vol.12 No.1) that anyone who is in the business of doing 
serious probability calculations should regard a function LNNCR which 
computes the logarithm of nCr as essential. The function given by Alan can be 
appreciably shortened by using “each”. First a basic function for scalar д and г 
is: 

[0] Zen inncr г 

[1] 2++/(e1+n-ır)-eırerln-r 


which can be generalised to Alan's LNNCR, where the right argument is an array 
of any shape, by: 


[0] Z«N LNNCR R 
[1] Z«(pR)pN Inncr ",R 


If only a small number of values are required, which is often likely to be the case 
in practice, this is the simplest way to go about things. 


However, as Alan points out, if the right argument is an array with a largish 
number of items, an efficiency gain is achieved by computing many nCrs, and 
then selecting the ones which are actually wanted. An extension of 2nncr to do 
this is 

[0] Zen 1nncri ritiu 


[1] Z«0,«X(oden-iu)-eisuc[/teuLn-ue,r 
[2] Z«tpr)pZ(1st1 


Alan says that extending this to array left arguments as for the “shriek” primitive 
is tricky and calls for a volunteer. I suggest: 


Сој Zen LNNCRi r 
[1] Zen Іппсгі`єг 


Unless Alan has some subtlety in mind which is escaping me, it seems I have just 
volunteered! Turning to HYPERGEOM, the function which evaluates hyper- 
geometric probabilities, I prefer the form: 

[0] Zer EYPERGEOM Wnm;¥;n,m 


C1] (W n m)-Nnm 
£2] Z-*-/(n,N-0,n) inncr^r,m-0,r 
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This is not only shorter than Alan's version, but also emphasises the essential 
symmetry of the calculation, which is not apparent in the mathematical notation, 
namely that the n,N-0,n which makes up the left argument to 2nncr refer to the 
overall sizes of population and sample, while r,m-0,r refer to the respective 
numbers of “successes”. 


Continuing Alan’s illustration of the National Lottery odds, the inverse 
probabilities of gaining 3, 4, 5, 6 successes, that is the odds expressed as "1 in...” 
are given by 


[+3 4 5 6 HYPERGEOM c49 6 6 
57 1033 54201 13983816 


Next, Poisson probabilities and cumulative probabilities are given by the 
following chains of functions: 


L0) 2+# pois MU 
Сај Zex(8/ (00) -e 1) -MU 


[0] Z-N poisi MU 
[1] Z-N pois MU 


[0] ZeN POIS MU 
[1] Ze(cN)po1si1"MU 


8 4¥>0 1 2 3 4 5 POIS 12 
«3679 — .3679  .1839  .0613  .0153 .0031 
.1853 .2707  ,2707 .180у  .0902  .0361 


[0] Z«N poiscdf MU 
Са] 2e4+/((4)+0,1-/-¥) pois" MU 


[0] Z«N poiscdfi MU 
[1] ZesNCCO 147 10148) ," N)poiscd£^ MU 


[0] ZeN POISCDF MU 
[1] Z-(cN)poiscdfi^MU 
1000 POISCDF 1000 
0.50840937 
B 92100 120 140 POISCDF 100 120 
+5266 +9773 .9999 
0347 +5243 +3669 


The functions LNNCRi, POIS and POISCDF all work with lower rank 
arguments, so that the subsidiary functions are not required at user level. 


Sooner or later, as the size of the integer arguments increases, even the above 
functions will hit limits, either doma£n errors, if even the logarithms exceed 
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the largest representable number, or ws full if the number of integers in the 
numerator and denominator of the recursive formula become too large. 


The first problem is addressed by increasing the base value of the logarithms (the 
value ¢ is, after all, arbitrary). 


The second problem is addressed by splitting the recursively defined fraction 
into blocks of a size which nearly fills the workspace, and using iteration. (On 
my computer 10000 is a desirable size.) 


When computing very large binomial coefficients, it might well be the case that 
the number of decimal digits in the result was the primary object of interest, in 
which case it is appropriate to use: 


[0] Zen Inncrio r 
[1] Z++/(1081+n-1r)-100:ırerln-r 


and the iterative function indicated above is: 


[0] Zen ler г;ё;и;і 

[1] Zen Inncri0 10000| terin-r o 1+0 
[2] Li:e(t«ue10000x1«141)/0 

Сэ] 2«2««/(10e1«n-u)-10eu*1:10000 о +11 


As an example, the number of decimal digits in 1000000 C 500000 is: 


11000000 icr 500000 
302464 


Floating-Point Precision 
From: Joseph De Kerf 23 August 1995 


In [1], it was shown how the so-called common mean concept may be used to 
calculate the complete elliptic integrals of the first and second kind. APL user- 
defined functions were given and illustrated for p = 0.0(0.1)0.9/1.0. Results were 
displayed with 10 and 15 digits after the decimal point. Calculations were done 
with VAX APL under VMS (DEC) and for comparison tolerance its default value 
1E15. 


In a companion note however, the production manager observes he found 
slightly different results when run in Dyalog/W. So I did my homework again, 
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however with comparison tolerance 1E16 (instead of 1E15). In addition I did the 
same calculations but with APL user-defined functions based on: 


1. gaussian quadrature (Newton-Cotes formulae) and 


2. series expansion (MacLaurin's series integrated, Addams-Hippisley series and 
the series dedicated to A. Caley). 


For all of these runs, carefully avoiding the accumulation of errors as far as 
possible, 1 found the same results as published in [1]. So, I suppose those results 
are correct. 


On the other hand, 1 did the same calculations on other APL systems and found, 
for all of them, some differences similar to or even greater than those observed 
by the production manager. 


1 suppose that the origin of the difference lies in the floating-point mantissa 
precision. For VAX APL under VMS this is more than 16 decimal digits (about 
16.86D), such that the 16th decimal digit may be correct. However, for all the 
other APL systems I checked, floating point representation is based on the 
ANSI/IEEE Std 754-1985 [2]. This means that for those systems the floating-point 
mantissa precision is somewhat less than 16 decimal digits (about 15.95D), such 
that the 16th decimal digit of the intermediate and a fortiori of the end results 
becomes more or less unreliable. 


This is the reason why, in some implementations such as IBM's APL2/PC and 
MicroAPL's APL.68000, the upper limit of the domain of the print precision 
variable is set to 15. On the other hand, the advantage of the standard is that the 
largest number representable is about 1.80E308, while in practice, who cares 
about the accuracy of the 16th decimal digit? 


Joseph De Kerf 
Rooienberg 72 
B-2570 Duffel 
Belgium 
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Microsoft Announces APL for the 8088 
- almost 20 years late! 


From: Adrian Smith Sept 7th 1995 


References from "Gates" by Stephen Manes and Paul Andrews 
Touchstone, 1993-4 


I came across this excessively flabby tome (mostly gossip and hearsay — not 
recommended for general reading) in Germany last year, and casually looked up 
APL in the index, as one does. To my surprise, there were several entries and 
some quite thorough references from the early microcomputer press. Here are 
the more interesting ones: 


Page 91 — "An Open Letter to Hobbyists" from early 1976 ... 


"To me, the most critical thing in the hobby market right now is the lack of 
good software courses, books and software itself. Without good software and 
an owner who understands programming, a hobby computer is wasted. Will 
quality software be written for the home market? ... 


... the fact is, no one besides us has invested a lot of money in hobby software. 
We have written 6800 BASIC, and are writing 8080 APL and 6800 APL, but 
there is very little incentive to make this software available to hobbyists." 


Pages 97-98 


Later that summer Bill headed off to Seattle to work on a pet project he had 
casually mentioned in both open letters and the MicroKid ad had described as 
"Upcoming bout: APL for 8080,” APL was the acronym for an offbeat 
programming language dubbed with terminal feyness “A Programming 
Language,” and Bill had become infatuated with it. Languages often create 
cults around them, and APL had one of the most vociferous. 


An interpreted language, APL was extremely condensed: In a couple of well- 
thought-out APL statements, you could do what would take line after line of 
code in other languages. APL used special symbols, its own goofy alphabet: 
Full of arrows and deltas and squiggles, the code appeared to be written in 
hieroglyphs or Greek or something. All this gave APL high marks among 
mathematicians and scientists who liked its elegant, powerful solutions to 
complex vector and matrix problems. APL code was almost impossible to read 
and therefore hard for a third party to understand and maintain, but it was 
popular enough to be one of two languages (BASIC was the other) available 


113 


VECTOR Vol.12 No.2 


for the IBM 5100, the first in the IBM series of small computers — some would 
later say personal computers — that eventually would lead to the IBM PC. 


A year or so before, Gates had been introduced to the language by Mike 
Courtney, a Seattle APL specialist who had read about Altair BASIC, sent for 
the manual, and couldn't believe a version that good had actually been been 
implemented on a microprocessor. A casual reference to APL in a phone call 
to Gates led to a meeting that summer in Seattle, where Courtney had fired up 
Bill's enthusiasm for the language. Over the following year, Cates delved 
deeper into APL, to the point of examining at least one version of its source 
code. 


At least one person saw Bill's affection for the language as a hopeless crush. 
Back in Seattle for a visit, Paul Allen lunched with Courtney to discuss the 
possibility of his coming work with MITS. The two had never met before, and, 
sensing Allen’s coolness, Courtney asked what the problem was. The problem, 
Allen told him, was that Courtney had Bill working on this APL thing, and it 
was a bunch of crap. Allen wanted Gates to work on something real, like 
FORTRAN or COBOL — something they could actually sell. 


But Weiland agreed with Bill. FORTRAN and COBOL seemed sort of passé; 
APL was avante-garde. Now all that was left was to write a version of it for 
the 8080. [n Seattle for the summer of 1976, Gates seemed to be making 
headway, and by August he was telling the Northwest Computer Club that it 
ought to be finished in the fall. In the fall Bill returned to Harvard with APL 
uncompleted, By January the club newsletter was asking “Whatever 
happened to MicroSoft's 8080 APL?” 


Page 124 — dated mid 1975 


And by all accounts, Gates's killer schedule and corrosive social skills kept 
him out of whatever dating scene Alberquerque had to offer. In this era, 
Microsoft was his only mistress. 


It could be a cruel one. Bill's APL continued to exist only in an incomplete 
state on yellow legal pads in a desk drawer somewhere. .... The lack of APL 
was hardly a make-or-break issue. Bill Gates really was busy — not just 
programming, but making contacts with customers. Chipmakers came calling. 
Intel wanted BASIC. National Semiconductor took BASIC and FORTRAN for 
its development systems. And COBOL was finally ready. 


Given the computing power available, and the complexity of the problem, the 
outcome was probably inevitable. However it is interesting to speculate how 
different things might have been if Bill had made it work on the 6800 back in 
1976! 
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At Work and Play with J 


The Bauer-Mengelberg Problem 
by Eugene McDonnell 


This paper discusses a combinatorial problem arising in the field of music, and 
shows the importance of the A. primitive discussed in my last column. 


The problem was told to me many years ago by Ken Iverson, who had heard it 
from Adin Falkoff, who in turn had heard it from Stephen Bauer-Mengelberg, a 
conductor / programmer who was a colleague of Ken and Adin's at IBM's 
Systems Research Institute at UN Plaza in New York City in the early 1960s. 
[Picturesque but irrelevant detail: Adin tells of asking Bauer Mengelberg how one 
of the pieces he conducted at a concert the night before had gone. The answer was 
"The first movement went only so-so, but with the second movement I floated off 
the podium."] 


The problem deals with the twelve-tone music associated with the composer 
Arnold Schoenberg. I am not a musician, so I shall only briefly describe it 
musically, and then convert it into a problem in combinatorial mathematics. 


The problem is to describe all the ways in which the twelve semitones of the 
octave can be written so that each is used exactly once, and so that each interval 
possible within the octave occurs exactly once. The Penguin book A Dictionary of 
Music, by Robert Illing (1950) gives an example of such a piece in figure (f) on 
page 297. 


pop 

Xr be © u 

> iu 3——5 

The notes begin with A natural, and then alternately rise and fall, in the sequence 
B flat, G sharp, В natural, G natural, C natural, F sharp, C sharp, F natural, D 
natural E natural and D sharp. l find it convenient to number these notes 


according to their signed distances from A natural, which I number as 0. The 
twelve notes are then seen as 
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And it simplifies things if we take these mod 12, giving 


0111210394857 6 (41 


I have found it helpful visually to write these numbers as the hours on a clock 
face (using 0 in place of 12), and to connect the hours by lines in the order given, 
that is, draw a line connecting 0 to 1, 1 to 11, 11 to 2, and so on, ending with a line 
drawn from 7 to 6. 


This clock figure makes more apparent various symmetries that reduce the 
number of permutations that need to be considered. 


If we take the first difference of (A], we get the following: 


i16 9.8. 7 & 5 À 25; 


and if we take this mod 12, we get 


11038567492 11 (51 


and it is easy to see that the list [A] is a 0-origin permutation having a first 
difference, mod 12 [B] which is a 1-origin permutation. Thus we have transformed 
the musical problem, having to do with twelve-tone rows, into the combinatorial 
problem of determining all the permutations of i. 12 having a first difference 
which is a permutation of >: i. 11. That is, we want to know how many such 
permutations there are, and what they are. To make it easier to discuss “а 
permutation having a first difference mod permutation length also a 
permutation”. TU call such an object a ‘dil’ (from Distinct Interval List). 


There are 479,001,600 permutations of i. 12, so it is a large problem to sift 
through these permutations looking for dils. For example, to load the table of all 
permutations of order 12 would take 4*12*12, or 22,992,076,800 bytes. I believe 
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that this would be impossible to load in real memory on the largest contemporary 
machine. This paper explores ways to cut it down to a more manageable size. 


I heard the problem in the early 1960s when Iverson notation was available only 
on the printed page, and worked at it by hand for several months without making 
much progress. Recently I decided to tackle it once more, beginning by studying 
the permutations of smaller order. I found that dils occur only among even length 
permutations. The order two permutations are easy: both are dils: 0 1 and 1 0, 
having an interval of 1. These can be done mentally, but it quickly becomes 
necessary to develop programming tools to aid in the exploration: 


NB. permutation table 
NB. modular first difference 
NB. distinct items? 


NB. a dil? 
dils-. dil # ] NB. all dils 
pt 3 
012 
021 
102 
120 
201 
210 
mfd 015243 
14325 
mn mfd 0 1 5 2 ¥ 3 


dilois 243 
1 


Studying the dils of order 4 give us some insight into the problem: 


ils pt 4 NB. dils of length 4 
3 


d 
1 
3 
0 
2 
1 
3 
0 
2 


vunneroo 
onnvone 
Hooven 


Some symmetries are present that will let us cut the problem down in size. Only 
permutations beginning with 0 need be considered, since the others can be 
obtained by clock face rotations: 
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o=. $09] | + NB. rotate y by x 


roe 
1ro0132 
1203 

2r0 0132 
231 


and similarly for the others. I call the dils beginning with zeros ‘basic dils’, since 
all the others can be obtained from them by rotation, or, in musical terms, by 
transposing. By looking for dils only among permutations beginning with 0, our 
order !12 problem has been reduced reduced to an order !11 problem, or 
39,916,800. Here are the basic dils of orders 4, 6, and 8: 


o 
o 
о 
о 


wear 


еоеоооеооооооооооооооооооосо 
зэ чача ча ше оз з їл (л (л (л сш مه‎ ш шю юююк кюе 
л кум مر‎ зал ما ف‎ ч @ ш о лл کر دم‎ чл шю чотыш 
ю ө тал ю کو ي پم یه م‎ хз کر ب‎ зз د‎ © юк e e шю юшл D 
م يم م ې که‎ ш ю сә ю م ي ي ډه‎ кә к ي‎ ө сэл ш с Oe ON 
FRC RE REE EE RE ER EER EC EE EEE 


Further efficiencies are possible. Notice that all of these dils not only begin with 
the constant 0, but end with a constant that is half of the order: 2, 3, and 4 for 
orders 4, 6, and 8, respectively. This means that in searching for dils we only have 
to look at those permutations beginning with 0 and ending with a constant, with 
some permutation between them. 
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The desired inner permutation is given by: 


Sis. i. -. 0: , -: NB. integers thro n-1, less 0 and -:n 
si 2 
sis 
13 
si 6 
1245 
siè 
123567 
si 10 
123456789 
si 12 
12345789 10 41 


By having to consider only inner permutations of order n-2, we have now reduced 
our problem to one of order !10, or 3,628,800. Furthermore, looking carefully again 
at the tables a4, a6, and a8 above, we see that only the first half of the basic dils 
need to be tested, since the rest can be found by clock face reflections in the y-axis. 
That is, any one of the rows in the lower half of any of these tables is obtainable 
from one of the rows in the upper half. The verb ry reflects a dil about the y-axis: 


This means that to find the dils of order 12, we have to test only - : !10, or 
1,814,400 permutations. This is a reduction from !12 by a factor of 264. 


Since we can always retrieve a dil if we know its atomic number and its length, 
we don’t need to exhibit the complete row. It suffices to obtain only its atomic 
number. For example, the dils of order 4 can be obtained using only 8 integers, 
rather than the 32 required by the display of the four atoms of each permutation 
form of the dil. We can define a verb dan to give us the dils in atomic number 
form: 


den-. (dil # A.) NB. dil atomic number 
dan pt 4 NB. atomic numbers of dils of order 4 
14 6 8 15 17 19 22 


There are two additional clock face reflective symmetries in these dils. In addition 
to the y-axis symmetry mentioned above, there are reflections possible in the x- 
axis, and in both the x and y axes. For example, the dil: 


га. 0 13 2 7 10 8 & 11 5 9 6 
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can be reflected in the x-axis by: 


rzs, [i |. wt 
rX t 
09172108114356 


and in the x-y axes by: 


rxy- 1. # | iO + ] 
rxy с 


0311 51024186976 


I haven't found a way to use these further symmetries to reduce the work 
necessary to solve the dil problem. The program I use to find the primitive dils of 
ordern is: 


pdon=. 3 : 0 
NB. argument is -item list, e.g. pdon 12 5040 0 1814500 
‘nibm'=.y. 
is length of permutation 
is size of batch (depends on memory size and n) 
is base index (usually 0 initially) 
maximum item number (usually -:'n-2) 
result, list of indices of primitive dils of order n 


NB. for example, si 8 is 1 2 3 5 6 7 
NB. for n-8, h is 4 
bem do. 


.0,.((bti.1)A. $),.h NB. provide another batch 
-z,dan t NB. append primitive dil atomic #s to z 
bz.bex NB. step base by batch size 


end. 
z 
) 


The line assigning t shows the utility of being able to specify the right argument 
to the А. primitive. On my computer, it took about 10 minutes to compute the 
dils of order 10. I don't know how long it took to do those of order 12. I started it 
going just before I went to bed, and it was ready in the morning. 


For the record, the number of dils of orders 2 through 12 are: 


order primitive dils basic dils all dils 
2 1 1 2 
4 1 2 8 
6 2 4 24 
а 12 24 192 
10 ius 288 2880 
12 1928 3856 46272 
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Here are a few nicely symmetrical dils of order 12: 


ptyi2s=.646517 3154657 4275293 5762095 7289175 9306655 
ptyi2s=. ptyi2s, 11633649 12187013 13754599 14826363 
16823821 


ptyi2s А. 1.12 


01 310 2 5118 4 9 76 
0110 8 311 59 24 76 
02 310 1 5117 + 9 86 
02 71011 3 95 а 1 86 
03 1 210 5115 8 7 96 
03 7 810 5115 2 1 96 
08 3 1 в 5112 7 9106 
ow 5 а 3 1 79 211106 
оч 911 2 1 78 5 3106 
05 110 8 9 32 4 7116 
05 8 ыу 3 1 7910 2116 


Tf you're a musician you might try playing these. They also make interesting clock 
face patterns. If you have a current version of J on your computer you can see 
them drawn using the graphics facilities available. The functions sogwin and 
sline are available if you have profilejs in the command line as advised in 
installing the system. Additional information about using the J graphics facilities 
are described in the book ‘Fractals Visualization and J by Clifford Reiter, available 
from Iverson Software, Inc. 


Here is the beginning of a sample session of visualizing dils on a clock face to 
help you get started: 


]ri2s: 12 X: 1 NB. 12th root of negative 1. 
0.965926j0.258819 
alle. ri2^2*i.12 NB. first 12 powers of this root 
Jeoords=. +.а11 NB. real & imaginary parts 
1 ° 
0.866025 0.5 
0.5 0.866025 
6.12574е 17 1 
_0.5 0.866025 
0.866025 0.5 
Li 1.22515e 16 
0.866025 ED 
.0.5 0.865025 
11.83772e 16 E! 
0.5 0.865025 
0.866025 _0.5 


scaled=. 500«1+coords NB. scale to screen coordinates 
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1000 500 
933,013 750 
750 933.013 

500 1000 

250 933.013 
86.9873 750 
o 500 
66.9873 250 
250 66.9873 

500 ° 

750 66.9873 
935.013 250 


With these defined you can create a graphics window with: 


sogwin 'scaled' 
9 sline scaled 


And display the lines for a given permutation on the clock face with 


perm-. 12 i 3+cy 0 1 11 2 103 9u 8 57 6 
p-. perm(scaled 
0 sline p 


The definitions of some of the graphics verbs needed are given below: 


sogwin =, 3 : 0 
3 3 500 500 sogvin y. : 
Xa. <.X.#2.5 
z=.'pe ',у.,';хуыһ 5.0: х),';сс g isigraph;pas ',”:2{.х 
wd z.';pcenter;pscale;pcloseok;pshow sw, showna;" 
) 


sline =. 3 : 
о 0 0 sline у. 


0"1 2 


wd ' 


wd 'gshow;' 


) 


spoly =. 3 : 0712 
wd 'gpolygon '.(,' ', 


+), igshow; ! 
sfill x. 


spoly y. 
) 
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Elegant Programming 
by Chris Burke 


Part1 &. dfh 


J defines inverses for many functions, and provides various ways of making use 
of them. A recent addition to the set of inverses in J2.06, namely the inverse to 
ng#., enables the elegant title expression, and prompted this note. We will look 
first at the title expression, and then examine how it works. 


First define: 
dfíh-. 168. @ ('0123456789ABCDEF'&i.) NB. decimal from hex 


hex-. &. dfh 


Then: 

'FEED' + hex 'B' 
FEF& 

"FE! ж hex '101' 
FFFF 


Thus, hex is an adverb which returns a verb that works in hexadecimal. 


Under 


The definition of hex uses the conjunction &. (under). Given verbs u and v where 
the inverse v * 3 


is defined, then u& . v is equivalent to у‘ u v. 
Here are some examples: 

inverse of natural log is the exponential: 

12 

inverse of reciprocal is itself: 


NB. arithmetic mean 
NB. harmonic mean 
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inverse of open is box: 


n=. twinston';(i.3 4);10 20 
# &.› n 


7|32 


$ &.>n 


7|3 {2 


each-. &.» 


|. each n 
notsniw|8 9 10 11|20 10 
45 6 7 
0 i 2:3 


inverse of the binary representation is the base-2 value: 


bitwises. &.#: 


5 +. bitwise 6 
7 

5 +. bitwise 6 
4 

5 -: bitwise 6 
3 


inverse of transpose is itself: 


+/\ 1.3% 
0 12 3 
4 6 8 10 
12 15 18 21 


Inverse 


You can access the inverse directly using ^ 
define: 
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NB. bitwise OR 
NB. bitwise AND 


NB. bitwise XOR 


NB. accumulate along columns 


NB. accumulate along rows 


:_1 (power of minus 1). For example, 
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invs. ^:_1 
inverse of add 2 is subtract 2: 


482 inv 


-|&|2 


482 inv 1 2 3 
Aor 


inverse of sum scan is first differences: 


+/\ inv 23 5 7 11 
21224 


inverse of product scan is rate of increase: 


*/N inv 23 5 7 11 
2 1.5 1.66667 1.4 1.57143 


({..}.%}:) 23 5 7 11 
2 1.5 1.66667 1.4 1.57143 


inverse of p: determines number of smaller primes: 


p: 100000 NB. 100000'th prime 
1299721 
p: inv 1299721 NB. number of primes less than 1299721 
100000 
Obverse 


You can define inverses for use with the conjunctions 2. and ^: directly, using 
the conjunction : . (obverse). The result of u :. visa verb that is equivalent to u 
with an assigned obverse v. 


In general, the term obverse is used instead of inverse, since the defined obverse 
need not be a true inverse, indeed it may be an unrelated verb. 


For example: 


X: (8:2) + (4:5) NB. square root of 2? + 5* 
5.38516 
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2+ 8. ж: 5 NB. same 
5.38516 

2 *8E.(*: 2, (^&1r2)) 5 МВ. same (inverse of + : is square root) 
5.38516 

2 *&.(*: :. (^E1r3)) 5 NB. using cube root as obverse 
3.07232 


Hex 
Now let's take a closer look at the definition of hex: 


dfhz. 1688. 6 ('0123456789ABCDEF'&i.) 
hex-. &. dfh 

The inverse of '0123u56789ABCDEF'&i, is: 
{ & '0123456789ABCDEF! 

while the inverse of 1684. is: 


16 16 ... 16 &#: 
with as many 16's as required. 


Also, the inverse of f € g is g` е f^ , hence the inverse of dfh can be 
calculated, 


For any verb f: 


f hex x <=> f &. dfh x dfh' f 
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Pat2 #- 1: = #04: 


This elegant expression uses the verb q: introduced іп J2.06, and provides an 
interesting exercise for the newcomer to J. The expression contains a hook, a fork, 
ап adverb, a conjunction, and a constant function. It is fair to say that once you 
understand this expression, then you also understand the essence of functional 
programming in J. 


First let's use the definition, as follows: 
p=. #- 1: = #@q: 


р 2+1.30 
23 5 7 11 13 17 19 23 29 31 


p 123%56+1.50 
1234057 123479 123491 123493 123499 123503 


Thus, p selects the primes in a list of positive integers. How does it work? 


The definition of p can be read from left to right as: select where 1 is the number of 
prime factors. 


Let’s build up this definition step by step. 
q: returns the prime factors of its argument, for example: 


q: 123456 
2222223 643 


The verb #@q: returns the count of the number of prime factors. Here, the 
conjunction @ (atop) creates a new verb that applies # to the result of q:. 


ро 123456 
8 


We now want to generate a boolean where a 1 indicates a prime, i.e. where the 
count of the number of prime factors is 1. This is achieved by the following fork: 


pis. 1: = ро 


pi 241.12 
110101000101 
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The 1: needs some explanation. A fork is a sequence of three verbs, f g h 
where: 


(f gh) x <=> (f x) g (hx) 


Note that the three elements of a fork are verbs. In the definition of p1, the 
leftmost verb is 1:, which is a verb that returns the value 1, given any argument. 
Here, J makes clear the difference between a number, and a verb that returns that 
number. This distinction is not made in standard mathematical notation, which 
can be confusing. 


We now use the boolean to select the primes. The selection verb is #, which takes 
à boolean left argument: 


1101021423u567 
2357 


In this case, however, we want to use # with the boolean as a right argument. To 
do so, we create a new verb with the adverb - (passive), that swaps its arguments: 


p2-. #- 


234567p2110101 
2357 


(You can read # as select, and #- as select where.) 


Finally, we define p as the hook p2 p1, giving us the original definition of p. 
A hook is a sequence of two verbs f. g, where: 


(f g)x <=> xfg*x 


Box Display 

Box display helps clarify the structure of functional expressions, and should be 
used by default as a learning aid. You graduate from the school of functional 
programming when you find that you no longer need box display to read 
functional expressions! For example: 


P 
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In box display, elements are grouped into twos and threes, as follows: 


2 elements: if the right element is an adverb, then the two elements represent 
the verb formed by applying the adverb to the argument on its left, 
otherwise, the two elements are a train (in the case of 2 verbs, a 
hook). 


3elements: if the center element is a conjunction, then the three elements 
represent the verb formed by applying the conjunction to its two 
arguments, otherwise, the three elements are a train (in the case of 3 
verbs, a fork). 


With this in mind, we can see from the above box display that p is a hook. The left 
element is the result of applying the adverb ~ to the verb #. The right element is a 
fork, whose rightmost element is the result of applying the conjunction @ to the 
verbs # and q:. 


Here are a couple more examples to illustrate: 


to=. +/ @ (1: = (+, i.)) 


Here to is the result of applying the conjunction 9, with a left argument of sum 
(+/), and a right argument of a fork, whose rightmost element is the hook 
(+. 1.). 


The expression +. i. computes the GCD (+.) of an integer and all the integers 
below it. For example, 


(+. i.) 12 
12123u1615321 


Thus, t1 defined below takes an integer argument and computes which integers 
below it are relatively prime to it (i.e. the GCD is 1): 


tis. 1: = (+, 1.) 
ti 12 
010001010001 
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Therefore to, which is defined as «/ @ t1, is Euler's totient function, i.e. the 
number of integers below a number that are relatively prime to it. 


to 12 
4 


A more efficient version is: 
t3-«. + -.0x8-.&.q: 


t3 12 


" 
Can you read its box display? 
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A Fractal Verb in] 


by Richard Oates 


At Christmas I sent a copy of J to a professor at Cornell College in Iowa where I 
went to school. I also sent this verb. It prints one line at a time like Norman 
Thomson's program in Vector Vol.11 No.3 page 17. 


Primitives like + and programs like F are verbs. F makes a square fractal. 


F 16 


Consider the sentence "She walks fast.”. The adverb “fast” changes the meaning 
of “walk”, English doesn’t need a “walkfast” verb. Consider the composition 
(+/). The Insert adverb (/) sticks the Plus verb (+) between each item of the 
argument. J doesn’t need a SUM built-in function. 


*/ 123 NB. +/ changes 1 2 3 to 1+2+3 
6 


J programs can be composed of nothing but verbs. You don't have to refer to 
application data. Iverson calls this “tacit definition”. F is a tacit program. The 
execution sequence is mapped by Track. In a track map each composed subverb 
is pegged with a horizontal line. 


F=. Kill@Draw@({.81) 
Kills. (1.0 0)"_ 

Drav=. Next@NEIPa:# 
Wext=. [ Screen 

бет ft (11:2) 2: 
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Track'F Kill Drev! 
Kill Draw 


ij 


^u 


e 
Kill@Draw сај Next@NEIP 


Track'Next Screen NEIP' 
Screen NEIP 


Tt Screen 
а — i52 ы E 


A verb applies to a noun on each side (like +) or to а noun on the right (like Р). 
An adverb applies to an argument on the left (like Insert). The argument may be 
verb or noun. A conjunction applies to an argument on each side. The result of 
an adverb or conjunction is a verb. The result of a verb is a noun. Verbs, adverbs, 
conjunctions and nouns can all be specified with the copula (= .). 


Atop (8) and Bond (8) are conjunctions. In F, Atop runs the (K111GDrav) verb 
after or "atop" the ((. £1) verb. The latter is composed of Take ({.), Bond, and 1. 
Names of primitives are constructed from one or two characters. If two, the 
second is period ({ .) or colon. This is what (t. 81) does. 


4{. 150323 НВ. No composition, tvo arguments for (. 
"(G1 NB. J supplies zeroes if necessary 


4.21 (4) NB. Composed verb ((.£1) applied to 4 


The argument of F is the argument of (€. £1). If the argument is 4 the result of 
({-81) is (1 © o 0) and (1 о о o) is the argument of (Ki11GDrav). In Draw the 
Power (^:) conjunction runs (Next@NEIP) # (Tally) times. The argument of # is 
(1 0 о 0) as is the initial argument of (Next@NEIP). 


#1000 
4 


HEIP (-:/\) makes the fractal. Not Equal (~:) is a verb. Insert (/) and Prefix (V) 
are adverbs. NEI P applies Not Equal Insert (~: /) to each prefix. Prefix structures 
its argument as a bunch of prefixes: the first item, the first two, the first three, etc. 


<\1000 NB. Вох (<) each prefix (X) 


> 


1 о]: o oT: ооо] 
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This is what J does when Not Equal Insert is applied to the last prefix (1 0 0 0) 
in place of Box. 


insert verb 
group verbs right 
07:0 is false 
07:0 is false 
1~:0 is true 


The result is also 1 when Not Equal Insert is applied to any other prefix. The 

result of NEIP on (1 0 0 0) is (1 1 1 1) In Next, Screen displays ‘HHY for 

(1 4 1 1), but Left (Г) returns (1 1 1 1) when applied to (1 1 1 1) and HH’. 
NEIP 1000 


114141 
Next@NEIP 1 0 0 0 


1111 
E NEIP 1000 


10 
"extenz1P Next@NEIP 1 0 0 0 


1010 


Next@NEIP ^: # 10 0 0 


1090 


A scalar verb like Minus (-) is applied to each item of a list. An underbar (_) 
touching a number {в a negative sign, not a verb. Scalar technique is extended to 
all verbs with the Rank (") conjunction. Rank applies a non-scalar verb to each 
part of its argument. 
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Tacit definition sharpens thought. Also, it facilitates the kind of pre-chip 
manipulation that is increasingly required for parallel and non-parallel 
hardware. The Fix (f.) adverb substitutes definitions for names. When Track is 
applied to a single name the linear definition appears below the map. 
Fractals. F f. 
Track'Fractal" 
— e 


0" 


{£ — 1: 


(1.0 0)" G(CE COE! +' (11:2) 2:))8(7:/N)5:6)8(C £1) 


The right argument of many conjunctions is executed before the left. In Fractal, 
(С — is Next and (— — 2:) is Screen. Two isolated verbs are a “hook” 
and three are a "fork". A fork applies the first and last verbs to its argument, and 
then applies the middle verb to their results. A hook applies the last verb to its 
argument, and then applies the first verb to its argument and to the result of the 
last verb. 


Vocabulary 
Verbs Adverbs 
+ / Insert 
2 \ Prefix 
i. f. Fix 
. 
t conjunctions 
‹ е Atop 
~: & Bond 
2: "Rank 
F Kap fractal п 
Xii — Display nothing 
Draw Kain Nouns 
Next — Skip ‘+! and blank 1 One 
Screen 1 0 9 Zero zero 
NEIP Infinity 
Fractal Map fractal тото 


Track Simplify boxed representation 


Two English words like A and Ad usually have nothing in common. Two J 
words like { and 4. frequently do. From selects any. Take selects from the top or 
bottom. Rank is (verb"noun). Constant is (noun"noun). Foreign links the 
operating system. Write is 11:2. Boxed representation is 51:2. Underbar ( ) is 
infinity only when isolated. 


134 


VECTOR Vol.12 No.2 


Tilting at Windmills: a New Attack 
on Nested Arrays 


by Douglas R. Bohrer 


Abstract 


This article is a proposal for a new notation for nested arrays, a notation the author feels 
is more consistent with the simple array APL concepts used before nested arrays were 
implemented. The proposal is to have ail simple atray APL functions work at the simple 
array level just as scalar functions work at the scalar level, An operator would be used 
to make simple array functions into nested array functions which would work on the 
simple array elements of nested arrays in ways similar to the way the simple artay 
functions worked on the scalar elements of simple arrays. Motivations for the proposal 
are discussed. 


Simple Arrays 

Back in the days before nested arrays, there were only objects which I will refer 
to as simple arrays. The elements of simple arrays were scalars, usually of uniform 
data type. Scalar functions were defined for arrays as working element-by- 
element through the arrays, that is at the lowest level of the array. Scalar 
expansion was defined by saying that the scalar argument would be used for each 
of the scalar elements in the array. The functions which manipulate simple 
arrays, like rho or ravel or iota, Į will call simple array functions. Simple array 
functions are defined to include operator-derived functions like inner and outer 
product. 


Nested Arrays 


I will define nested arrays as arrays in which each element is a simple array or 
scalar, Just as scalar functions work at the scalar level, I define simple array 
functions to work at the simple array level. This means that simple array 
functions work element-by-element through nested array arguments. Below, I 
will refer to this definition as Rule 1. 


A simple array argument could be expanded by saying that the simple array 
argument would be used for each of the simple array elements in the nested 
array. Below, I will refer to this definition as Rule 2. 


I need to expand the definition of nested arrays by saying that nested arrays may 
have elements that are themselves nested arrays. Rule 1, simple array functions, 
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and Rule 2, simple array expansion, would then be applied recursively until a 
simple array was reached. 


For ease in illustrating, I will display nested arrays as enclosed in braces, with 
simple arrays separated by semicolons. For example, if А,В,С and D are simple 
arrays, then (4;B) + (C;D) is defined as (A+C;B+D} according to Rule 1. 
Similarly, Rule 1 indicates that +/{4;8;С) is defined as (+/А;+/8;+/С) 
because +/ is a simple array function. Simple array expansion means that 
А,С11(В;С) is defined as ( А,[11В ; A,[1]C ) according to Rule 2. 


To illustrate the recursive use of Rules 1 and 2 I will use (4; (B; C) ) to indicate 
a nested array in which the first element is the simple array A, and the second 
element is a nested array of simple arrays B and C. Then if D and E are simple 
arrays (4; (B: C) ) -D is defined as (4-D; (B-D; C-D)) by applying Rule 2 
recursively, Similarly, (4; (B; C) ) -£Di E) is defined as (A-D; {B-E;C-E}} 
using both rules. 


Nested Array Functions 


Nested array functions will be derived from simple array functions with a left 
operator, the dollar sign. 1 think $ is a good choice because it is a widely 
available character that currently isn't being used for anything by APL. A 
derived nested array function is defined as working on nested arrays in a way 
Similar to the way the simple array function works on simple arrays. For 
example: 

{A;B} $, {C;D} is defined as (4; B; C; D) 

justas 2, 3 &is1 2 3 4. 


A little bit more complex, 


{A;B} $+.= (C;D) resolves to (A=C)+(B=D) 
justasi 24.-3 4is (1=3)+(2=4). 


It should be noted that $p (А; B;C) is defined as 3 but 
0 (4; B; C) is defined as (04; pB; oC) which is not the same. 


The definition of $ above allows making a nested array containing two nested 
array elements, that is of depth greater than 2 using enclosure, that is 


(<{А;В)) $, c{C;D) defined as ((4; B) (C; D) ). 
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A New Null Element 


One of the implications of this scheme is that there will be a null nested array 
with a nested length of 0. The notation for this would be defined so that 


(4;B) S, 0 $p X gives {A;B} asa result where X is any scalar, simple 
array or nested array. 


The depth of the null nested array is 2. It’s not clear what the simple array rho of 
the null nested array should be. Perhaps zero would be convenient. 


Simple Array Indexing 


it may be redundant but I feel compelled to start my discussion of indexing by 
reminding the reader that {A;8}[C] is defined as (A4(C]; BLC1) because 
indexing is a simple array function which according to Rule 1 applies element by 
element to a nested array. For a nested array index to a simple array Rule 1 
implies that 4((B;C)] is defined as {A(C];B(D1}. I think that similar 
methods can easily be used to resolve simple indexing in higher dimension 
objects for nested arrays. I will leave this elaboration to the student as an 
exercise. (I always wanted to say that.) It then seems logical to look at simple 
indexed assignment in similar fashion: 


{A;B}(C] Dis defined as (ALC] D;BLC] D) using both Rule 1, ће 
indexed assignment applying to the elements of the nested array, and Rule 2, 
the simple array being expanded. 


For the more complex case (4; B) [C] (D; E) is defined as (ACC 
D;B(C] E) where 4, B, C, D and E are simple arrays. 


For the trickiest case, if F is also a simple array then 

{A;B} [{C;D}1], CE; F) is defined as (ALC] , E; BCD] , F) using Rule 
1element-by-element for the nested array index as well as the other nested 
array arguments, 


Nested Array Indexing 


Nested array indexing works on nested arrays in a fashion similar to the way 
simple indexing works on simple arrays. Indexing a nested array with a simple 
array index gives nested array results whose simple array elements are 
positioned like the scalar elements would be for a simple array indexed by a 
simple array. For example, (4;B;C;D)S[1 3 1 4] is defined as 
{A;C3A;D) justas 'ABCD' [1 3 1 4] would be ' ACAD' . In a special case, 
the result of a nested array indexed by a scalar is a simple array just as a simple 
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array indexed by a scalar is a scalar. For example, {А;В; С}$ [1] is defined as А 
a simple array result, just as the result of ' ABC' [1] isa scalar 'А'. 


For a nested array index to a nested array, the most useful definition is not 
obvious. I think that the best thing to do is to define it as a nested array of the 
results from each simple array element of the index applied to the nested array. 
This definition would mean, for example, that (4:8: C;D)$[((1 3); (2 
4))) would be ((4:) B; 3). 


For the multi-dimension case, the result would be a nested array which 
exhausted every combination of simple array elements. If we define, for 
example, NN as a nested matrix then NN$L(4;B) ; (C;D)] would be 
2 28p(NN$ULA; C) ; WN$LA; D) ; NN$CB; С); NN$ LB; D1) where each of the 
indexing operations on NN would follow the method for simple array indexes of 
nested arrays. This definition is consistent with the case where if 4 is a simple 
matrix then AC1 2;3 4] is2 204(1;31,4(1;91,4[2;3],4(2;5. 


User Functions 


Should user-defined functions be simple array functions or nested array 
functions? Probably the user should have the option to define them either way. If 
a dollar sign appeared in the header just before the function name, the function 
would always be a nested array function. If called with nested array arguments, 
such a function would be passed the arguments as nested arrays. 


The default should be that all user functions would be assumed simple array 
functions. Such functions would use Rule 1 to process nested arguments. If UFUN 
is а user-defined function without a $ in the header, then UFUN {A;B) would 
Бе {ОРОН A;UFUN B) and UFUN itself would not be able to see that it had 
been called with a nested argument. 


The nesting operator should work for user functions as well as native interpreter 
functions. The operator $ placed before a user function name would feed the 
nested array argument into the function instead of calling it repeatedly using 
Rule 1. It would then be possible to define functions that could work either as 
simple array functions or nested array functions. 


Why Bother? 


I think this scheme is "more APL-ike" than current nested array 
implementations and would therefore be less confusing and easier to learn. Let 
me expand on this point. 
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T have felt for some time that the shift from simple array APL to nested array 
APL is needlessly confusing. I think the confusion comes from the difference 
between the behaviour of scalar functions on arrays and the need to disclose 
nested array elements to work on them. If scalar functions worked on simple 
arrays the way most functions work on nested arrays, then you would write 
A+”B to add all the elements of simple arrays A and B together. Unconsciously, 
the student expects that the "for each" will be assumed for nested arrays just as it 
is for scalar functions working on simple arrays. 


This confusion is NOT a figment of my imagination nor the result of a unique 
personal learning, difficulty. Since I first wrote about this in 1982, I have had the 
opportunity to talk to lots of APLers about nested arrays. Most thought learning 
them difficult. Some found them so confusing they don't use them at all. Even 
several APL fanatics at the APL91 conference admitted to me that they don't use 
nested arrays because they're too confusing to be worth the trouble. 


The scheme I propose has a lot of educational economy built into it. It builds on 
the methods of handling simple arrays to introduce nested arrays. You already 
know how nested array index or nested array catenate is going to work because 
it works just like the simple array function does. A similar method was used to 
add networking commands to UNIX. The UNIX remote copy command, rcp, 
works just like the ordinary copy command, cp, except that it works using 
somebody else's machine. Think of the leading "r" as an operator. 


In contrast with the scheme I propose, current methods have negative educational 
economy. All current methods of nested array implementation have formerly 
predictable simple array functions doing seemingly arbitrary things with nested 
arrays. These behaviours when applied to unintentional nested arrays created by 
Strand notation can get even experienced APL programmers into deep 
debugging doo-doo. Here I speak from personal experience of working with 5- to 
20-person APL teams for over 5 years. As the "debugger of last resort" I tracked 
down a lot of these problems for everyone from the beginners to APLers who 
had years’ more practical experience than 1 did. 


T think this state of affairs is unfortunate because it is limiting our ability to teach 
APL to the masses. It used to be that APL required the student to learn only as 
much as was required for the problem at hand. What he didn't know couldn't 
hurt him because he wouldn't use it accidentally. The current implementations of 
nested arrays with strand notation make what the student doesn't know 
dangerous, causing problems with objects he will be unable to identify, let alone 
fix on his own. 
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Where Do We Stand Now? 


1 have great confidence that this proposal has merit. | first wrote about it in March of 
1982. While I have refined the idea a lot since then, | have not changed the basic scheme 
much at all. 


I do not have any confidence at all that this proposal will receive any serious 
consideration from implementers, the APL Standards Committee, or anybody 
influential or powerful. The reactions of all of these people when I discussed this with 
them have not changed in the last 13 years. It is uniformly assumed that (a) 1 have an 
obvious learning disability because І can't learn to love nested arrays as they are; (b) I 
am extremely silly to suspect that people of superior intellect, such as those who 
designed nested array implementations, could be as fallible as any ordinary mortal, 
such as myself; and (с) | am naive in the extreme to expect that the greedy capitalists 
who pay the bills for implementers would ever write off the investment made in 
implementing the current design. 


Standards Deadlock Broken 


At APL91, the Standards Committee representatives jubilantly reported that they had 
reached agreement on a standard for nested arrays. All previous disagreement had 
been resolved miraculously. How did this miracle happen? Did the Lord come to all the 
members in a dream and tell them what the standard should be? Unfortunately, the 
result seems to have been far closer to the “Barebones” Parliament of Cromwellian 
England. It seems that all of the representatives who had argued in past deliberations 
against the position the committee approved were not present at the meeting where the 
proposal was adopted. Their companies were no longer interested, having left the APL 
language field. 


I'm not convinced that exhaustion necessarily means that the issue has been optimally 
settled. While I realize that it is commercially a dead issue for now, 1 hope that the 
process of “creative destruction" capitalism is famous for may eventually yield a better 
result. Until then, tilting at windmills is good exercise, 
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] Locales 


by Richard Oates 


Explicit Definition isolates names with the Local Copula (-.). Locales isolate 
names with spelling, possibly expressed, usually implied. My trip may have 
been more worst case than typical. I wrote three verbs. 


The name table bob , say, is table in locale bob, while table z is table in 
locale z. The names table bob and table z are locatives. One locale, the base, 
has no name. One locale, z, is special: its names are known in all locales unless 
overridden. Locales are usually populated by script. The application is scripted 
to the base. Secondary verbs script to named locales. If a file called add is: 


table-. 1: 0 

bys. ' '5;8,.0L .. 2 
overs, ((. ; }.)@":@, 
tbls. 1 : '(by) over х./' 
x, tbl 

) 


^. then 01:0«' Vadd' scripts the adverb to the base. This scripts the adverb to 
locale z: 


SC 2 -. 01:0 
Sc z  «'Nadd' 


I script most of my secondaries without change to locale z. If the first sentence of 
table is removed, and the last two as well, by, over and tbl become global. 
That is how I script my verb tracking program, and also the DOS editor 1 
described in Vector 11.3. The former is scripted to locale t and the latter to locale 
e, in each case to keep names isolated by locale after Explicit (:) is removed. 


I happened to have cut defined one way in the base and another way in locale z. 
Locales are populated by definition as well as script. In the next figure 51:59< 
represents the tracking program in locale t. In the figure but not for real, 51 :58« 
is also defined in locale z. The verbs trk and zing are distinct. The first box 
surprised me. I expected trk'cut' to be <;.1. Then I realized locale t knows 
nothing about the base, just as the base knows nothing about t. After trk t. is 
plugged into z it can be called trk but it still runs in t. On the other hand 
zing z runs in locale 2 and zing runs in the base. 
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Cute. «i.i 
Cut z «. «;.2 
zing 2 -. trk t «. 51:50« 

NB. plug into z 
ut'');(x.''cut  !'!');(x.''cut z !!)! 
(trk d),(zing.z. 4), :(21пд d) 


After 1 saw this I added a verb to trk that makes a locative for any name that is 
not full (cut__ from cut). 


My DOS editor in locale e scripts its result to the locale determined by the 
extension of the file name. I substituted Locale for 01:0 in the definition and 
added two verbs. For a calendar program in locale c, say, the first verb makes 
"вес." from ' Nj Vs Vcal . c! , and the second verb is almost: 


los. ('locales. '" ":@, 9) ; (1 ": 


‚ re: 01:0'°_) 
1о'всс_' 


For each Format (" :) 10 really has Do (".). Do is Шке Execute іл APL. I get a 
nonce error if вс_с_ is not global (=: }. 


T made no change to applications. І wrote three verbs that construct locatives for 
two secondaries not in т. All other locatives appear in a profile script. In an 
application session I never key a locative or see one. Names have grown up in J. 
Like the best people, they have lost their innocence but kept their simplicity. 
Locales are a major improvement. 


Richard H. Oates 

333 East 30 Street, Apt. 18E 
New York, 

NY 10016 USA 
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Index to Advertisers 
Causeway Graphical Systems Ltd 9 
Dyadic Systems Ltd 4 
MicroAPL 16 
APL Booklist (Renaissance Data Systems) 2 
Vector Back Numbers 32 


All queries regarding advertising in VECTOR should be made to Gill Smith, 
at 01439-788385, Compuserve: 100331,644. 


Submitting Material to Vector 


The Vector working group meets towards the end of the month in which Vector 
appears; we review material for issue n+1 and discuss themes for issues n*2 
onwards. Please send the text of submitted articles (with diskette as appropriate) 
to the Editor: 


Anthony Camacho, 

11 Auburn Road, Redland, 

BRISTOL, BS6 6LS 

Tel: 0117-973 0036 

Email: acamachoGcix.compulink.co.uk 


Authors wishing to use Windows Write or Word for Windows should contact 


Vector Production for a copy of the Vector APL TrueType font and Vector APL 
typebox. 


Camera-ready artwork (e.g. advertisements) and diskettes of ‘standard’ material 
(e.g. sustaining members’ news) should be sent to Vector Production, Brook 
House, Gilling East, YORK YO6 4JJ. 


Tel: 01439-788385 
Compuserve: 100331,644. 
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British APL Association: Membership Form 


Membership is open to anyone interested in APL. The membership year 
normally runs from 1st May to 30th April, but new members may join from 1st 
August, November or February if preferred. The British APL Association is a 
special interest group of the British Computer Society, Reg. Charity No. 292,786 


Name: 
Address: 


Postcode / Country: 
Telephone Number: 
Email Address: 


Category (please tick box) to run from: 1st Му Û August Û Nov Û Feb 

UK private membership ........... EY 20) 

Overseas private membership . 
Airmail supplement (not needed for Europe) 

UK Corporate membership ............ 

Corporate membership overseas . . 


f A35 
Sustaining membership ......... s s sisse. £430 
Non-voting UK member (student/OAP/ unemployed only) £6 


DODOCCODO 


PAYMENT — in Sterling or by Visa/Mastercard/JCB only 

Payment should be enclosed with membership applications in the form of a UK 
Sterling cheque to "The British APL Association", or you may quote your 
Mastercard, Visa or JCB number. 


I authorise you to debit my Visa/Mastercard/JCB account 


Number: L1 13 tirada besas Lir. Expiry date: ti iiid 


for the membership category indicated above, 


Data Protection Act 
"m 5 The information supplied may be 
a annually, at the Prevailing rate, until further notice | „== d 
one year's subscription only in accordance with the registration 
(ofthe Brits) Computer Society. 
(please tick the required option above) 
Signature: Send the completed form to: 


British APL Association, c/o Rowena Small, 8 Cardigan Road, LONDON E3 5HU, UK 
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The British APL Association is a Specialist Group of the British Computer Society. It is administered by a Committee 
of officers who are elected by a postal ballot of Association members prior to the Annual General Meeting. Working 
groups are also established in areas such as activity planning and journal production. Offers of assistance and 
involvement with any Association matters аге welcomed and should be addressed in the first instance to the Secretary. 


1995/96 Committee 
Chairman: Dr Alan Mayer European Business Management School, 
01792-205678x4274 ‘Swansea University, 
a.d.mayer @swansea.ac.uk Singleton Park, SWANSEA SA2 8PP 
Secretary: Sylvia Camacho 11 Auburn Road, Redland, 
0117-973 0036 BRISTOL, BS6 61.5 
100612.1057@compuserve.com 
Treasurer: Nicholas Small 8 Cardigan Road, 
0181-980 7870 LONDON E35HU 
Journal Editor: Anthony Camacho 11 Auburn Road, Redland, 
0117-973 0036 BRISTOL, 
acamachocixcompulinkcouk ^ BS66LS 
Activities: Vacant Post 
Education: Dr lan Clark 9 Hill End, Frosterley 
01388-527190 Bishop Auckland 
100021.3073 @compuserve.com Со. Durham DL132SX 
Technical: Vacant Post 
Publicity: David Eastwood MicroAPL Ltd., 
0171-922 8866 South Bank Te 
MicroAPL @microapl.demon.co.uk 90 London Road, LONDON SEI 6LN 
Recruitment: Jon Sandles 138 Burton Stone Lane, 
01904-612882 York YO3 6DF 
100257.1756 6 compuserve.com 
Administration: Rowena Small 8 Cardigan Road, 
0181-980 7870 LONDON E3 SHU 
Journal Working Group 
Editor: Anthony Camacho 0117-973 0036 
Production: Adrian & Gill Smith Brook House, Gilling East, YORK (01439-788385) 
Advertising: Gill Smith Brook House, Gilling East, YORK (01439-788385) 


Support Team: Jonathan Barman (01488-648575), Duncan Pearson (01653-618900), 
Richard and Adam Weber (01302-539761), Sylvia Camacho, Ray Cannon (01252-874697), 
John Searle (0181-858 6811), David Ziemann (0181-348 4039), Jon Sandles 


‘Typeset by APL-385 with MS Word 5.0 and GoScript 
Printed in England by Short-Run Press Lid, Exeter 
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VECTOR is the quarterly Journal of the British APL Association and is distributed to Association 
members in the UK and overseas. The British APL Association is a Specialist Group of the British 
Computer Society. APL stands for "A Programming Language” — an interactive computer 
language noted for its elegance, conciseness and fast development speed. It is supported on 
most mainframes, workstations and personal computers. 


SUSTAINING MEMBERS 


The Committee of the British APL Association wish to acknowledge the generous financial 
support of the following Association Sustaining Members. In many cases these organisations also 
provide manpower and administrative assistance to the Association at their own cost. 


Causeway Graphical Systems Ltd 
5 The Maltings, Castlegate, 

MALTON, North Yorks YO17 ODP 

Tel: 01653-696760 

Fax: 01653-697719 

Email: 100265. 1564@ compuserve.com 


Dyadic Systems Ltd 
Riverside View, Basing Road, 
Old Basing, BASINGSTOKE, 
Hants, RG24 OAL 
Tel:01256-811125 
Fax:01256-811130 
Email:Sales@dyadic.com 


Insight Systems ApS 
Nordre Strandvej 119A 
DK-3150 Hellebæk 
Denmark 

Tel: +45 42 10 70 22 

Fax: +45 42 10 75 74 
Email: insight & ineLuni-c.dk 


MicroAPL Ltd 

South Bank Technopark 

90 London Road 

LONDON SE1 6LN 

Tel:0171-922 8866 

Fax:0171-928 1006 

Email:microap!@ microapl.demon.co.uk 


Dutch APL Association 
Postbus 1341 

3430BH Nieuwegein 
Netherlands 
Те:03474-2337 


Compass R&D Ltd 
10 Frederick Sanger Road 
Surrey Research Park. 
GUILDFORD, Surrey GU2 SYD 
Tel.01483-302249 
Ғах:01483-302279 


HMW Trading Systems Ltd 

Hamilton House, 

1 Temple Avenue, 

LONDON EC4Y OHA 

Tel:0171-353 8900 

Ғах:0171-353 3325 
Email:100020.2632@ Compuserve.com 


Manugistics 

2115 East Jefferson St 
Rockville 

MARYLAND 20852 USA 
Tel: +1 (301) 984-5412 


Soliton Associates Lid 
Groot Blankenberg 53 
1082 AC Amsterdam 
Netherlands 

Tel:+31 20 646 4475 
Fax:+31 20 644 1206 
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